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.     

Archaeal nucleosome positioning sequence from Methanothermus fervidus.

DNA in Methanothermus fervidus, a hyperthermophilic archaeon, is constrained into archaeal nucleosomes in vivo by the archaeal histones HMfA and HMfB. Here, we document the translational and rotational positioning of archaeal nucleosome assembly in vitro by a sequence from the 7S RNA encoding region of the M. fervidus genome. The minor groove of the DNA at the center of the DNA sequence, protected from micrococcal nuclease digestion by incorporation into a positioned archaeal nucleosome, faces away from the archaeal histone core.     

Casposon integration shows strong target site preference and recapitulates protospacer integration by CRISPR-Cas systems

Casposons are a recently discovered group of large DNA transposons present in diverse bacterial and archaeal genomes. For integration into the host chromosome, casposons employ an endonuclease that is homologous to the Cas1 protein involved in protospacer integration by the CRISPR-Cas adaptive immune system. Here we describe the site-preference of integration by the Cas1 integrase (casposase) encoded by the casposon of the archaeon Aciduliprofundum boonei. Oligonucleotide duplexes derived from the terminal inverted repeats (TIR) of the A. boonei casposon as well as mini-casposons flanked by the TIR inserted preferentially at a site reconstituting the original A. boonei target site. As in the A. boonei genome, the insertion was accompanied by a 15-bp direct target site duplication (TSD). The minimal functional target consisted of the 15-bp TSD segment and the adjacent 18-bp sequence which comprises the 3′ end of the tRNA-Pro gene corresponding to the TΨC loop. The functional casposase target site bears clear resemblance to the leader sequence-repeat junction which is the target for protospacer integration catalyzed by the Cas1–Cas2 adaptation module of CRISPR-Cas. These findings reinforce the mechanistic similarities and evolutionary connection between the casposons and the adaptation module of the prokaryotic adaptive immunity systems.     

Gene decay in archaea

The gene-dense chromosomes of archaea and bacteria were long thought to be devoid of pseudogenes, but with the massive increase in available genome sequences, whole genome comparisons between closely related species have identified mutations that have rendered numerous genes inactive. Comparative analyses of sequenced archaeal genomes revealed numerous pseudogenes, which can constitute up to 8.6% of the annotated coding sequences in some genomes. The largest proportion of pseudogenes is created by gene truncations, followed by frameshift mutations. Within archaeal genomes, large numbers of pseudogenes contain more than one inactivating mutation, suggesting that pseudogenes are deleted from the genome more slowly in archaea than in bacteria. Although archaea seem to retain pseudogenes longer than do bacteria, most archaeal genomes have unique repertoires of pseudogenes.     

Genomics of bacteria and archaea: the emerging dynamic view of the prokaryotic world

The first bacterial genome was sequenced in 1995, and the first archaeal genome in 1996. Soon after these breakthroughs, an exponential rate of genome sequencing was established, with a doubling time of approximately 20 months for bacteria and approximately 34 months for archaea. Comparative analysis of the hundreds of sequenced bacterial and dozens of archaeal genomes leads to several generalizations on the principles of genome organization and evolution. A crucial finding that enables functional characterization of the sequenced genomes and evolutionary reconstruction is that the majority of archaeal and bacterial genes have conserved orthologs in other, often, distant organisms. However, comparative genomics also shows that horizontal gene transfer (HGT) is a dominant force of prokaryotic evolution, along with the loss of genetic material resulting in genome contraction. A crucial component of the prokaryotic world is the mobilome, the enormous collection of viruses, plasmids and other selfish elements, which are in constant exchange with more stable chromosomes and serve as HGT vehicles. Thus, the prokaryotic genome space is a tightly connected, although compartmentalized, network, a novel notion that undermines the ‘Tree of Life’ model of evolution and requires a new conceptual framework and tools for the study of prokaryotic evolution.     

Archaeal genomics.

Four euryarchaeal genomes have been completely sequenced and are publicly available: Methanococcus jannaschii, Methanobacterium thermoautotrophicum, Pyrococcus horikoshii and Archaeoglobus fulgidus. Four more genome sequences, two crenarchaeal and two pyrococci, will soon be released. In addition, seven more archaeal genome sequencing projects are under way, including two halophiles, two Thermoplasma, and a methanogen. These projects cover all branches of the archaeal domain and will lead to new insights into archaeal metabolism, DNA processing, and evolutionary relationships with the Bacteria and Eukarya.     

An ssDNA virus infecting archaea: a new lineage of viruses with a membrane envelope

Archaeal organisms are generally known as diverse extremophiles, but they play a crucial role also in moderate environments. So far, only about 50 archaeal viruses have been described in some detail. Despite this, unusual viral morphotypes within this group have been reported. Interestingly, all isolated archaeal viruses have a double-stranded DNA (dsDNA) genome. To further characterize the diversity of archaeal viruses, we screened highly saline water samples for archaea and their viruses. Here, we describe a new haloarchaeal virus, Halorubrum pleomorphic virus 1 (HRPV-1) that was isolated from a solar saltern and infects an indigenous host belonging to the genus Halorubrum . Infection does not cause cell lysis, but slightly retards growth of the host and results in high replication of the virus. The sequenced genome (7048 nucleotides) of HRPV-1 is single-stranded DNA (ssDNA), which makes HRPV-1 the first characterized archaeal virus that does not have a dsDNA genome. In spite of this, similarities to another archaeal virus were observed. Two major structural proteins were recognized in protein analyses, and by lipid analyses it was shown that the virion contains a membrane. Electron microscopy studies indicate that the enveloped virion is pleomorphic (approximately 44 × 55 nm). HRPV-1 virion may represent commonly used virion architecture, and it seems that structure-based virus lineages may be extended to non-icosahedral viruses.     

Horizontal Transfer of Archaeal Genes into the Deinococcaceae: Detection by Molecular and Computer-Based Approaches

Members of the Deinococcaceae (e.g., Thermus, Meiothermus, Deinococcus) contain A/V-ATPases typically found in Archaea or Eukaryotes which were probably acquired by horizontal gene transfer. Two methods were used to quantify the extent to which archaeal or eukaryotic genes have been acquired by this lineage. Screening of a Meiothermus ruber library with probes made against Thermoplasma acidophilum DNA yielded a number of clones which hybridized more strongly than background. One of these contained the prolyl tRNA synthetase (RS) gene. Phylogenetic analysis shows the M. ruber and D. radiodurans prolyl RS to be more closely related to archaeal and eukaryal forms of this gene than to the typical bacterial type. Using a bioinformatics approach, putative open reading frames (ORFs) from the prerelease version of the D. radiodurans genome were screened for genes more closely related to archaeal or eukaryotic genes. Putative ORFs were searched against representative genomes from each of the three domains using automated BLAST. ORFs showing the highest matches against archaeal and eukaryotic genes were collected and ranked. Among the top-ranked hits were the A/V-ATPase catalytic and noncatalytic subunits and the prolyl RS genes. Using phylogenetic methods, ORFs were analyzed and trees assessed for evidence of horizontal gene transfer. Of the 45 genes examined, 20 showed topologies in which D. radiodurans homologues clearly group with eukaryotic or archaeal homologues, and 17 additional trees were found to show probable evidence of horizontal gene transfer. Compared to the total number of ORFs in the genome, those that can be identified as having been acquired from Archaea or Eukaryotes are relatively few (approximately 1%), suggesting that interdomain transfer is rare.     

Calculating orthologs in bacteria and Archaea: a divide and conquer approach

Among proteins, orthologs are defined as those that are derived by vertical descent from a single progenitor in the last common ancestor of their host organisms. Our goal is to compute a complete set of protein orthologs derived from all currently available complete bacterial and archaeal genomes. Traditional approaches typically rely on all-against-all BLAST searching which is prohibitively expensive in terms of hardware requirements or computational time (requiring an estimated 18 months or more on a typical server). Here, we present xBASE-Orth, a system for ongoing ortholog annotation, which applies a "divide and conquer" approach and adopts a pragmatic scheme that trades accuracy for speed. Starting at species level, xBASE-Orth carefully constructs and uses pan-genomes as proxies for the full collections of coding sequences at each level as it progressively climbs the taxonomic tree using the previously computed data. This leads to a significant decrease in the number of alignments that need to be performed, which translates into faster computation, making ortholog computation possible on a global scale. Using xBASE-Orth, we analyzed an NCBI collection of 1,288 bacterial and 94 archaeal complete genomes with more than 4 million coding sequences in 5 weeks and predicted more than 700 million ortholog pairs, clustered in 175,531 orthologous groups. We have also identified sets of highly conserved bacterial and archaeal orthologs and in so doing have highlighted anomalies in genome annotation and in the proposed composition of the minimal bacterial genome. In summary, our approach allows for scalable and efficient computation of the bacterial and archaeal ortholog annotations. In addition, due to its hierarchical nature, it is suitable for incorporating novel complete genomes and alternative genome annotations. The computed ortholog data and a continuously evolving set of applications based on it are integrated in the xBASE database, available at http://www.xbase.ac.uk/.     

An integrated analysis of the genome of the hyperthermophilic archaeon Pyrococcus abyssi

The hyperthermophilic euryarchaeon Pyrococcus abyssi and the related species Pyrococcus furiosus and Pyrococcus horikoshii, whose genomes have been completely sequenced, are presently used as model organisms in different laboratories to study archaeal DNA replication and gene expression and to develop genetic tools for hyperthermophiles. We have performed an extensive re-annotation of the genome of P. abyssi to obtain an integrated view of its phylogeny, molecular biology and physiology. Many new functions are predicted for both informational and operational proteins. Moreover, several candidate genes have been identified that might encode missing links in key metabolic pathways, some of which have unique biochemical features. The great majority of Pyrococcus proteins are typical archaeal proteins and their phylogenetic pattern agrees with its position near the root of the archaeal tree. However, proteins probably from bacterial origin, including some from mesophilic bacteria, are also present in the P. abyssi genome.     

Remarkable sequence signatures in archaeal genomes

Complete archaeal genomes were probed for the presence of long (> or = 25 bp) oligonucleotide repeats (words). We detected the presence of many words distributed in tandem with narrow ranges of periodicity (i.e., spacer length between repeats). Similar words were not identified in genomes of non-archaeal species, namely Escherichia coli, Bacillus subtilis, Haemophilus influenzae, Mycoplasma genitalium and Mycoplasma pneumoniae. BLAST similarity searches against the GenBank nucleotide sequence database revealed that these words were archaeal species-specific, indicating that they are of a signature character. Sequence analysis and genome viewing tools showed these repeats to be restricted to non-coding regions. Thus, archaea appear to possess a non-coding genomic signature that is absent in bacterial species. The identification of a species-specific genomic signature would be of great value to archaeal genome mapping, evolutionary studies and analyses of genome complexity.     

Archaeal DNA replication and repair

Since the first archaeal genome was sequenced, much attention has been focused on the study of these unique microorganisms. We have learnt that although archaeal DNA metabolic processes (replication, recombination and repair) are more similar to the metabolic processes of Eukarya than those of Bacteria, Archaea are not simply 'mini Eukarya'. They are, in fact, a mosaic of the eukaryal and bacterial systems that also possess archaeal-specific features. Recent biochemical and structural studies of the proteins that participate in archaeal DNA replication and repair have increased our understanding of these processes.     

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

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

Archaeal histone selection of nucleosome positioning sequences and the procaryotic origin of histone-dependent genome evolution

Archaeal histones and the eucaryal (eucaryotic) nucleosome core histones have almost identical histone folds. Here, we show that DNA molecules selectively incorporated by rHMfB (recombinant archaeal histone B from Methanothermus fervidus) into archaeal nucleosomes from a mixture of approximately 10(14) random sequence molecules contain sequence motifs shown previously to direct eucaryal nucleosome positioning. The dinucleotides GC, AA (=TT) and TA are repeated at approximately 10 bp intervals, with the GC harmonic displaced approximately 5 bp from the AA and TA harmonics [(GCN(3)AA or TA)(n)]. AT and CG were not strongly selected, indicating that TA not equalAT and GC not equalCG in terms of facilitating archaeal nucleosome assembly. The selected molecules have affinities for rHMfB ranging from approximately 9 to 18-fold higher than the level of affinity of the starting population, and direct the positioned assembly of archaeal nucleosomes. Fourier-transform analyses have revealed that AA dinucleotides are much enriched at approximately 10. 1 bp intervals, the helical repeat of DNA wrapped around a nucleosome, in the genomes of Eucarya and the histone-containing Euryarchaeota, but not in the genomes of Bacteria and Crenarchaeota, procaryotes that do not have histones. Facilitating histone packaging of genomic DNA has apparently therefore imposed constraints on genome sequence evolution, and since archaeal histones have no structure in addition to the histone fold, these constraints must result predominantly from histone fold-DNA contacts. Based on the three-domain universal phylogeny, histones and histone-dependent genome sequence evolution most likely evolved after the bacterial-archaeal divergence but before the archaeal-eucaryal divergence, and were subsequently lost in the Crenarchaeota. However, with lateral gene transfer, the first histone fold could alternatively have evolved after the archaeal-eucaryal divergence, early in either the euryarchaeal or eucaryal lineages.     

Assembly and function of the archaeal flagellum

Motility is a common behaviour in prokaryotes. Both bacteria and archaea use flagella for swimming motility, but it has been well documented that structures of the flagellum from these two domains of life are completely different, although they contribute to a similar function. Interestingly, information available to date has revealed that structurally archaeal flagella are more similar to bacterial type?IV pili rather than to bacterial flagella. With the increasing genome sequence information and advancement in genetic tools for archaea, identification of the components involved in the assembly of the archaeal flagellum is possible. A subset of these components shows similarities to components from type?IV pilus-assembly systems. Whereas the molecular players involved in assembly of the archaeal flagellum are being identified, the mechanics and dynamics of the assembly of the archaeal flagellum have yet to be established. Recent computational analysis in our laboratory has identified conserved highly charged loop regions within one of the core proteins of the flagellum, the membrane integral protein FlaJ, and predicted that these are involved in the interaction with the assembly ATPase FlaI. Interestingly, considerable variation was found among the loops of FlaJ from the two major subkingdoms of archaea, the Euryarchaeota and the Crenarchaeota. Understanding the assembly pathway and creating an interaction map of the molecular players in the archaeal flagellum will shed light on the details of the assembly and also the evolutionary relationship to the bacterial type?IV pili-assembly systems.     

Identification and distribution of new insertion sequences in the genome of the extremely halotolerant and alkaliphilic Oceanobacillus iheyensis HTE831.

Six kinds of new insertion sequences (ISs), IS667 to IS672, a group II intron (Oi.Int), and an incomplete transposon (Tn852loi) were identified in the 3,630,528-bp genome of the extremely halotolerant and alkaliphilic Oceanobacillus iheyensis HTE831. Of 19 ISs identified in the HTE831 genome, 7 were truncated, indicating the occurrence of internal rearrangement of the genome. All ISs except IS669 generated a 4- to 8-bp duplication of the target site sequence, and these ISs carried 23- to 28-bp inverted repeats (IRs). Sequence analysis revealed that four ISs (IS669, IS670, IS671, and IS672) were newly identified as belonging to separate IS families (IS200/IS605, IS30, IS5, and IS3, respectively). IS667 and IS668 were also characterized as new members of the ISL3 family. Tn8521oi, which belongs to the Tn3 family as a new member, generated a 5-bp duplication of the target site sequence and carried complete 38-bp IRs. Of the eight protein-coding sequences (CDSs) identified in Tn8521oi, three CDSs (OB481, OB482, and OB483) formed a ger gene cluster, and two other paralogous gene clusters were found in the HTE831 genome. Most of the ISs and the group II intron widely distributed throughout the genome were inserted in noncoding regions, while two ISs (IS667-08 and IS668-02) and Oi.Int-04 were inserted in the coding regions.