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

Disrupted tRNA Genes and tRNA Fragments: A Perspective on tRNA Gene Evolution

Transfer RNAs (tRNAs) are small non-coding RNAs with lengths of approximately 70-100 nt. They are directly involved in protein synthesis by carrying amino acids to the ribosome. In this sense, tRNAs are key molecules that connect the RNA world and the protein world. Thus, study of the evolution of tRNA molecules may reveal the processes that led to the establishment of the central dogma: genetic information flows from DNA to RNA to protein. Thanks to the development of DNA sequencers in this century, we have determined a huge number of nucleotide sequences from complete genomes as well as from transcriptomes in many species. Recent analyses of these large data sets have shown that particular tRNA genes, especially in Archaea, are disrupted in unique ways: some tRNA genes contain multiple introns and some are split genes. Even tRNA molecules themselves are fragmented post-transcriptionally in many species. These fragmented small RNAs are known as tRNA-derived fragments (tRFs). In this review, I summarize the progress of research into the disrupted tRNA genes and the tRFs, and propose a possible model for the molecular evolution of tRNAs based on the concept of the combination of fragmented tRNA halves.

Cultivation of a human-associated TM7 phylotype reveals a reduced genome and epibiotic parasitic lifestyle

The candidate phylum TM7 is globally distributed and often associated with human inflammatory mucosal diseases. Despite its prevalence, the TM7 phylum remains recalcitrant to cultivation, making it one of the most enigmatic phyla known. In this study, we cultivated a TM7 phylotype (TM7x) from the human oral cavity. This extremely small coccus (200-300 nm) has a distinctive lifestyle not previously observed in human-associated microbes. It is an obligate epibiont of an Actinomyces odontolyticus strain (XH001) yet also has a parasitic phase, thereby killing its host. This first completed genome (705 kb) for a human-associated TM7 phylotype revealed a complete lack of amino acid biosynthetic capacity. Comparative genomics analyses with uncultivated environmental TM7 assemblies show remarkable conserved gene synteny and only minimal gene loss/gain that may have occurred as TM7x adapted to conditions within the human host. Transcriptomic and metabolomic profiles provided the first indications, to our knowledge, that there is signaling interaction between TM7x and XH001. Furthermore, the induction of TNF-α production in macrophages by XH001 was repressed in the presence of TM7x, suggesting its potential immune suppression ability. Overall, our data provide intriguing insights into the uncultivability, pathogenicity, and unique lifestyle of this previously uncharacterized oral TM7 phylotype.

Rooting the domain archaea by phylogenomic analysis supports the foundation of the new kingdom Proteoarchaeota

The first 16S rRNA-based phylogenies of the Archaea showed a deep division between two groups, the kingdoms Euryarchaeota and Crenarchaeota. This bipartite classification has been challenged by the recent discovery of new deeply branching lineages (e.g., Thaumarchaeota, Aigarchaeota, Nanoarchaeota, Korarchaeota, Parvarchaeota, Aenigmarchaeota, Diapherotrites, and Nanohaloarchaeota) which have also been given the same taxonomic status of kingdoms. However, the phylogenetic position of some of these lineages is controversial. In addition, phylogenetic analyses of the Archaea have often been carried out without outgroup sequences, making it difficult to determine if these taxa actually define lineages at the same level as the Euryarchaeota and Crenarchaeota. We have addressed the question of the position of the root of the Archaea by reconstructing rooted archaeal phylogenetic trees using bacterial sequences as outgroup. These trees were based on commonly used conserved protein markers (32 ribosomal proteins) as well as on 38 new markers identified through phylogenomic analysis. We thus gathered a total of 70 conserved markers that we analyzed as a concatenated data set. In contrast with previous analyses, our trees consistently placed the root of the archaeal tree between the Euryarchaeota (including the Nanoarchaeota and other fast-evolving lineages) and the rest of archaeal species, which we propose to class within the new kingdom Proteoarchaeota. This implies the relegation of several groups previously classified as kingdoms (e.g., Crenarchaeota, Thaumarchaeota, Aigarchaeota, and Korarchaeota) to a lower taxonomic rank. In addition to taxonomic implications, this profound reorganization of the archaeal phylogeny has also consequences on our appraisal of the nature of the last archaeal ancestor, which most likely was a complex organism with a gene-rich genome.

Comparative genomics highlights the unique biology of Methanomassiliicoccales, a Thermoplasmatales-related seventh order of methanogenic archaea that encodes pyrrolysine

BACKGROUND

A seventh order of methanogens, the Methanomassiliicoccales, has been identified in diverse anaerobic environments including the gastrointestinal tracts (GIT) of humans and other animals and may contribute significantly to methane emission and global warming. Methanomassiliicoccales are phylogenetically distant from all other orders of methanogens and belong to a large evolutionary branch composed by lineages of non-methanogenic archaea such as Thermoplasmatales, the Deep Hydrothermal Vent Euryarchaeota-2 (DHVE-2, Aciduliprofundum boonei) and the Marine Group-II (MG-II). To better understand this new order and its relationship to other archaea, we manually curated and extensively compared the genome sequences of three Methanomassiliicoccales representatives derived from human GIT microbiota, "Candidatus Methanomethylophilus alvus", "Candidatus Methanomassiliicoccus intestinalis" and Methanomassiliicoccus luminyensis.

RESULTS

Comparative analyses revealed atypical features, such as the scattering of the ribosomal RNA genes in the genome and the absence of eukaryotic-like histone gene otherwise present in most of Euryarchaeota genomes. Previously identified in Thermoplasmatales genomes, these features are presently extended to several completely sequenced genomes of this large evolutionary branch, including MG-II and DHVE2. The three Methanomassiliicoccales genomes share a unique composition of genes involved in energy conservation suggesting an original combination of two main energy conservation processes previously described in other methanogens. They also display substantial differences with each other, such as their codon usage, the nature and origin of their CRISPRs systems and the genes possibly involved in particular environmental adaptations. The genome of M. luminyensis encodes several features to thrive in soil and sediment conditions suggesting its larger environmental distribution than GIT. Conversely, "Ca. M. alvus" and "Ca. M. intestinalis" do not present these features and could be more restricted and specialized on GIT. Prediction of the amber codon usage, either as a termination signal of translation or coding for pyrrolysine revealed contrasted patterns among the three genomes and suggests a different handling of the Pyl-encoding capacity.

CONCLUSIONS

This study represents the first insights into the genomic organization and metabolic traits of the seventh order of methanogens. It suggests contrasted evolutionary history among the three analyzed Methanomassiliicoccales representatives and provides information on conserved characteristics among the overall methanogens and among Thermoplasmata.

Insights into the metabolism, lifestyle and putative evolutionary history of the novel archaeal phylum 'Diapherotrites'

The archaeal phylum 'Diapherotrites' was recently proposed based on phylogenomic analysis of genomes recovered from an underground water seep in an abandoned gold mine (Homestake mine in Lead, SD, USA). Here we present a detailed analysis of the metabolic capabilities and genomic features of three single amplified genomes (SAGs) belonging to the 'Diapherotrites'. The most complete of the SAGs, Candidatus 'Iainarchaeum andersonii' (Cand. IA), had a small genome (∼1.24 Mb), short average gene length (822 bp), one ribosomal RNA operon, high coding density (∼90.4%), high percentage of overlapping genes (27.6%) and low incidence of gene duplication (2.16%). Cand. IA genome possesses limited catabolic capacities that, nevertheless, could theoretically support a free-living lifestyle by channeling a narrow range of substrates such as ribose, polyhydroxybutyrate and several amino acids to acetyl-coenzyme A. On the other hand, Cand. IA possesses relatively well-developed anabolic capabilities, although it remains auxotrophic for several amino acids and cofactors. Phylogenetic analysis suggests that the majority of Cand. IA anabolic genes were acquired from bacterial donors via horizontal gene transfer. We thus propose that members of the 'Diapherotrites' have evolved from an obligate symbiotic ancestor by acquiring anabolic genes from bacteria that enabled independent biosynthesis of biological molecules previously acquired from symbiotic hosts. 'Diapherotrites' 16S rRNA genes exhibit multiple mismatches with the majority of archaeal 16S rRNA primers, a fact that could be responsible for their observed rarity in amplicon-generated data sets. The limited substrate range, complex growth requirements and slow growth rate predicted could be responsible for its refraction to isolation.

Life on the edge: functional genomic response of Ignicoccus hospitalis to the presence of Nanoarchaeum equitans

The marine hyperthermophilic crenarchaeon Ignicoccus hospitalis supports the propagation on its surface of Nanoarchaeum equitans, an evolutionarily enigmatic archaeon that resembles highly derived parasitic and symbiotic bacteria. The cellular and molecular mechanisms that enable this interarchaea relationship and the intimate physiologic consequences to I. hospitalis are unknown. Here, we used concerted proteomic and transcriptomic analyses to probe into the functional genomic response of I. hospitalis as N. equitans multiplies on its surface. The expression of over 97% of the genes was detected at mRNA level and over 80% of the predicted proteins were identified and their relative abundance measured by proteomics. These indicate that little, if any, of the host genomic information is silenced during growth in the laboratory. The primary response to N. equitans was at the membrane level, with increases in relative abundance of most protein complexes involved in energy generation as well as that of several transporters and proteins involved in cellular membrane stabilization. Similar upregulation was observed for genes and proteins involved in key metabolic steps controlling nitrogen and carbon metabolism, although the overall biosynthetic pathways were marginally impacted. Proliferation of N. equitans resulted, however, in selective downregulation of genes coding for transcription factors and replication and cell cycle control proteins as I. hospitalis shifted its physiology from its own cellular growth to that of its ectosymbiont/parasite. The combination of these multiomic approaches provided an unprecedented level of detail regarding the dynamics of this interspecies interaction, which is especially pertinent as these organisms are not genetically tractable.

Genome-wide miRNA seeds prediction in Archaea

Growing evidence indicates that miRNA genes exist in the archaeal genome, though the functional role of such noncoding RNA remains unclear. Here, we integrated the phylogenetic information of available archaeal genomes to predict miRNA seeds (typically defined as the 2-8 nucleotides of mature miRNAs) on the genomic scale. Finally, we found 2649 candidate seeds with significant conservation signal. Eleven of 29 unique seeds from previous study support our result (P value <0.01), which demonstrates that the pipeline is suitable to predict experimentally detectable miRNA seeds. The statistical significance of the overlap between the detected archaeal seeds and known eukaryotic seeds shows that the miRNA may evolve before the divergence of these two domains of cellular life. In addition, miRNA targets are enriched for genes involved in transcriptional regulation, which is consistent with the situation in eukaryote. Our research will enhance the regulatory network analysis in Archaea.

Protein Ser/Thr/Tyr phosphorylation in the Archaea

The third domain of life, the Archaea (formerly Archaebacteria), is populated by a physiologically diverse set of microorganisms, many of which reside at the ecological extremes of our global environment. Although ostensibly prokaryotic in morphology, the Archaea share much closer evolutionary ties with the Eukarya than with the superficially more similar Bacteria. Initial genomic, proteomic, and biochemical analyses have revealed the presence of "eukaryotic" protein kinases and phosphatases and an intriguing set of serine-, threonine-, and tyrosine-phosphorylated proteins in the Archaea that may offer new insights into this important regulatory mechanism.

Identification and codon reading properties of 5-cyanomethyl uridine, a new modified nucleoside found in the anticodon wobble position of mutant haloarchaeal isoleucine tRNAs

Most archaea and bacteria use a modified C in the anticodon wobble position of isoleucine tRNA to base pair with A but not with G of the mRNA. This allows the tRNA to read the isoleucine codon AUA without also reading the methionine codon AUG. To understand why a modified C, and not U or modified U, is used to base pair with A, we mutated the C34 in the anticodon of Haloarcula marismortui isoleucine tRNA (tRNA2(Ile)) to U, expressed the mutant tRNA in Haloferax volcanii, and purified and analyzed the tRNA. Ribosome binding experiments show that although the wild-type tRNA2(Ile) binds exclusively to the isoleucine codon AUA, the mutant tRNA binds not only to AUA but also to AUU, another isoleucine codon, and to AUG, a methionine codon. The G34 to U mutant in the anticodon of another H. marismortui isoleucine tRNA species showed similar codon binding properties. Binding of the mutant tRNA to AUG could lead to misreading of the AUG codon and insertion of isoleucine in place of methionine. This result would explain why most archaea and bacteria do not normally use U or a modified U in the anticodon wobble position of isoleucine tRNA for reading the codon AUA. Biochemical and mass spectrometric analyses of the mutant tRNAs have led to the discovery of a new modified nucleoside, 5-cyanomethyl U in the anticodon wobble position of the mutant tRNAs. 5-Cyanomethyl U is present in total tRNAs from euryarchaea but not in crenarchaea, eubacteria, or eukaryotes.

Functional implications of ribosomal RNA methylation in response to environmental stress

The study of post-transcriptional RNA modifications has long been focused on the roles these chemical modifications play in maintaining ribosomal function. The field of ribosomal RNA modification has reached a milestone in recent years with the confirmation of the final unknown ribosomal RNA methyltransferase in Escherichia coli in 2012. Furthermore, the last 10 years have brought numerous discoveries in non-coding RNAs and the roles that post-transcriptional modification play in their functions. These observations indicate the need for a revitalization of this field of research to understand the role modifications play in maintaining cellular health in a dynamic environment. With the advent of high-throughput sequencing technologies, the time is ripe for leaps and bounds forward. This review discusses ribosomal RNA methyltransferases and their role in responding to external stress in Escherichia coli, with a specific focus on knockout studies and on analysis of transcriptome data with respect to rRNA methyltransferases.

Close encounters of the third domain: the emerging genomic view of archaeal diversity and evolution

The Archaea represent the so-called Third Domain of life, which has evolved in parallel with the Bacteria and which is implicated to have played a pivotal role in the emergence of the eukaryotic domain of life. Recent progress in genomic sequencing technologies and cultivation-independent methods has started to unearth a plethora of data of novel, uncultivated archaeal lineages. Here, we review how the availability of such genomic data has revealed several important insights into the diversity, ecological relevance, metabolic capacity, and the origin and evolution of the archaeal domain of life.

Predominant Acidilobus-like populations from geothermal environments in yellowstone national park exhibit similar metabolic potential in different hypoxic microbial communities

High-temperature (>70°C) ecosystems in Yellowstone National Park (YNP) provide an unparalleled opportunity to study chemotrophic archaea and their role in microbial community structure and function under highly constrained geochemical conditions. Acidilobus spp. (order Desulfurococcales) comprise one of the dominant phylotypes in hypoxic geothermal sulfur sediment and Fe(III)-oxide environments along with members of the Thermoproteales and Sulfolobales. Consequently, the primary goals of the current study were to analyze and compare replicate de novo sequence assemblies of Acidilobus-like populations from four different mildly acidic (pH 3.3 to 6.1) high-temperature (72°C to 82°C) environments and to identify metabolic pathways and/or protein-encoding genes that provide a detailed foundation of the potential functional role of these populations in situ. De novo assemblies of the highly similar Acidilobus-like populations (>99% 16S rRNA gene identity) represent near-complete consensus genomes based on an inventory of single-copy genes, deduced metabolic potential, and assembly statistics generated across sites. Functional analysis of coding sequences and confirmation of gene transcription by Acidilobus-like populations provide evidence that they are primarily chemoorganoheterotrophs, generating acetyl coenzyme A (acetyl-CoA) via the degradation of carbohydrates, lipids, and proteins, and auxotrophic with respect to several external vitamins, cofactors, and metabolites. No obvious pathways or protein-encoding genes responsible for the dissimilatory reduction of sulfur were identified. The presence of a formate dehydrogenase (Fdh) and other protein-encoding genes involved in mixed-acid fermentation supports the hypothesis that Acidilobus spp. function as degraders of complex organic constituents in high-temperature, mildly acidic, hypoxic geothermal systems.

Microbial syntrophy: interaction for the common good

Classical definitions of syntrophy focus on a process, performed through metabolic interaction between dependent microbial partners, such as the degradation of complex organic compounds under anoxic conditions. However, examples from past and current scientific discoveries suggest that a new, simple but wider definition is necessary to cover all aspects of microbial syntrophy. We suggest the term 'obligately mutualistic metabolism', which still focuses on microbial metabolic cooperation but also includes an ecological aspect: the benefit for both partners. By the combined metabolic activity of microorganisms, endergonic reactions can become exergonic through the efficient removal of products and therefore enable a microbial community to survive with minimal energy resources. Here, we explain the principles of classical and non-classical syntrophy and illustrate the concepts with various examples. We present biochemical fundamentals that allow microorganism to survive under a range of environmental conditions and to drive important biogeochemical processes. Novel technologies have contributed to the understanding of syntrophic relationships in cultured and uncultured systems. Recent research highlights that obligately mutualistic metabolism is not limited to certain metabolic pathways nor to certain environments or microorganisms. This beneficial microbial interaction is not restricted to the transfer of reducing agents such as hydrogen or formate, but can also involve the exchange of organic, sulfurous- and nitrogenous compounds or the removal of toxic compounds.

The impact of genomics on research in diversity and evolution of archaea

Since the definition of archaea as a separate domain of life along with bacteria and eukaryotes, they have become one of the most interesting objects of modern microbiology, molecular biology, and biochemistry. Sequencing and analysis of archaeal genomes were especially important for studies on archaea because of a limited availability of genetic tools for the majority of these microorganisms and problems associated with their cultivation. Fifteen years since the publication of the first genome of an archaeon, more than one hundred complete genome sequences of representatives of different phylogenetic groups have been determined. Analysis of these genomes has expanded our knowledge of biology of archaea, their diversity and evolution, and allowed identification and characterization of new deep phylogenetic lineages of archaea. The development of genome technologies has allowed sequencing the genomes of uncultivated archaea directly from enrichment cultures, metagenomic samples, and even from single cells. Insights have been gained into the evolution of key biochemical processes in archaea, such as cell division and DNA replication, the role of horizontal gene transfer in the evolution of archaea, and new relationships between archaea and eukaryotes have been revealed.

Yellowstone lake nanoarchaeota

Considerable Nanoarchaeota novelty and diversity were encountered in Yellowstone Lake, Yellowstone National Park (YNP), where sampling targeted lake floor hydrothermal vent fluids, streamers and sediments associated with these vents, and in planktonic photic zones in three different regions of the lake. Significant homonucleotide repeats (HR) were observed in pyrosequence reads and in near full-length Sanger sequences, averaging 112 HR per 1349 bp clone and could confound diversity estimates derived from pyrosequencing, resulting in false nucleotide insertions or deletions (indels). However, Sanger sequencing of two different sets of PCR clones (110 bp, 1349 bp) demonstrated that at least some of these indels are real. The majority of the Nanoarchaeota PCR amplicons were vent associated; however, curiously, one relatively small Nanoarchaeota OTU (71 pyrosequencing reads) was only found in photic zone water samples obtained from a region of the lake furthest removed from the hydrothermal regions of the lake. Extensive pyrosequencing failed to demonstrate the presence of an Ignicoccus lineage in this lake, suggesting the Nanoarchaeota in this environment are associated with novel Archaea hosts. Defined phylogroups based on near full-length PCR clones document the significant Nanoarchaeota 16S rRNA gene diversity in this lake and firmly establish a terrestrial clade distinct from the marine Nanoarcheota as well as from other geographical locations.

Synthesis and dissolution of hemicatenanes by type IA DNA topoisomerases

Type IA DNA topoisomerases work with a unique mechanism of strand passage through an enzyme-bridged, ssDNA gate, thus enabling them to carry out diverse reactions in processing structures important for replication, recombination, and repair. Here we report a unique reaction mediated by an archaeal type IA topoisomerase, the synthesis and dissolution of hemicatenanes. We cloned, purified, and characterized an unusual type IA enzyme from a hyperthermophilic archaeum, Nanoarchaeum equitans, which is split into two pieces. The recombinant heterodimeric enzyme has the expected activities in its preference of relaxing negatively supercoiled DNA. Its amino acid sequence and cleavage site sequence analysis suggest that it is topoisomerase III, and therefore we named it "NeqTop3." At high enzyme concentrations, NeqTop3 can generate high-molecular-weight DNA networks. Biochemical and electron microscopic data indicate that the DNA networks are connected through hemicatenane linkages. The hemicatenane formation likely is mediated by the single-strand passage through denatured bubbles in the substrate DNA under high temperature. NeqTop3 at lower concentrations can reverse hemicatenanes. A complex of human topoisomerase 3α, Bloom helicase, and RecQ-mediated genome instability protein 1 and 2 can partially disentangle the hemicatenane network. Both the formation and dissolution of hemicatenanes by type IA topoisomerases demonstrate that these enzymes have an important role in regulating intermediates from replication, recombination, and repair.

A brief history of the discovery of hyperthermophilic life

Hyperthermophiles, growing optimally at 80°C and above were first discovered in 1981. They represent the upper temperature border of life and are found within water-containing terrestrial and submarine environments of active volcanism and geothermally heated subterranean rocks. The energy-yielding reactions represent mainly anaerobic and aerobic types of respiration rather than fermentation. Within the ss (single-stranded) rRNA phylogenetic tree, hyperthermophiles occupy all of the short deep branches closest to the root. Members of the deepest branch-offs are represented by the newly found Nanoarchaeota and Korarchaeota.

C/D box sRNA, CRISPR RNA and tRNA processing in an archaeon with a minimal fragmented genome

The analysis of deep sequencing data allows for a genome-wide overview of all the small RNA molecules (the 'sRNome') that are present in a single organism. In the present paper, we review the processing of CRISPR (clustered regularly interspaced short palindromic repeats) RNA, C/D box sRNA (small non-coding RNA) and tRNA in Nanoarchaeum equitans. The minimal and fragmented genome of this tiny archaeon permits a sequencing depth that enables the identification of processing intermediates in the study of RNA processing pathways. These intermediates include circular C/D box sRNA molecules and tRNA half precursors.

Increments and duplication events of enzymes and transcription factors influence metabolic and regulatory diversity in prokaryotes

In this work, the content of enzymes and DNA-binding transcription factors (TFs) in 794 non-redundant prokaryotic genomes was evaluated. The identification of enzymes was based on annotations deposited in the KEGG database as well as in databases of functional domains (COG and PFAM) and structural domains (Superfamily). For identifications of the TFs, hidden Markov profiles were constructed based on well-known transcriptional regulatory families. From these analyses, we obtained diverse and interesting results, such as the negative rate of incremental changes in the number of detected enzymes with respect to the genome size. On the contrary, for TFs the rate incremented as the complexity of genome increased. This inverse related performance shapes the diversity of metabolic and regulatory networks and impacts the availability of enzymes and TFs. Furthermore, the intersection of the derivatives between enzymes and TFs was identified at 9,659 genes, after this point, the regulatory complexity grows faster than metabolic complexity. In addition, TFs have a low number of duplications, in contrast to the apparent high number of duplications associated with enzymes. Despite the greater number of duplicated enzymes versus TFs, the increment by which duplicates appear is higher in TFs. A lower proportion of enzymes among archaeal genomes (22%) than in the bacterial ones (27%) was also found. This low proportion might be compensated by the interconnection between the metabolic pathways in Archaea. A similar proportion was also found for the archaeal TFs, for which the formation of regulatory complexes has been proposed. Finally, an enrichment of multifunctional enzymes in Bacteria, as a mechanism of ecological adaptation, was detected.

Domain structures and inter-domain interactions defining the holoenzyme architecture of archaeal d-family DNA polymerase

Archaea-specific D-family DNA polymerase (PolD) forms a dimeric heterodimer consisting of two large polymerase subunits and two small exonuclease subunits. According to the protein-protein interactions identified among the domains of large and small subunits of PolD, a symmetrical model for the domain topology of the PolD holoenzyme is proposed. The experimental evidence supports various aspects of the model. The conserved amphipathic nature of the N-terminal putative α-helix of the large subunit plays a key role in the homodimeric assembly and the self-cyclization of the large subunit and is deeply involved in the archaeal PolD stability and activity. We also discuss the evolutional transformation from archaeal D-family to eukaryotic B-family polymerase on the basis of the structural information.

Phylogenetic and Functional Analysis of Metagenome Sequence from High-Temperature Archaeal Habitats Demonstrate Linkages between Metabolic Potential and Geochemistry

Geothermal habitats in Yellowstone National Park (YNP) provide an unparalleled opportunity to understand the environmental factors that control the distribution of archaea in thermal habitats. Here we describe, analyze, and synthesize metagenomic and geochemical data collected from seven high-temperature sites that contain microbial communities dominated by archaea relative to bacteria. The specific objectives of the study were to use metagenome sequencing to determine the structure and functional capacity of thermophilic archaeal-dominated microbial communities across a pH range from 2.5 to 6.4 and to discuss specific examples where the metabolic potential correlated with measured environmental parameters and geochemical processes occurring in situ. Random shotgun metagenome sequence (∼40–45 Mb Sanger sequencing per site) was obtained from environmental DNA extracted from high-temperature sediments and/or microbial mats and subjected to numerous phylogenetic and functional analyses. Analysis of individual sequences (e.g., MEGAN and G + C content) and assemblies from each habitat type revealed the presence of dominant archaeal populations in all environments, 10 of whose genomes were largely reconstructed from the sequence data. Analysis of protein family occurrence, particularly of those involved in energy conservation, electron transport, and autotrophic metabolism, revealed significant differences in metabolic strategies across sites consistent with differences in major geochemical attributes (e.g., sulfide, oxygen, pH). These observations provide an ecological basis for understanding the distribution of indigenous archaeal lineages across high-temperature systems of YNP.

Insights into archaeal evolution and symbiosis from the genomes of a nanoarchaeon and its inferred crenarchaeal host from Obsidian Pool, Yellowstone National Park

Background

A single cultured marine organism, Nanoarchaeum equitans, represents the Nanoarchaeota branch of symbiotic Archaea, with a highly reduced genome and unusual features such as multiple split genes.

Results

The first terrestrial hyperthermophilic member of the Nanoarchaeota was collected from Obsidian Pool, a thermal feature in Yellowstone National Park, separated by single cell isolation, and sequenced together with its putative host, a Sulfolobales archaeon. Both the new Nanoarchaeota (Nst1) and N. equitans lack most biosynthetic capabilities, and phylogenetic analysis of ribosomal RNA and protein sequences indicates that the two form a deep-branching archaeal lineage. However, the Nst1 genome is more than 20% larger, and encodes a complete gluconeogenesis pathway as well as the full complement of archaeal flagellum proteins. With a larger genome, a smaller repertoire of split protein encoding genes and no split non-contiguous tRNAs, Nst1 appears to have experienced less severe genome reduction than N. equitans. These findings imply that, rather than representing ancestral characters, the extremely compact genomes and multiple split genes of Nanoarchaeota are derived characters associated with their symbiotic or parasitic lifestyle. The inferred host of Nst1 is potentially autotrophic, with a streamlined genome and simplified central and energetic metabolism as compared to other Sulfolobales.

Conclusions

Comparison of the N. equitans and Nst1 genomes suggests that the marine and terrestrial lineages of Nanoarchaeota share a common ancestor that was already a symbiont of another archaeon. The two distinct Nanoarchaeota-host genomic data sets offer novel insights into the evolution of archaeal symbiosis and parasitism, enabling further studies of the cellular and molecular mechanisms of these relationships.

Reviewers

This article was reviewed by Patrick Forterre, Bettina Siebers (nominated by Michael Galperin) and Purification Lopez-Garcia

Ribonucleoproteins in archaeal pre-rRNA processing and modification

Given that ribosomes are one of the most important cellular macromolecular machines, it is not surprising that there is intensive research in ribosome biogenesis. Ribosome biogenesis is a complex process. The maturation of ribosomal RNAs (rRNAs) requires not only the precise cleaving and folding of the pre-rRNA but also extensive nucleotide modifications. At the heart of the processing and modifications of pre-rRNAs in Archaea and Eukarya are ribonucleoprotein (RNP) machines. They are called small RNPs (sRNPs), in Archaea, and small nucleolar RNPs (snoRNPs), in Eukarya. Studies on ribosome biogenesis originally focused on eukaryotic systems. However, recent studies on archaeal sRNPs have provided important insights into the functions of these RNPs. This paper will introduce archaeal rRNA gene organization and pre-rRNA processing, with a particular focus on the discovery of the archaeal sRNP components, their functions in nucleotide modification, and their structures.

Detection of genomic idiosyncrasies using fuzzy phylogenetic profiles

Phylogenetic profiles express the presence or absence of genes and their homologs across a number of reference genomes. They have emerged as an elegant representation framework for comparative genomics and have been used for the genome-wide inference and discovery of functionally linked genes or metabolic pathways. As the number of reference genomes grows, there is an acute need for faster and more accurate methods for phylogenetic profile analysis with increased performance in speed and quality. We propose a novel, efficient method for the detection of genomic idiosyncrasies, i.e. sets of genes found in a specific genome with peculiar phylogenetic properties, such as intra-genome correlations or inter-genome relationships. Our algorithm is a four-step process where genome profiles are first defined as fuzzy vectors, then discretized to binary vectors, followed by a de-noising step, and finally a comparison step to generate intra- and inter-genome distances for each gene profile. The method is validated with a carefully selected benchmark set of five reference genomes, using a range of approaches regarding similarity metrics and pre-processing stages for noise reduction. We demonstrate that the fuzzy profile method consistently identifies the actual phylogenetic relationship and origin of the genes under consideration for the majority of the cases, while the detected outliers are found to be particular genes with peculiar phylogenetic patterns. The proposed method provides a time-efficient and highly scalable approach for phylogenetic stratification, with the detected groups of genes being either similar to their own genome profile or different from it, thus revealing atypical evolutionary histories.

Multifunctional enzymes in archaea: promiscuity and moonlight

Enzymes from many archaea colonizing extreme environments are of great interest because of their potential for various biotechnological processes and scientific value of evolution. Many enzymes from archaea have been reported to catalyze promiscuous reactions or moonlight in different functions. Here, we summarize known archaeal enzymes of both groups that include different kinds of proteins. Knowledge of their biochemical properties and three-dimensional structures has proved invaluable in understanding mechanism, application, and evolutionary implications of this manifestation. In addition, the review also summarizes the methods to unravel the extra function which almost was discovered serendipitously. The study of these amazing enzymes will provide clues to optimize protein engineering applications and how enzymes might have evolved on Earth.

The evolution, function, structure, and expression of the plant sHSPs

Small heat shock proteins are a diverse, ancient, and important family of proteins. All organisms possess small heat shock proteins (sHSPs), indicating that these proteins evolved very early in the history of life prior to the divergence of the three domains of life (Archaea, Bacteria, and Eukarya). Comparing the structures of sHSPs from diverse organisms across these three domains reveals that despite considerable amino acid divergence, many structural features are conserved. Comparisons of the sHSPs from diverse organisms reveal conserved structural features including an oligomeric form with a β-sandwich that forms a hollow ball. This conservation occurs despite significant divergence in primary sequences. It is well established that sHSPs are molecular chaperones that prevent misfolding and irreversible aggregation of their client proteins. Most notably, the sHSPs are extremely diverse and variable in plants. Some plants have >30 individual sHSPs. Land plants, unlike other groups, possess distinct sHSP subfamilies. Most are highly up-regulated in response to heat and other stressors. Others are selectively expressed in seeds and pollen, and a few are constitutively expressed. As a family, sHSPs have a clear role in thermotolerance, but attributing specific effects to individual proteins has proved challenging. Considerable progress has been made during the last 15 years in understanding the sHSPs. However, answers to many important questions remain elusive, suggesting that the next 15 years will be at least equally rewarding.

Archaea in symbioses

During the last few years, the analysis of microbial diversity in various habitats greatly increased our knowledge on the kingdom Archaea. At the same time, we became aware of the multiple ways in which Archaea may interact with each other and with organisms of other kingdoms. The large group of euryarchaeal methanogens and their methane oxidizing relatives, in particular, take part in essential steps of the global methane cycle. Both of these processes, which are in reverse to each other, are partially conducted in a symbiotic interaction with different partners, either ciliates and xylophagous animals or sulfate reducing bacteria. Other symbiotic interactions are mostly of unknown ecological significance but depend on highly specific mechanisms. This paper will give an overview on interactions between Archaea and other organisms and will point out the ecological relevance of these symbiotic processes, as long as these have been already recognized.

Translation in giant viruses: a unique mixture of bacterial and eukaryotic termination schemes

Mimivirus and Megavirus are the best characterized representatives of an expanding new family of giant viruses infecting Acanthamoeba. Their most distinctive features, megabase-sized genomes carried in particles of size comparable to that of small bacteria, fill the gap between the viral and cellular worlds. These giant viruses are also uniquely equipped with genes coding for central components of the translation apparatus. The presence of those genes, thought to be hallmarks of cellular organisms, revived fundamental interrogations on the evolutionary origin of these viruses and the link they might have with the emergence of eukaryotes. In this work, we focused on the Mimivirus-encoded translation termination factor gene, the detailed primary structure of which was elucidated using computational and experimental approaches. We demonstrated that the translation of this protein proceeds through two internal stop codons via two distinct recoding events: a frameshift and a readthrough, the combined occurrence of which is unique to these viruses. Unexpectedly, the viral gene carries an autoregulatory mechanism exclusively encountered in bacterial termination factors, though the viral sequence is related to the eukaryotic/archaeal class-I release factors. This finding is a hint that the virally-encoded translation functions may not be strictly redundant with the one provided by the host. Lastly, the perplexing occurrence of a bacterial-like regulatory mechanism in a eukaryotic/archaeal homologous gene is yet another oddity brought about by the study of giant viruses.

Giant viruses, such as Mimivirus and Megavirus, have huge near-micron-sized particles and possess more genes than several cellular organisms. Furthermore their genomes encode functions not supposed to be in a virus, such as components of the protein translation apparatus. Since Lwoff in 1957, viruses are defined as ultimate obligate intracellular parasites from their need to hijack the peptide synthesis machinery of their host to replicate. We looked at the Mimivirus and Megavirus proteins that recognize the stop codons, the translation termination factors. We found that these genes contain two internal stop codons, meaning that their translation bypasses two distinct stop codons to produce a functional translation termination factor. These types of autoregulatory mechanisms are found in bacterial termination factors, although it involves only a single internal stop codon and not two, and are absent from their eukaryotic and archaeal homologs. Despite these bacterial-like features, giant viruses' termination factors have sequences that do not resemble bacterial genes but are clearly related to the eukaryotic and archaeal termination factors. Thus, giant viruses' termination factors surprisingly combine elements from eukaryotes/archaea and bacteria.

Comparative genomics guided discovery of two missing archaeal enzyme families involved in the biosynthesis of the pterin moiety of tetrahydromethanopterin and tetrahydrofolate

C-1 carriers are essential cofactors in all domains of life, and in Archaea, these can be derivatives of tetrahydromethanopterin (H₄-MPT) or tetrahydrofolate (H₄-folate). Their synthesis requires 6-hydroxymethyl-7,8-dihydropterin diphosphate (6-HMDP) as the precursor, but the nature of pathways that lead to its formation were unknown until the recent discovery of the GTP cyclohydrolase IB/MptA family that catalyzes the first step, the conversion of GTP to dihydroneopterin 2′,3′-cyclic phosphate or 7,8-dihydroneopterin triphosphate [El Yacoubi, B.; et al. (2006) J. Biol. Chem., 281, 37586–37593 and Grochowski, L. L.; et al. (2007) Biochemistry46, 6658–6667]. Using a combination of comparative genomics analyses, heterologous complementation tests, and in vitro assays, we show that the archaeal protein families COG2098 and COG1634 specify two of the missing 6-HMDP synthesis enzymes. Members of the COG2098 family catalyze the formation of 6-hydroxymethyl-7,8-dihydropterin from 7,8-dihydroneopterin, while members of the COG1634 family catalyze the formation of 6-HMDP from 6-hydroxymethyl-7,8-dihydropterin. The discovery of these missing genes solves a long-standing mystery and provides novel examples of convergent evolutions where proteins of dissimilar architectures perform the same biochemical function.

Breaking through a phylogenetic impasse: a pair of associated archaea might have played host in the endosymbiotic origin of eukaryotes

For over a century, the origin of eukaryotes has been a topic of intense debate among scientists. Although it has become widely accepted that organelles such as the mitochondria and chloroplasts arose via endosymbiosis, the origin of the eukaryotic nucleus remains enigmatic. Numerous models for the origin of the nucleus have been proposed over the years, many of which use endosymbiosis to explain its existence. Proposals of microbes whose ancestors may have served as either a host or a guest in various endosymbiotic scenarios abound, none of which have been able to sufficiently incorporate the cell biological as well as phylogenetic data which links these organisms to the nucleus. While it is generally agreed that eukaryotic nuclei share more features in common with archaea rather than with bacteria, different studies have identified either one or the other of the two major groups of archaea as potential ancestors, leading to somewhat of a stalemate. This paper seeks to resolve this impasse by presenting evidence that not just one, but a pair of archaea might have served as host to the bacterial ancestor of the mitochondria. This pair may have consisted of ancestors of both Ignicoccus hospitalis as well as its ectosymbiont/ectoparasite ‘Nanoarchaeum equitans’.

RNase P branches out from RNP to protein: organelle-triggered diversification?

RNase P is the enzyme that removes 5' leader sequences from precursor tRNAs. Remarkably, in most organisms, RNase P is a ribonucleoprotein particle where the RNA component is responsible for catalysis. In this issue of Genes & Development, Gutmann and colleagues (pp. 1022-1027) report the first organism, Arabidopsis thaliana, to employ protein-only RNase P in both its nucleus and organelles. An intriguing possibility is that replacement of RNase P ribonucleoprotein particles (RNPs) by proteins may have been triggered by the acquisition of organelles.

RNA processing in the minimal organism Nanoarchaeum equitans

BACKGROUND

The minimal genome of the tiny, hyperthermophilic archaeon Nanoarchaeum equitans contains several fragmented genes and revealed unusual RNA processing pathways. These include the maturation of tRNA molecules via the trans-splicing of tRNA halves and genomic rearrangements to compensate for the absence of RNase P.

RESULTS

Here, the RNA processing events in the N. equitans cell are analyzed using RNA-Seq deep sequencing methodology. All tRNA half precursor and tRNA termini were determined and support the tRNA trans-splicing model. The processing of CRISPR RNAs from two CRISPR clusters was verified. Twenty-seven C/D box small RNAs (sRNAs) and a H/ACA box sRNA were identified. The C/D box sRNAs were found to flank split genes, to form dicistronic tRNA-sRNA precursors and to be encoded within the tRNAMet intron.

CONCLUSIONS

The presented data provide an overview of the production and usage of small RNAs in a cell that has to survive with a highly reduced genome. N. equitans lost many essential metabolic pathways but maintains highly active CRISPR/Cas and rRNA modification systems that appear to play an important role in genome fragmentation.

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.

The unusual cell biology of the hyperthermophilic Crenarchaeon Ignicoccus hospitalis

The Crenarchaeon Ignicoccus hospitalis is an anaerobic, obligate chemolithoautotrophic hyperthermophile, growing by reduction of elemental sulfur using molecular hydrogen as electron donor. Together with Nanoarchaeum equitans it forms a unique, archaeal biocoenosis, in which I. hospitalis serves as host for N. equitans. Both organisms can be cultivated in a stable coculture which is mandatory for N. equitans but not for I. hospitalis. This strong dependence is affirmed by the fact that N. equitans obtains its lipids and amino acids from the host. I. hospitalis cells exhibit several unique features: they can adhere to surfaces by extracellular appendages ('fibers') which are not used for motility; they use a novel CO(2) fixation pathway, the dicarboxylate/4-hydroxybutyrate pathway; and they exhibit a unique cell envelope for Archaea consisting of two membranes but lacking an S-layer. These membranes form two cell compartments, a tightly packed cytoplasm surrounded by a weakly staining intermembrane compartment (IMC) with a variable width from 20 to 1,000 nm. In this IMC, many round or elongated vesicles are found which may function as carriers of lipids or proteins out of the cytoplasm. Based on immuno-EM analyses and immuno-fluorescence experiments it was demonstrated recently that the A(1)A(O) ATP synthase, the H(2):sulfur oxidoreductase complex and the acetyl-CoA synthetase (ACS) of I. hospitalis are located in its outermost membrane. Therefore, this membrane is energized and is here renamed as "outer cellular membrane" (OCM). Among all prokaryotes possessing two membranes in their cell envelope, I. hospitalis is the first organism with an energized outermost membrane and ATP synthesis outside the cytoplasm. Since DNA and ribosomes are localized in the cytoplasm, energy conservation is separated from information processing and protein biosynthesis in I. hospitalis. This raises questions concerning the function and characterization of the two membranes, the two cell compartments and of a possible ATP transfer to N. equitans.

Archaeal chromatin proteins

Archaea, along with Bacteria and Eukarya, are the three domains of life. In all living cells, chromatin proteins serve a crucial role in maintaining the integrity of the structure and function of the genome. An array of small, abundant and basic DNA-binding proteins, considered candidates for chromatin proteins, has been isolated from the Euryarchaeota and the Crenarchaeota, the two major phyla in Archaea. While most euryarchaea encode proteins resembling eukaryotic histones, crenarchaea appear to synthesize a number of unique DNA-binding proteins likely involved in chromosomal organization. Several of these proteins (e.g., archaeal histones, Sac10b homologs, Sul7d, Cren7, CC1, etc.) have been extensively studied. However, whether they are chromatin proteins and how they function in vivo remain to be fully understood. Future investigation of archaeal chromatin proteins will lead to a better understanding of chromosomal organization and gene expression in Archaea and provide valuable information on the evolution of DNA packaging in cellular life.

Structural and functional analysis of the archaeal endonuclease Nob1

Eukaryotic ribosome biogenesis requires the concerted action of numerous ribosome assembly factors, for most of which structural and functional information is currently lacking. Nob1, which can be identified in eukaryotes and archaea, is required for the final maturation of the small subunit ribosomal RNA in yeast by catalyzing cleavage at site D after export of the preribosomal subunit into the cytoplasm. Here, we show that this also holds true for Nob1 from the archaeon Pyrococcus horikoshii, which efficiently cleaves RNA-substrates containing the D-site of the preribosomal RNA in a manganese-dependent manner. The structure of PhNob1 solved by nuclear magnetic resonance spectroscopy revealed a PIN domain common with many nucleases and a zinc ribbon domain, which are structurally connected by a flexible linker. We show that amino acid residues required for substrate binding reside in the PIN domain whereas the zinc ribbon domain alone is sufficient to bind helix 40 of the small subunit rRNA. This suggests that the zinc ribbon domain acts as an anchor point for the protein on the nascent subunit positioning it in the proximity of the cleavage site.

Positive, Neutral and Negative Interactions in Cocultures between Pyrococcus furiosus and Different Methanogenic Archaea

The model organism Pyrococcus furiosus has recently been reported to interact with Methanopyrus kandleri in coculture, suggesting a H2 symbiosis. In the current study we further investigated this hypothesis by growing P. furiosus with four other hyperthermophilic methanogens providing evidence that the organisms did not only exert positive effects ( P. furiosus/ Methanocaldococcus villosus and P. furiosus/ Methanocaldococcus infernus) on each other, but also neutral ( P. furiosus/ Methanocaldococcus jannaschii) and even inhibitory interactions ( P. furiosus/ Methanotorris igneus) were detected suggesting interspecies relationships not only based on H2 symbiosis. Using various microscopic techniques we further analyzed the coculture with the highest positive interactions ( P. furiosus/ M. villosus) concerning its growth behavior on various surfaces, which turned out to be in stark contrast to the previous reported coculture of P. furiosus/ M. kandleri. This communication provides new insights into possible interactions of extremophilic Archaea in cocultures and again raises the question if and how hyperthermophilic Archaea communicate besides metabolic intermediates like H2.

Modeling a minimal cell

One important aim of synthetic biology is to develop a self-replicating biological system capable of performing useful tasks. A mathematical model of a synthetic organism would greatly enhance its value by providing a platform in which proposed modifications to the system could be rapidly prototyped and tested. Such a platform would allow the explicit connection of genomic sequence information to physiological predictions. As an initial step toward this aim, a minimal cell model (MCM) has been formulated. The MCM is defined as a model of a hypothetical cell with the minimum number of genes necessary to grow and divide in an optimally supportive culture environment. It is chemically detailed in terms of genes and gene products, as well as physiologically complete in terms of bacterial cell processes (e.g., DNA replication and cell division). A mathematical framework originally developed for modeling Escherichia coli has been used to build the platform MCM. A MCM with 241 product-coding genes (those which produce protein or stable RNA products) is presented. This gene set is genomically complete in that it codes for all the functions that a minimal chemoheterotrophic bacterium would require for sustained growth and division. With this model, the hypotheses behind a minimal gene set can be tested using a chemically detailed, dynamic, whole-cell modeling approach. Furthermore, the MCM can simulate the behavior of a whole cell that depends on the cell's (1) metabolic rates and chemical state, (2) genome in terms of expression of various genes, (3) environment both in terms of direct nutrient starvation and competitive inhibition leading to starvation, and (4) genomic sequence in terms of the chromosomal locations of genes.

Purine biosynthesis in archaea: variations on a theme

Background

The ability to perform de novo biosynthesis of purines is present in organisms in all three domains of life, reflecting the essentiality of these molecules to life. Although the pathway is quite similar in eukaryotes and bacteria, the archaeal pathway is more variable. A careful manual curation of genes in this pathway demonstrates the value of manual curation in archaea, even in pathways that have been well-studied in other domains.

Results

We searched the Integrated Microbial Genome system (IMG) for the 17 distinct genes involved in the 11 steps of de novo purine biosynthesis in 65 sequenced archaea, finding 738 predicted proteins with sequence similarity to known purine biosynthesis enzymes. Each sequence was manually inspected for the presence of active site residues and other residues known or suspected to be required for function.

Many apparently purine-biosynthesizing archaea lack evidence for a single enzyme, either glycinamide ribonucleotide formyltransferase or inosine monophosphate cyclohydrolase, suggesting that there are at least two more gene variants in the purine biosynthetic pathway to discover. Variations in domain arrangement of formylglycinamidine ribonucleotide synthetase and substantial problems in aminoimidazole carboxamide ribonucleotide formyltransferase and inosine monophosphate cyclohydrolase assignments were also identified.

Manual curation revealed some overly specific annotations in the IMG gene product name, with predicted proteins without essential active site residues assigned product names implying enzymatic activity (21 proteins, 2.8% of proteins inspected) or Enzyme Commission (E. C.) numbers (57 proteins, 7.7%). There were also 57 proteins (7.7%) assigned overly generic names and 78 proteins (10.6%) without E.C. numbers as part of the assigned name when a specific enzyme name and E. C. number were well-justified.

Conclusions

The patchy distribution of purine biosynthetic genes in archaea is consistent with a pathway that has been shaped by horizontal gene transfer, duplication, and gene loss. Our results indicate that manual curation can improve upon automated annotation for a small number of automatically-annotated proteins and can reveal a need to identify further pathway components even in well-studied pathways.

Reviewers

This article was reviewed by Dr. Céline Brochier-Armanet, Dr Kira S Makarova (nominated by Dr. Eugene Koonin), and Dr. Michael Galperin.