A hot story from comparative genomics: reverse gyrase is the only hyperthermophile-specific protein

Adenosine 5'-O-(3-thio)triphosphate (ATPgammaS) promotes positive supercoiling of DNA by T. maritima reverse gyrase

Reverse gyrases are topoisomerases that catalyze ATP-dependent positive supercoiling of circular covalently closed DNA. They consist of an N-terminal helicase-like domain, fused to a C-terminal topoisomerase I-like domain. Most of our knowledge on reverse gyrase-mediated positive DNA supercoiling is based on studies of archaeal enzymes. To identify general and individual properties of reverse gyrases, we set out to characterize the reverse gyrase from a hyperthermophilic eubacterium. Thermotoga maritima reverse gyrase relaxes negatively supercoiled DNA in the presence of ADP or the non-hydrolyzable ATP-analog ADPNP. Nucleotide binding is necessary, but not sufficient for the relaxation reaction. In the presence of ATP, positive supercoils are introduced at temperatures above 50 degrees C. However, ATP hydrolysis is stimulated by DNA already at 37 degrees C, suggesting that reverse gyrase is not frozen at this temperature, but capable of undergoing inter-domain communication. Positive supercoiling by reverse gyrase is strictly coupled to ATP hydrolysis. At the physiological temperature of 75 degrees C, reverse gyrase binds and hydrolyzes ATPgammaS. Surprisingly, ATPgammaS hydrolysis is stimulated by DNA, and efficiently promotes positive DNA supercoiling, demonstrating that inter-domain communication during positive supercoiling is fully functional with both ATP and ATPgammaS. These findings support a model for communication between helicase-like and topoisomerase domains in reverse gyrase, in which an ATP and DNA-induced closure of the cleft in the helicase-like domain initiates a cycle of conformational changes that leads to positive DNA supercoiling.

Comparative genome analysis of four magnetotactic bacteria reveals a complex set of group-specific genes implicated in magnetosome biomineralization and function

Magnetotactic bacteria (MTB) are a heterogeneous group of aquatic prokaryotes with a unique intracellular organelle, the magnetosome, which orients the cell along magnetic field lines. Magnetotaxis is a complex phenotype, which depends on the coordinate synthesis of magnetosomes and the ability to swim and orient along the direction caused by the interaction with the Earth's magnetic field. Although a number of putative magnetotaxis genes were recently identified within a conserved genomic magnetosome island (MAI) of several MTB, their functions have remained mostly unknown, and it was speculated that additional genes located outside the MAI might be involved in magnetosome formation and magnetotaxis. In order to identify genes specifically associated with the magnetotactic phenotype, we conducted comparisons between four sequenced magnetotactic Alphaproteobacteria including the nearly complete genome of Magnetospirillum gryphiswaldense strain MSR-1, the complete genome of Magnetospirillum magneticum strain AMB-1, the complete genome of the magnetic coccus MC-1, and the comparative-ready preliminary genome assembly of Magnetospirillum magnetotacticum strain MS-1 against an in-house database comprising 426 complete bacterial and archaeal genome sequences. A magnetobacterial core genome of about 891 genes was found shared by all four MTB. In addition to a set of approximately 152 genus-specific genes shared by the three Magnetospirillum strains, we identified 28 genes as group specific, i.e., which occur in all four analyzed MTB but exhibit no (MTB-specific genes) or only remote (MTB-related genes) similarity to any genes from nonmagnetotactic organisms and which besides various novel genes include nearly all mam and mms genes previously shown to control magnetosome formation. The MTB-specific and MTB-related genes to a large extent display synteny, partially encode previously unrecognized magnetosome membrane proteins, and are either located within (18 genes) or outside (10 genes) the MAI of M. gryphiswaldense. These genes, which represent less than 1% of the 4,268 open reading frames of the MSR-1 genome, as yet are mostly of unknown functions but are likely to be specifically involved in magnetotaxis and, thus, represent prime targets for future experimental analysis.

Reverse gyrase: an insight into the role of DNA-topoisomerases

Reverse gyrase was discovered more than twenty years ago. Recent biochemical and structural results have greatly enhanced our understanding of their positive supercoiling mechanism. In addition to new biochemical properties, a fine tuning of reverse gyrase regulation in response to DNA damaging agents has been recently described. These data give us a new insight in the cellular role of reverse gyrase. Moreover, it has been proposed that reverse gyrase has been implicated in genome stability.

Origin and evolution of DNA topoisomerases

The DNA topoisomerases are essential for DNA replication, transcription, recombination, as well as for chromosome compaction and segregation. They may have appeared early during the formation of the modern DNA world. Several families and subfamilies of the two types of DNA topoisomerases (I and II) have been described in the three cellular domains of life (Archaea, Bacteria and Eukarya), as well as in viruses infecting eukaryotes or bacteria. The main families of DNA topoisomerases, Topo IA, Topo IB, Topo IC (Topo V), Topo IIA and Topo IIB (Topo VI) are not homologous, indicating that they originated independently. However, some of them share homologous modules or subunits that were probably recruited independently to produce different topoisomerase activities. The puzzling phylogenetic distribution of the various DNA topoisomerase families and subfamilies cannot be easily reconciled with the classical models of early evolution describing the relationships between the three cellular domains. A possible scenario is based on a Last Universal Common Ancestor (LUCA) with a RNA genome (i.e. without the need for DNA topoisomerases). Different families of DNA topoisomerases (some of them possibly of viral origin) would then have been independently introduced in the different cellular domains. We review here the main characteristics of the different families and subfamilies of DNA topoisomerases in a historical and evolutionary perspective, with the hope to stimulate further works and discussions on the origin and evolution of these fascinating enzymes.

Type IA topoisomerases: a simple puzzle?

Type IA topoisomerases are enzymes that can modify DNA topology. They form a distinct family of proteins present in all domains of life, from bacteria to archaea and higher eukaryotes. They are composed of two domains: a core domain containing all the conserved motifs involved in the trans-esterification reactions, and a carboxyl-terminal domain that is highly variable in size and sequence. The latter appears to interact with other proteins, defining the physiological use of the topoisomerase activity. The evolutionary relevance of this topoisomerase-cofactor complex, also known as the "toposome", as well as its enzymatic consequences are discussed in this review.

Characterization of DNA transport in the thermophilic bacterium Thermus thermophilus HB27

Horizontal gene transfer has been a major force for genome plasticity over evolutionary history, and is largely responsible for fitness-enhancing traits, including antibiotic resistance and virulence factors. In particular, for adaptation of prokaryotes to extreme environments, lateral gene transfer seems to have played a crucial role. Recently, by performing a genome-wide mutagenesis approach with Thermus thermophilus HB27, we identified the first genes in a thermophilic bacterium for the uptake of free DNA, a process called natural transformation. Here, we present the first data on the biochemistry and bioenergetics of the DNA transport process in this thermophile. We report that linear and circular plasmid DNA are equally well taken up with a high maximal velocity of 1.5 microg DNA.(mg protein)(-1).min(-1), demonstrating an extremely efficient binding and uptake rate of 40 kb.s(-1).cell(-1). Uncouplers and ATPase inhibitors immediately inhibited DNA uptake, providing clear evidence that DNA translocation in HB27 is an energy-dependent process. DNA uptake studies with genomic DNA of Bacteria, Archaea and Eukarya revealed that Thermus thermophilus HB27 takes up DNA from members of all three domains of life. We propose that the extraordinary broad substrate specificity of the highly efficient Thermus thermophilus HB27 DNA uptake system may contribute significantly to thermoadaptation of Thermus thermophilus HB27 and to interdomain DNA transfer in hot environments.

The origin and evolution of Archaea: a state of the art

Environmental surveys indicate that the Archaea are diverse and abundant not only in extreme environments, but also in soil, oceans and freshwater, where they may fulfil a key role in the biogeochemical cycles of the planet. Archaea display unique capacities, such as methanogenesis and survival at temperatures higher than 90 degrees C, that make them crucial for understanding the nature of the biota of early Earth. Molecular, genomics and phylogenetics data strengthen Woese's definition of Archaea as a third domain of life in addition to Bacteria and Eukarya. Phylogenomics analyses of the components of different molecular systems are highlighting a core of mainly vertically inherited genes in Archaea. This allows recovering a globally well-resolved picture of archaeal evolution, as opposed to what is observed for Bacteria and Eukarya. This may be due to the fact that no rapid divergence occurred at the emergence of present-day archaeal lineages. This phylogeny supports a hyperthermophilic and non-methanogenic ancestor to present-day archaeal lineages, and a profound divergence between two major phyla, the Crenarchaeota and the Euryarchaeota, that may not have an equivalent in the other two domains of life. Nanoarchaea may not represent a third and ancestral archaeal phylum, but a fast-evolving euryarchaeal lineage. Methanogenesis seems to have appeared only once and early in the evolution of Euryarchaeota. Filling up this picture of archaeal evolution by adding presently uncultivated species, and placing it back in geological time remain two essential goals for the future.

Comparative genomics of Thermus thermophilus: Plasticity of the megaplasmid and its contribution to a thermophilic lifestyle

The bacterium Thermus thermophilus grows at temperatures up to 85 degrees C and is equipped with thermostable enzymes of biotechnological interest. The recently decoded genomes of two strains of T. thermophilus, HB27 and HB8, each composed of a chromosome and a megaplasmid, must certainly encode specific strategies to encounter the thermophile challenge. Here, a genome comparison was undertaken to distinguish common functions from the flexible gene pool, which gave some clues about the biological traits involved in a thermophile lifestyle. The chromosomes were highly conserved, with about 100 strain-specific genes probably reflecting adaptations to the corresponding biological niche, such as metabolic specialities and distinct cell surface determinates including type IV pili. The two megaplasmids showed an elevated plasticity. Upon comparison and re-examination of their gene content, both megaplasmids seem to be implicated in assisting thermophilic growth: a large portion of their genes are apparently involved in DNA repair functions. About 30 plasmid-encoded genes exhibit sequence and domain composition similarity to a predicted DNA repair system specific for thermophilic Archaea and bacteria. Moreover, the plasmid-encoded carotenoid biosynthesis gene cluster is interlocked with genes involved in UV-induced DNA damage repair. This illustrates the importance of DNA protection and repair at elevated growth temperatures.

Oligonucleotide cleavage and rejoining by topoisomerase III from the hyperthermophilic archaeon Sulfolobus solfataricus: temperature dependence and strand annealing-promoted DNA religation

Topoisomerase III from the hyperthermophilic archaeon Sulfolobus solfataricus (Sso topo III) is optimally active in DNA relaxation at 75 degrees C. We report here that Sso topo III-catalysed DNA cleavage and religation differed significantly in temperature dependence: the enzyme was most active in cleaving ssDNA containing a cleavage site at 25-50 degrees C, but was efficient in rejoining the cleaved DNA strand only at higher temperatures (e.g. > or = 45 degrees C). The failure of Sso topo III to rejoin the cleaved DNA strand efficiently appeared to be responsible for the inability of the enzyme to relax negatively supercoiled DNA at low temperature (e.g. 25 degrees C). Intriguingly, Sso topo III facilitated DNA annealing although it showed higher affinity for ssDNA than for dsDNA. Religation of the DNA strand cleaved by Sso topo III was drastically enhanced when the DNA was allowed to anneal to a complementary non-cleaved oligonucleotide, presumably as a result of destabilization of the interaction between the enzyme and the cleaved strand through the formation of duplex DNA. A region in the non-cleaved strand corresponding to a sequence containing six bases on the 5' side and two bases on the 3' side of the cleavage site in the cleaved strand was crucial to the annealing-promoted religation. However, the annealing-promoted religation was relatively insensitive to mismatches in this region and the region conserved for oligonucleotide cleavage, except for that at the 5' end of the broken strand. These results suggest that Sso topo III is well suited for a role in DNA rewinding, whether it leads to homoduplex or heteroduplex formation.

Hot crenarchaeal viruses reveal deep evolutionary connections

The discovery of archaeal viruses provides insights into the fundamental biochemistry and evolution of the Archaea. Recent studies have identified a wide diversity of archaeal viruses within the hot springs of Yellowstone National Park and other high-temperature environments worldwide. These viruses are often morphologically unique and code for genes with little similarity to other known genes in the biosphere, a characteristic that has complicated efforts to trace their evolutionary history. Comparative genomics combined with structural analysis indicate that spindle-shaped virus lineages might be unique to the Archaea, whereas other icosahedral viruses might share a common lineage with viruses of Bacteria and Eukarya. These studies provide insights into the evolutionary history of viruses in all three domains of life.

Use of a multi-way method to analyze the amino acid composition of a conserved group of orthologous proteins in prokaryotes

BACKGROUND

Amino acids in proteins are not used equally. Some of the differences in the amino acid composition of proteins are between species (mainly due to nucleotide composition and lifestyle) and some are between proteins from the same species (related to protein function, expression or subcellular localization, for example). As several factors contribute to the different amino acid usage in proteins, it is difficult both to analyze these differences and to separate the contributions made by each factor.

RESULTS

Using a multi-way method called Tucker3, we have analyzed the amino composition of a set of 64 orthologous groups of proteins present in 62 archaea and bacteria. This dataset corresponds to essential proteins such as ribosomal proteins, tRNA synthetases and translational initiation or elongation factors, which are common to all the species analyzed. The Tucker3 model can be used to study the amino acid variability within and between species by taking into consideration the tridimensionality of the data set. We found that the main factor behind the amino acid composition of proteins is independent of the organism or protein function analyzed. This factor must be related to the biochemical characteristics of each amino acid. The difference between the non-ribosomal proteins and the ribosomal proteins (which are rich in arginine and lysine) is the main factor behind the differences in amino acid composition within species, while G+C content and optimal growth temperature are the main factors behind the differences in amino acid usage between species.

CONCLUSION

We show that a multi-way method is useful for comparing the amino acid composition of several groups of orthologous proteins from the same group of species. This kind of dataset is extremely useful for detecting differences between and within species.

Selective degradation of reverse gyrase and DNA fragmentation induced by alkylating agent in the archaeon Sulfolobus solfataricus

Reverse gyrase is a peculiar DNA topoisomerase, specific of hyperthermophilic Archaea and Bacteria, which has the unique ability of introducing positive supercoiling into DNA molecules. Although the function of the enzyme has not been established directly, it has been suggested to be involved in DNA protection and repair. We show here that the enzyme is degraded after treatment of Sulfolobus solfataricus cells with the alkylating agent MMS. MMS-induced reverse gyrase degradation is highly specific, since (i) neither hydroxyurea (HU) nor puromycin have a similar effect, and (ii) topoisomerase VI and two chromatin components are not degraded. Reverse gyrase degradation does not depend on protein synthesis. Experiments in vitro show that direct exposure of cell extracts to MMS does not induce reverse gyrase degradation; instead, extracts from MMS-treated cells contain some factor(s) able to degrade the enzyme in extracts from control cells. In vitro, degradation is blocked by incubation with divalent metal chelators, suggesting that reverse gyrase is selectively degraded by a metal-dependent protease in MMS-treated cells. In addition, we find a striking concurrence of extensive genomic DNA degradation and reverse gyrase loss in MMS-treated cells. These results support the hypothesis that reverse gyrase plays an essential role in DNA thermoprotection and repair in hyperthermophilic organisms.

On the correlation between genomic G+C content and optimal growth temperature in prokaryotes: data quality and confounding factors

The correlation between genomic G+C content and optimal growth temperature in prokaryotes has gained renewed interest after Musto et al. [H. Musto, H. Naya, A. Zavala, H. Romero, F. Alvarex-Valin, G. Bernardi, Correlations between genomic GC levels and optimal growth temperatures in prokaryotes, FEBS Lett. 573 (2004) 73-77], reported that positive correlations exist in 15 families studied. We have reanalyzed their data and found that when genome size and data quality were adjusted for, there was no significant evidence of relationship between optimal temperature and GC content for two of the families that had previously shown strongly significant correlations. Using updated temperature optima for Halobacteriaceae species we found the correlation is insignificant in this family. For the family Enterobacteriaceae when genome size and optimal temperature are included in a multiple linear regression, only genome size is significant as a predictor of GC content. We showed that more profound statistical methods than simple two factor correlation analysis should be used for analyzing complex intrinsic and extrinsic factors that affect genomic GC content. We further found that a positive correlation between temperature and genomic GC is only evident in free-living species of low optimal growth temperatures.

An atypical type II topoisomerase from Mycobacterium smegmatis with positive supercoiling activity

Topoisomerases are essential ubiquitous enzymes, falling into two distinct classes. A number of eubacteria including Escherichia coli, typically contain four topoisomerases, two type I topoisomerases and two type II topoisomerases viz. DNA gyrase and topoisomerase IV. In contrast several other bacterial genomes including mycobacteria, encode for one type I topoisomerase and a DNA gyrase. Here we describe a new type II topoisomerase from Mycobacterium smegmatis which is different from DNA gyrase or topoisomerase IV in its characteristics and origin. The topoisomerase is distinct with respect to domain organization, properties and drug sensitivity. The enzyme catalyses relaxation of negatively supercoiled DNA in an ATP-dependent manner and also introduces positive supercoils to both relaxed and negatively supercoiled substrates. The genes for this additional topoisomerase are not found in other sequenced mycobacterial genomes and may represent a distant lineage.

Reverse gyrase functions as a DNA renaturase: annealing of complementary single-stranded circles and positive supercoiling of a bubble substrate

Reverse gyrase is a hyperthermophile-specific enzyme that can positively supercoil DNA concomitant with ATP hydrolysis. However, the DNA supercoiling activity is inefficient and requires an excess amount of enzyme relative to DNA. We report here several activities that reverse gyrase can efficiently mediate with a substoichiometric amount of enzyme. In the presence of a nucleotide cofactor, reverse gyrase can readily relax negative supercoils, but not the positive ones, from a plasmid DNA substrate. Reverse gyrase can completely relax positively supercoiled DNA, provided that the DNA substrate contains a single-stranded bubble. Reverse gyrase efficiently anneals complementary single-stranded circles. A substoichiometric amount of reverse gyrase can insert positive supercoils into DNA with a single-stranded bubble, in contrast to plasmid DNA substrate. We have designed a novel method based on phage-mid DNA vectors to prepare a circular DNA substrate containing a single-stranded bubble with defined length and sequence. With these bubble DNA substrates, we demonstrated that efficient positive supercoiling by reverse gyrase requires a bubble size larger than 20 nucleotides. The activities of annealing single-stranded DNA circles and positive supercoiling of bubble substrate demonstrate that reverse gyrase can function as a DNA renaturase. These biochemical activities also suggest that reverse gyrase can have an important biological function in sensing and eliminating unpaired regions in the genome of a hyperthermophilic 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.

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.

Human topoisomerase IIalpha rapidly relaxes positively supercoiled DNA: implications for enzyme action ahead of replication forks

Movement of the DNA replication machinery through the double helix induces acute positive supercoiling ahead of the fork and precatenanes behind it. Because topoisomerase I and II create transient single- and double-stranded DNA breaks, respectively, it has been assumed that type I enzymes relax the positive supercoils that precede the replication fork. Conversely, type II enzymes primarily resolve the precatenanes and untangle catenated daughter chromosomes. However, studies on yeast and bacteria suggest that type II topoisomerases may also function ahead of the replication machinery. If this is the case, then positive DNA supercoils should be the preferred relaxation substrate for topoisomerase IIalpha, the enzyme isoform involved in replicative processes in humans. Results indicate that human topoisomerase IIalpha relaxes positively supercoiled plasmids >10-fold faster than negatively supercoiled molecules. In contrast, topoisomerase IIbeta, which is not required for DNA replication, displays no such preference. In addition to its high rates of relaxation, topoisomerase IIalpha maintains lower levels of DNA cleavage complexes with positively supercoiled molecules. These properties suggest that human topoisomerase IIalpha has the potential to alleviate torsional stress ahead of replication forks in an efficient and safe manner.

Horizontal transfer of two operons coding for hydrogenases between bacteria and archaea

Using a phylogenetic approach, we discovered three putative horizontal transfers between bacterial and archaeal species involving large clusters of genes. One transfer involves an operon of 13 genes, called mbx, which probably was transferred into the genome of Thermotoga maritima from a species belonging or close to the Pyrococcus genus. The two others implied an operon of six genes, called ech, transferred independently to the genomes of Thermoanaerobacter tengcongensis and Desulfovibrio gigas, from a species belonging or close to the Methanosarcina genus. All these transfers affected operons coding for multisubunit membrane-bound (NiFe) hydrogenases involved in the energy metabolism of the donor genomes. The functionality of the transferred operons has not been experimentally demonstrated for T. maritima, whereas in D. gigas and T. tengcongensis the encoded multisubunit hydrogenase could have a role in energy conservation. This report adds several cases of horizontal gene transfers among hydrogenases already described.

The genomics of disulfide bonding and protein stabilization in thermophiles

Thermophilic organisms flourish in varied high-temperature environmental niches that are deadly to other organisms. Recently, genomic evidence has implicated a critical role for disulfide bonds in the structural stabilization of intracellular proteins from certain of these organisms, contrary to the conventional view that structural disulfide bonds are exclusively extracellular. Here both computational and structural data are presented to explore the occurrence of disulfide bonds as a protein-stabilization method across many thermophilic prokaryotes. Based on computational studies, disulfide-bond richness is found to be widespread, with thermophiles containing the highest levels. Interestingly, only a distinct subset of thermophiles exhibit this property. A computational search for proteins matching this target phylogenetic profile singles out a specific protein, known as protein disulfide oxidoreductase, as a potential key player in thermophilic intracellular disulfide-bond formation. Finally, biochemical support in the form of a new crystal structure of a thermophilic protein with three disulfide bonds is presented together with a survey of known structures from the literature. Together, the results provide insight into biochemical specialization and the diversity of methods employed by organisms to stabilize their proteins in exotic environments. The findings also motivate continued efforts to sequence genomes from divergent organisms.

Certain thermophiles are found to stabilize their proteins in extreme environments with additional disulfide bonds. A phylogenetic profile identifies a protein disulfide oxidoreductase critical to the stabilization process.

Nanoarchaea: representatives of a novel archaeal phylum or a fast-evolving euryarchaeal lineage related to Thermococcales?

An analysis of the position of Nanoarcheum equitans in the archaeal phylogeny using a large dataset of concatenated ribosomal proteins from 25 archaeal genomes suggests that N. equitans is likely to be the representative of a fast-evolving euryarchaeal lineage.

Background

Cultivable archaeal species are assigned to two phyla - the Crenarchaeota and the Euryarchaeota - by a number of important genetic differences, and this ancient split is strongly supported by phylogenetic analysis. The recently described hyperthermophile Nanoarchaeum equitans, harboring the smallest cellular genome ever sequenced (480 kb), has been suggested as the representative of a new phylum - the Nanoarchaeota - that would have diverged before the Crenarchaeota/Euryarchaeota split. Confirming the phylogenetic position of N. equitans is thus crucial for deciphering the history of the archaeal domain.

Results

We tested the placement of N. equitans in the archaeal phylogeny using a large dataset of concatenated ribosomal proteins from 25 archaeal genomes. We indicate that the placement of N. equitans in archaeal phylogenies on the basis of ribosomal protein concatenation may be strongly biased by the coupled effect of its above-average evolutionary rate and lateral gene transfers. Indeed, we show that different subsets of ribosomal proteins harbor a conflicting phylogenetic signal for the placement of N. equitans. A BLASTP-based survey of the phylogenetic pattern of all open reading frames (ORFs) in the genome of N. equitans revealed a surprisingly high fraction of close hits with Euryarchaeota, notably Thermococcales. Strikingly, a specific affinity of N. equitans and Thermococcales was strongly supported by phylogenies based on a subset of ribosomal proteins, and on a number of unrelated molecular markers.

Conclusion

We suggest that N. equitans may more probably be the representative of a fast-evolving euryarchaeal lineage (possibly related to Thermococcales) than the representative of a novel and early diverging archaeal phylum.

Nucleotide- and stoichiometry-dependent DNA supercoiling by reverse gyrase

Reverse gyrase is a unique type IA topoisomerase that can introduce positive supercoils into DNA. We have investigated some of the biochemical properties of Archaeoglobus fulgidus reverse gyrase. It can mediate three distinct supercoiling reactions depending on the adenine nucleotide cofactor that is present in the reaction. Besides the ATP-driven positive supercoiling reaction, the enzyme can introduce negative supercoils with a nonhydrolyzable analog, adenylyl imidodiphosphate. In the presence of ADP the plasmid DNA is relaxed almost completely, leaving a very low level of positive supercoiling. Surprisingly, the final supercoiling extent for all three distinct reactions depends on the stoichiometry of enzyme to DNA. This dependence is not due to the difference of reaction rate, suggesting that the amount of enzyme bound to DNA is an important determinant for the final supercoiling state of the reaction product. Reverse gyrase also displays exquisite sensitivity toward temperature. Raising the reaction temperatures from 80 to 85 degrees C, both of which are within the optimal growth temperature of A. fulgidus, greatly increases enzyme activity for all the supercoiling reactions. For the reaction with AMPPNP, the product is a hypernegatively supercoiled DNA. This dramatic enhancement of the reverse gyrase activity is also correlated with the appearance of DNA in a pre-melting state at 85 degrees C, likely due to the presence of extensively unwound regions in the plasmid. The possible mechanistic insights from these findings will be presented here.

Functional interaction of reverse gyrase with single-strand binding protein of the archaeon Sulfolobus

Reverse gyrase is a unique hyperthermophile-specific DNA topoisomerase that induces positive supercoiling. It is a modular enzyme composed of a topoisomerase IA and a helicase domain, which cooperate in the ATP-dependent positive supercoiling reaction. Although its physiological function has not been determined, it can be hypothesized that, like the topoisomerase–helicase complexes found in every organism, reverse gyrase might participate in different DNA transactions mediated by multiprotein complexes. Here, we show that reverse gyrase activity is stimulated by the single-strand binding protein (SSB) from the archaeon Sulfolobus solfataricus. Using a combination of in vitro assays we analysed each step of the complex reverse gyrase reaction. SSB stimulates all the steps of the reaction: binding to DNA, DNA cleavage, strand passage and ligation. By co-immunoprecipitation of cell extracts we show that reverse gyrase and SSB assemble a complex in the presence of DNA, but do not make stable protein–protein interactions. In addition, SSB stimulates reverse gyrase positive supercoiling activity on DNA templates associated with the chromatin protein Sul7d. Furthermore, SSB enhances binding and cleavage of UV-irradiated substrates by reverse gyrase. The results shown here suggest that these functional interactions may have biological relevance and that the interplay of different DNA binding proteins might modulate reverse gyrase activity in DNA metabolic pathways.

Archaeal genetics - the third way

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

Cold shock of a hyperthermophilic archaeon: Pyrococcus furiosus exhibits multiple responses to a suboptimal growth temperature with a key role for membrane-bound glycoproteins

The hyperthermophilic archaeon, Pyrococcus furiosus, was grown on maltose near its optimal growth temperature, 95 degrees C, and at the lower end of the temperature range for significant growth, 72 degrees C. In addition, cultures were shocked by rapidly dropping the temperature from 95 to 72 degrees C. This resulted in a 5-h lag phase, during which time little growth occurred. Transcriptional analyses using whole-genome DNA microarrays representing 2,065 open reading frames (ORFs) in the P. furiosus genome showed that cells undergo three very different responses at 72 degrees C: an early shock (1 to 2 h), a late shock (5 h), and an adapted response (occurring after many generations at 72 degrees C). Each response involved the up-regulation in the expression of more than 30 ORFs unique to that response. These included proteins involved in translation, solute transport, amino acid biosynthesis, and tungsten and intermediary carbon metabolism, as well as numerous conserved-hypothetical and/or membrane-associated proteins. Two major membrane proteins were evident after one-dimensional sodium dodecyl sulfate-gel analysis of cold-adapted cells, and staining revealed them to be glycoproteins. Their cold-induced expression evident from the DNA microarray analysis was confirmed by quantitative PCR. Termed CipA (PF0190) and CipB (PF1408), both appear to be solute-binding proteins. While the archaea do not contain members of the bacterial cold shock protein (Csp) family, they all contain homologs of CipA and CipB. These proteins are also related phylogenetically to some cold-responsive genes recently identified in certain bacteria. The Cip proteins may represent a general prokaryotic-type cold response mechanism that is present even in hyperthermophilic archaea.

Thermoadaptation trait revealed by the genome sequence of thermophilic Geobacillus kaustophilus

We present herein the first complete genome sequence of a thermophilic Bacillus-related species, Geobacillus kaustophilus HTA426, which is composed of a 3.54 Mb chromosome and a 47.9 kb plasmid, along with a comparative analysis with five other mesophilic bacillar genomes. Upon orthologous grouping of the six bacillar sequenced genomes, it was found that 1257 common orthologous groups composed of 1308 genes (37%) are shared by all the bacilli, whereas 839 genes (24%) in the G.kaustophilus genome were found to be unique to that species. We were able to find the first prokaryotic sperm protamine P1 homolog, polyamine synthase, polyamine ABC transporter and RNA methylase in the 839 unique genes; these may contribute to thermophily by stabilizing the nucleic acids. Contrasting results were obtained from the principal component analysis (PCA) of the amino acid composition and synonymous codon usage for highlighting the thermophilic signature of the G.kaustophilus genome. Only in the PCA of the amino acid composition were the Bacillus-related species located near, but were distinguishable from, the borderline distinguishing thermophiles from mesophiles on the second principal axis. Further analysis revealed some asymmetric amino acid substitutions between the thermophiles and the mesophiles, which are possibly associated with the thermoadaptation of the organism.

A novel nuclease-ATPase (Nar71) from archaea is part of a proposed thermophilic DNA repair system

We have identified a novel structure-specific nuclease in highly fractionated extracts of the thermophilic archaeon Methanothermobacter thermautotrophicus (Mth). The 71 kDa protein product of open reading frame mth1090 is a nuclease with ATPase activity, which we call Nar71 (Nuclease-ATPase in Repair, 71 kDa). The nar71 gene is located in a gene neighbourhood proposed by genomics to encode a novel DNA repair system conserved in thermophiles. The biochemical characterization of Nar71 presented here is the first analysis from within this neighbourhood, and it supports the insight from genomics. Nuclease activity of Nar71 is specific for 3' flaps and flayed duplexes, targeting single-stranded DNA (ssDNA) regions. This activity requires Mg2+ or Mn2+ and is greatly reduced in ATP. In ATP, Nar71 displaces ssDNA, also with high specificity for 3' flap and flayed duplex DNA. Strand displacement is weak compared with nuclease activity, but in ATPS it is abolished, suggesting that Nar71 couples ATP hydrolysis to DNA strand separation. ATPase assays confirmed that Nar71 is stimulated by ssDNA, though not double-stranded DNA. Mutation of Lys-117 in Nar71 abolished ATPase and nuclease activity, and we describe a separation-of-function mutant (K68A) that has lost ATPase activity but retains nuclease activity. A model of possible Nar71 function in DNA repair is presented.

PCOGR: phylogenetic COG ranking as an online tool to judge the specificity of COGs with respect to freely definable groups of organisms

Background

The rapidly increasing number of completely sequenced genomes led to the establishment of the COG-database which, based on sequence homologies, assigns similar proteins from different organisms to clusters of orthologous groups (COGs). There are several bioinformatic studies that made use of this database to determine (hyper)thermophile-specific proteins by searching for COGs containing (almost) exclusively proteins from (hyper)thermophilic genomes. However, public software to perform individually definable group-specific searches is not available.

Results

The tool described here exactly fills this gap. The software is accessible at and is linked to the COG-database. The user can freely define two groups of organisms by selecting for each of the (current) 66 organisms to belong either to groupA, to the reference groupB or to be ignored by the algorithm. Then, for all COGs a specificity index is calculated with respect to the specificity to groupA, i. e. high scoring COGs contain proteins from the most of groupA organisms while proteins from the most organisms assigned to groupB are absent. In addition to ranking all COGs according to the user defined specificity criteria, a graphical visualization shows the distribution of all COGs by displaying their abundance as a function of their specificity indexes.

Conclusions

This software allows detecting COGs specific to a predefined group of organisms. All COGs are ranked in the order of their specificity and a graphical visualization allows recognizing (i) the presence and abundance of such COGs and (ii) the phylogenetic relationship between groupA- and groupB-organisms. The software also allows detecting putative protein-protein interactions, novel enzymes involved in only partially known biochemical pathways, and alternate enzymes originated by convergent evolution.

Genomic and proteomic adaptations to growth at high temperature

Most positively selected mutations cause changes in metabolism, resulting in a better-adapted phenotype. But as well as acting on the information content of genes, natural selection may also act directly on nucleic acid and protein molecules. We review the evidence for direct temperature-dependent natural selection acting on genomes, transcriptomes and proteomes.

Reverse gyrase is not a prerequisite for hyperthermophilic life

We disrupted the reverse gyrase gene from a hyperthermophilic archaeon, Thermococcus kodakaraensis KOD1. An apparent positive supercoiling activity that was observed in the host strain was not found in the disruptant strain. We found that a lack of reverse gyrase led to a retardation in growth that was more striking at higher temperatures. However, the disruption of the reverse gyrase gene did not lead to a lethal phenotype at 90 degrees C. This study provides experimental evidence that reverse gyrase is not a prerequisite for hyperthermophilic life.

Reverse Gyrase Recruitment to DNA after UV Light Irradiation in Sulfolobus solfataricus

Induction of DNA damage triggers a complex biological response concerning not only repair systems but also virtually every cell function. DNA topoisomerases regulate the level of DNA supercoiling in all DNA transactions. Reverse gyrase is a peculiar DNA topoisomerase, specific to hyperthermophilic microorganisms, which contains a helicase and a topoisomerase IA domain that has the unique ability to introduce positive supercoiling into DNA molecules. We show here that reverse gyrase of the archaean Sulfolobus solfataricus is mobilized to DNA in vivo after UV irradiation. The enzyme, either purified or in cell extracts, forms stable covalent complexes with UV-damaged DNA in vitro. We also show that the reverse gyrase translocation to DNA in vivo and the stabilization of covalent complexes in vitro are specific effects of UV light irradiation and do not occur with the intercalating agent actinomycin D. Our results suggest that reverse gyrase might participate, directly or indirectly, in the cell response to UV light-induced DNA damage. This is the first direct evidence of the recruitment of a topoisomerase IA enzyme to DNA after the induction of DNA damage. The interaction between helicase and topoisomerase activities has been previously proposed to facilitate aspects of DNA replication or recombination in both Bacteria and Eukarya. Our results suggest a general role of the association of such activities in maintaining genome integrity and a mutual effect of DNA topology and repair.

Reverse gyrase has heat-protective DNA chaperone activity independent of supercoiling

Hyperthermophilic organisms must protect their constituent macromolecules from heat-induced degradation. A general mechanism for thermoprotection of DNA in active cells is unknown. We show that reverse gyrase, the only protein that is both specific and common to all hyperthermophiles, reduces the rate of double-stranded DNA breakage approximately 8-fold at 90 degrees C. This activity does not require ATP hydrolysis and is independent of the positive supercoiling activity of the enzyme. Reverse gyrase has a minor nonspecific effect on the rate of depurination, and a major specific effect on the rate of double-strand breakage. Using electron microscopy, we show that reverse gyrase recognizes nicked DNA and recruits a protein coat to the site of damage through cooperative binding. Analogously to molecular chaperones that assist unfolded proteins, we found that reverse gyrase prevents inappropriate aggregation of denatured DNA regions and promotes correct annealing. We propose a model for a targeted protection mechanism in vivo in which reverse gyrase detects damaged DNA and acts as a molecular splint to prevent DNA breakage in the vicinity of the lesion, thus maintaining damaged DNA in a conformation that is amenable to repair.

Lateral gene transfer and the origins of prokaryotic groups

Lateral gene transfer (LGT) is now known to be a major force in the evolution of prokaryotic genomes. To date, most analyses have focused on either (a) verifying phylogenies of individual genes thought to have been transferred, or (b) estimating the fraction of individual genomes likely to have been introduced by transfer. Neither approach does justice to the ability of LGT to effect massive and complex transformations in basic biology. In some cases, such transformation will be manifested as the patchy distribution of a seemingly fundamental property (such as aerobiosis or nitrogen fixation) among the members of a group classically defined by the sharing of other properties (metabolic, morphological, or molecular, such as small subunit ribosomal RNA sequence). In other cases, the lineage of recipients so transformed may be seen to comprise a new group of high taxonomic rank ("class" or even "phylum"). Here we review evidence for an important role of LGT in the evolution of photosynthesis, aerobic respiration, nitrogen fixation, sulfate reduction, methylotrophy, isoprenoid biosynthesis, quorum sensing, flotation (gas vesicles), thermophily, and halophily. Sometimes transfer of complex gene clusters may have been involved, whereas other times separate exchanges of many genes must be invoked.

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

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

Thermophile-specific proteins: the gene product of aq_1292 from Aquifex aeolicus is an NTPase

Background

To identify thermophile-specific proteins, we performed phylogenetic patterns searches of 66 completely sequenced microbial genomes. This analysis revealed a cluster of orthologous groups (COG1618) which contains a protein from every thermophile and no sequence from 52 out of 53 mesophilic genomes. Thus, COG1618 proteins belong to the group of thermophile-specific proteins (THEPs) and therefore we here designate COG1618 proteins as THEP1s. Since no THEP1 had been analyzed biochemically thus far, we characterized the gene product of aq_1292 which is THEP1 from the hyperthermophilic bacterium Aquifex aeolicus (aaTHEP1).

Results

aaTHEP1 was cloned in E. coli, expressed and purified to homogeneity. At a temperature optimum between 70 and 80°C, aaTHEP1 shows enzymatic activity in hydrolyzing ATP to ADP + P_i with k_cat = 5 × 10^-3 s^-1 and K_m = 5.5 × 10^-6 M. In addition, the enzyme exhibits GTPase activity (k_cat = 9 × 10^-3 s^-1 and K_m= 45 × 10^-6 M). aaTHEP1 is inhibited competitively by CTP, UTP, dATP, dGTP, dCTP, and dTTP. As shown by gel filtration, aaTHEP1 in its purified state appears as a monomer. The enzyme is resistant to limited proteolysis suggesting that it consists of a single domain. Although THEP1s are annotated as "predicted nucleotide kinases" we could not confirm such an activity experimentally.

Conclusion

Since aaTHEP1 is the first member of COG1618 that is characterized biochemically and functional information about one member of a COG may be transferred to the entire COG, we conclude that COG1618 proteins are a family of thermophilic NTPases.

OrthoMCL: identification of ortholog groups for eukaryotic genomes

The identification of orthologous groups is useful for genome annotation, studies on gene/protein evolution, comparative genomics, and the identification of taxonomically restricted sequences. Methods successfully exploited for prokaryotic genome analysis have proved difficult to apply to eukaryotes, however, as larger genomes may contain multiple paralogous genes, and sequence information is often incomplete. OrthoMCL provides a scalable method for constructing orthologous groups across multiple eukaryotic taxa, using a Markov Cluster algorithm to group (putative) orthologs and paralogs. This method performs similarly to the INPARANOID algorithm when applied to two genomes, but can be extended to cluster orthologs from multiple species. OrthoMCL clusters are coherent with groups identified by EGO, but improved recognition of "recent" paralogs permits overlapping EGO groups representing the same gene to be merged. Comparison with previously assigned EC annotations suggests a high degree of reliability, implying utility for automated eukaryotic genome annotation. OrthoMCL has been applied to the proteome data set from seven publicly available genomes (human, fly, worm, yeast, Arabidopsis, the malaria parasite Plasmodium falciparum, and Escherichia coli). A Web interface allows queries based on individual genes or user-defined phylogenetic patterns (http://www.cbil.upenn.edu/gene-family). Analysis of clusters incorporating P. falciparum genes identifies numerous enzymes that were incompletely annotated in first-pass annotation of the parasite genome.

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

Archaea comprise one of the three distinct domains of life (with bacteria and eukaryotes). With 16 complete archaeal genomes sequenced to date, comparative genomics has revealed a conserved core of 313 genes that are represented in all sequenced archaeal genomes, plus a variable 'shell' that is prone to lineage-specific gene loss and horizontal gene exchange. The majority of archaeal genes have not been experimentally characterized, but novel functional pathways have been predicted.

Potential genomic determinants of hyperthermophily

We searched for genes that could be important for hyperthermophily using a flexible approach to phyletic pattern analysis. We identified 290 clusters of orthologous groups of proteins (COGs) that are preferentially present in archaeal and bacterial hyperthermophiles. Of these, 58 COGs include proteins from at least one bacterium and two archaea, and these were considered to be the best candidates for a specific association with the hyperthermophilic phenotype. Detailed sequence and genome-context analysis of these COGs led to functional predictions for several previously uncharacterized protein families, including a novel group of putative molecular chaperones and a unique transcriptional regulator.

Phylogenomics of type II DNA topoisomerases

Type II DNA topoisomerases (Topo II) are essential enzymes implicated in key nuclear processes. The recent discovery of a novel kind of Topo II (DNA topoisomerase VI) in Archaea led to a division of these enzymes into two non-homologous families, (Topo IIA and Topo IIB) and to the identification of the eukaryotic protein that initiates meiotic recombination, Spo11. In the present report, we have updated the distribution of all Topo II in the three domains of life by a phylogenomic approach. Both families exhibit an atypical distribution by comparison with other informational proteins, with predominance of Topo IIA in Bacteria, Eukarya and viruses, and Topo IIB in Archaea. However, plants and some Archaea contain Topo II from both families. We confront this atypical distribution with current hypotheses on the evolution of the three domains of life and origin of DNA genomes.

Algorithms for computing parsimonious evolutionary scenarios for genome evolution, the last universal common ancestor and dominance of horizontal gene transfer in the evolution of prokaryotes

BACKGROUND

Comparative analysis of sequenced genomes reveals numerous instances of apparent horizontal gene transfer (HGT), at least in prokaryotes, and indicates that lineage-specific gene loss might have been even more common in evolution. This complicates the notion of a species tree, which needs to be re-interpreted as a prevailing evolutionary trend, rather than the full depiction of evolution, and makes reconstruction of ancestral genomes a non-trivial task.

RESULTS

We addressed the problem of constructing parsimonious scenarios for individual sets of orthologous genes given a species tree. The orthologous sets were taken from the database of Clusters of Orthologous Groups of proteins (COGs). We show that the phyletic patterns (patterns of presence-absence in completely sequenced genomes) of almost 90% of the COGs are inconsistent with the hypothetical species tree. Algorithms were developed to reconcile the phyletic patterns with the species tree by postulating gene loss, COG emergence and HGT (the latter two classes of events were collectively treated as gene gains). We prove that each of these algorithms produces a parsimonious evolutionary scenario, which can be represented as mapping of loss and gain events on the species tree. The distribution of the evolutionary events among the tree nodes substantially depends on the underlying assumptions of the reconciliation algorithm, e.g. whether or not independent gene gains (gain after loss after gain) are permitted. Biological considerations suggest that, on average, gene loss might be a more likely event than gene gain. Therefore different gain penalties were used and the resulting series of reconstructed gene sets for the last universal common ancestor (LUCA) of the extant life forms were analysed. The number of genes in the reconstructed LUCA gene sets grows as the gain penalty increases. However, qualitative examination of the LUCA versions reconstructed with different gain penalties indicates that, even with a gain penalty of 1 (equal weights assigned to a gain and a loss), the set of 572 genes assigned to LUCA might be nearly sufficient to sustain a functioning organism. Under this gain penalty value, the numbers of horizontal gene transfer and gene loss events are nearly identical. This result holds true for two alternative topologies of the species tree and even under random shuffling of the tree. Therefore, the results seem to be compatible with approximately equal likelihoods of HGT and gene loss in the evolution of prokaryotes.

CONCLUSIONS

The notion that gene loss and HGT are major aspects of prokaryotic evolution was supported by quantitative analysis of the mapping of the phyletic patterns of COGs onto a hypothetical species tree. Algorithms were developed for constructing parsimonious evolutionary scenarios, which include gene loss and gain events, for orthologous gene sets, given a species tree. This analysis shows, contrary to expectations, that the number of predicted HGT events that occurred during the evolution of prokaryotes might be approximately the same as the number of gene losses. The approach to the reconstruction of evolutionary scenarios employed here is conservative with regard to the detection of HGT because only patterns of gene presence-absence in sequenced genomes are taken into account. In reality, horizontal transfer might have contributed to the evolution of many other genes also, which makes it a dominant force in prokaryotic evolution.

Evolution of the Archaea

Archaea, members of the third domain of life, are bacterial-looking prokaryotes that harbour many unique genotypic and phenotypic properties, testifying for their peculiar evolutionary status. The archaeal ancestor was probably a hyperthermophilic anaerobe. Two archaeal phyla are presently recognized, the Euryarchaeota and the Crenarchaeota. Methanogenesis was the main invention that occurred in the euryarchaeal phylum and is now shared by several archaeal groups. Adaptation to aerobic conditions occurred several times independently in both Euryarchaeota and Crenarchaeota. Recently, many new groups of Archaea that have not yet been cultured have been detected by PCR amplification of 16S ribosomal RNA from environmental samples. The phenotypic and genotypic characterization of these new groups is now a top priority for further studies on archaeal evolution.