Understanding the adaptation of Halobacterium species NRC-1 to its extreme environment through computational analysis of its genome sequence

Post-genomics of the model haloarchaeon Halobacterium sp. NRC-1

Halobacteriumsp. NRC-1 is an extremely halophilic archaeon that is easily cultured and genetically tractable. Since its genome sequence was completed in 2000, a combination of genetic, transcriptomic, proteomic, and bioinformatic approaches have provided insights into both its extremophilic lifestyle as well as fundamental cellular processes common to all life forms. Here, we review post-genomic research on this archaeon, including investigations of DNA replication and repair systems, phototrophic, anaerobic, and other physiological capabilities, acidity of the proteome for function at high salinity, and role of lateral gene transfer in its evolution.

Dodecins: a family of lumichrome binding proteins

Dodecin is a small dodecameric flavoprotein from Halobacterium salinarum that contains two flavins stacked between two tryptophan residues to form an aromatic tetrade. The functional properties of heterologously expressed dodecin were investigated by fluorescence spectroscopy, which allowed the determination of dissociation constants for a number of protein-ligand complexes. The values obtained were in the nanomolar to micromolar range and correlate positively with the ligand size. These data were supplemented by X-ray crystal structures of the apododecin and holocomplexes with lumichrome, lumiflavin, riboflavin and FMN at resolutions between 1.55 to 1.95 A to unravel a gating mechanism as the structural basis for the preferential binding of the small ligands lumichrome and lumiflavin. The detailed analysis of the dodecin manifold for preferential binding of lumichrome and lumiflavin provides insight on a subatom level into a protein's strategy to gain selectivity for low molecular mass compounds by steric restrictions rather than specific interactions. Investigations on the ligand composition of a wild-type dodecin crystal (1.32 A resolution) support conclusions of functional and structural investigations on heterologously expressed dodecin, and strongly suggest that lumichrome, a molecule associated with the flavin metabolism, is a ligand of dodecin in vivo. Studies on mutant protein and a Halorhodospira halophila homologue spread the idea of a lumichrome binding system as a possible "waste"-trapping device, widely distributed in prokaryotes.

The genome of Salinibacter ruber: convergence and gene exchange among hyperhalophilic bacteria and archaea

Saturated thalassic brines are among the most physically demanding habitats on Earth: few microbes survive in them. Salinibacter ruber is among these organisms and has been found repeatedly in significant numbers in climax saltern crystallizer communities. The phenotype of this bacterium is remarkably similar to that of the hyperhalophilic Archaea (Haloarchaea). The genome sequence suggests that this resemblance has arisen through convergence at the physiological level (different genes producing similar overall phenotype) and the molecular level (independent mutations yielding similar sequences or structures). Several genes and gene clusters also derive by lateral transfer from (or may have been laterally transferred to) haloarchaea. S. ruber encodes four rhodopsins. One resembles bacterial proteorhodopsins and three are of the haloarchaeal type, previously uncharacterized in a bacterial genome. The impact of these modular adaptive elements on the cell biology and ecology of S. ruber is substantial, affecting salt adaptation, bioenergetics, and photobiology.

Comparison of theoretical proteomes: identification of COGs with conserved and variable pI within the multimodal pI distribution

BACKGROUND

Theoretical proteome analysis, generated by plotting theoretical isoelectric points (pI) against molecular masses of all proteins encoded by the genome show a multimodal distribution for pI. This multimodal distribution is an effect of allowed combinations of the charged amino acids, and not due to evolutionary causes. The variation in this distribution can be correlated to the organisms ecological niche. Contributions to this variation maybe mapped to individual proteins by studying the variation in pI of orthologs across microorganism genomes.

RESULTS

The distribution of ortholog pI values showed trimodal distributions for all prokaryotic genomes analyzed, similar to whole proteome plots. Pairwise analysis of pI variation show that a few COGs are conserved within, but most vary between, the acidic and basic regions of the distribution, while molecular mass is more highly conserved. At the level of functional grouping of orthologs, five groups vary significantly from the population of orthologs, which is attributed to either conservation at the level of sequences or a bias for either positively or negatively charged residues contributing to the function. Individual COGs conserved in both the acidic and basic regions of the trimodal distribution are identified, and orthologs that best represent the variation in levels of the acidic and basic regions are listed.

CONCLUSION

The analysis of pI distribution by using orthologs provides a basis for resolution of theoretical proteome comparison at the level of individual proteins. Orthologs identified that significantly vary between the major acidic and basic regions maybe used as representative of the variation of the entire proteome.

The cobweb of life revealed by genome-scale estimates of horizontal gene transfer

With the availability of increasing amounts of genomic sequences, it is becoming clear that genomes experience horizontal transfer and incorporation of genetic information. However, to what extent such horizontal gene transfer (HGT) affects the core genealogical history of organisms remains controversial. Based on initial analyses of complete genomic sequences, HGT has been suggested to be so widespread that it might be the “essence of phylogeny” and might leave the treelike form of genealogy in doubt. On the other hand, possible biased estimation of HGT extent and the findings of coherent phylogenetic patterns indicate that phylogeny of life is well represented by tree graphs. Here, we reexamine this question by assessing the extent of HGT among core orthologous genes using a novel statistical method based on statistical comparisons of tree topology. We apply the method to 40 microbial genomes in the Clusters of Orthologous Groups database over a curated set of 297 orthologous gene clusters, and we detect significant HGT events in 33 out of 297 clusters over a wide range of functional categories. Estimates of positions of HGT events suggest a low mean genome-specific rate of HGT (2.0%) among the orthologous genes, which is in general agreement with other quantitative of HGT. We propose that HGT events, even when relatively common, still leave the treelike history of phylogenies intact, much like cobwebs hanging from tree branches.

A stastical approach applied to 297 orthologous gene clusters in 40 microbial genomes suggests a low rate of interspecies gene transfer. Species relationships can therefore be modeled with a tree structure.

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.

Analysis of the cytosolic proteome of Halobacterium salinarum and its implication for genome annotation

The halophilic archaeon Halobacterium salinarum (strain R1, DSM 671) contains 2784 protein-coding genes as derived from the genome sequence. The cytosolic proteome containing 2042 proteins was separated by two-dimensional gel electrophoresis (2-DE) and systematically analyzed by a semi-automatic procedure. A reference map was established taking into account the narrow isoelectric point (pI) distribution of halophilic proteins between 3.5 and 5.5. Proteins were separated on overlapping gels covering the essential areas of pI and molecular weight. Every silver-stained spot was analyzed resulting in 661 identified proteins out of about 1800 different protein spots using matrix-assisted laser desorption/ionization time of flight mass spectrometry (MALDI-TOF MS) peptide mass fingerprinting (PMF). There were 94 proteins that were found in multiple spots, indicating post-translational modification. An additional 141 soluble proteins were identified on 2-D gels not corresponding to the reference map. Thus about 40% of the cytosolic proteome was identified. In addition to the 2784 protein-coding genes, the H. salinarum genome contains more than 6000 spurious open reading frames longer than 100 codons. Proteomic information permitted an improvement in genome annotation by validating and correcting gene assignments. The correlation between theoretical pI and gel position is exceedingly good and was used as a tool to improve start codon assignments. The fraction of identified chromosomal proteins was much higher than that of those encoded on the plasmids. In combination with analysis of the GC content this observation permitted an unambiguous identification of an episomal insert of 60 kbp ("AT-rich island") in the chromosome, as well as a 70 kbp region from the chromosome that has integrated into one of the megaplasmids and carries a series of essential genes. About 63% of the chromosomally encoded proteins larger than 25 kDa were identified, proving the efficacy of 2-DE MALDI-TOF MS PMF technology. The analysis of the integral membrane proteome by tandem mass spectrometric techniques added another 141 identified proteins not identified by the 2-DE approach (see following paper).

Crystallization and preliminary X-ray analysis of binary and ternary complexes of Haloferax mediterranei glucose dehydrogenase

Haloferax mediterranei glucose dehydrogenase (EC 1.1.1.47) belongs to the medium-chain alcohol dehydrogenase superfamily and requires zinc for catalysis. In the majority of these family members, the catalytic zinc is tetrahedrally coordinated by the side chains of a cysteine, a histidine, a cysteine or glutamate and a water molecule. In H. mediterranei glucose dehydrogenase, sequence analysis indicates that the zinc coordination is different, with the invariant cysteine replaced by an aspartate residue. In order to analyse the significance of this replacement and to contribute to an understanding of the role of the metal ion in catalysis, a range of binary and ternary complexes of the wild-type and a D38C mutant protein have been crystallized. For most of the complexes, crystals belonging to space group I222 were obtained using sodium/potassium citrate as a precipitant. However, for the binary and non-productive ternary complexes with NADPH/Zn, it was necessary to replace the citrate with 2-methyl-2,4-pentanediol. Despite the radical change in conditions, the crystals thus formed were isomorphous.

Identification of replication origins in archaeal genomes based on the Z-curve method

The Z-curve is a three-dimensional curve that constitutes a unique representation of a DNA sequence, i.e., both the Z-curve and the given DNA sequence can be uniquely reconstructed from the other. We employed Z-curve analysis to identify one replication origin in the Methanocaldococcus jannaschii genome, two replication origins in the Halobacterium species NRC-1 genome and one replication origin in the Methanosarcina mazei genome. One of the predicted replication origins of Halobacterium species NRC-1 is the same as a replication origin later identified by in vivo experiments. The Z-curve analysis of the Sulfolobus solfataricus P2 genome suggested the existence of three replication origins, which is also consistent with later experimental results. This review aims to summarize applications of the Z-curve in identifying replication origins of archaeal genomes, and to provide clues about the locations of as yet unidentified replication origins of the Aeropyrum pernix K1, Methanococcus maripaludis S2, Picrophilus torridus DSM 9790 and Pyrobaculum aerophilum str. IM2 genomes.

Weighted genome trees: refinements and applications

There are many ways to group completed genome sequences in hierarchical patterns (trees) reflecting relationships between their genes. Such groupings help us organize biological information and bear crucially on underlying processes of genome and organismal evolution. Genome trees make use of all comparable genes but can variously weight the contributions of these genes according to similarity, congruent patterns of similarity, or prevalence among genomes. Here we explore such possible weighting strategies, in an analysis of 142 prokaryotic and 5 eukaryotic genomes. We demonstrate that alternate weighting strategies have different advantages, and we propose that each may have its specific uses in systematic or evolutionary biology. Comparisons of results obtained with different methods can provide further clues to major events and processes in genome evolution.

An efficient proteomics based strategy for the functional characterization of a novel halophilic enzyme fromHalobacterium salinarum

The extremely halophilic archaeon, Halobacterium salinarum grows in environments containing over 25% NaCl. The enzymes of this organism have thus been adapted to be active and stable in hypersaline conditions, which makes them strong candidates as robust industrial enzymes. In this study, the proteomics approach was applied to screen novel halophilic enzymes. We focused initially on proteins that are differentially expressed under different salt concentrations in culture media. After two-dimensional gel electrophoresis over a pH 3.5-4.5 range, 29 differentially expressed protein spots were identified by tandem mass spectrometry and six of these had no similarity to preexisting genes of known function. To predict the function of them, we used various bioinformatic methods. Among other proteins, we selected Vng0487h, which showed a high similarity to acetyltransferases. As a step toward assaying the enzymatic activity of this protein, we cloned the Vng0487h gene of H. salinarum and expressed and purified the recombinant protein with a glutathione-S-transferase (GST) tag in Escherichia coli. Using a GST-pulldown assay, a protein fragment derived from E. coli could interact with recombinant Vng0487h, and was identified to be the ribosomal protein L3. This protein showed high sequence homology with ribosomal protein L7/12 from E. coli and ribosomal protein L13p from H. salinarum. This suggests that Vng0487h acetylates a subunit of ribosomal protein, possibly L13p, in H. salinarum. During the present study, an efficient procedure was established to screen novel halophilic enzymes, and to predict and assess their functions.

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.

Genome sequence of Haloarcula marismortui: a halophilic archaeon from the Dead Sea

We report the complete sequence of the 4,274,642-bp genome of Haloarcula marismortui, a halophilic archaeal isolate from the Dead Sea. The genome is organized into nine circular replicons of varying G+C compositions ranging from 54% to 62%. Comparison of the genome architectures of Halobacterium sp. NRC-1 and H. marismortui suggests a common ancestor for the two organisms and a genome of significantly reduced size in the former. Both of these halophilic archaea use the same strategy of high surface negative charge of folded proteins as means to circumvent the salting-out phenomenon in a hypersaline cytoplasm. A multitiered annotation approach, including primary sequence similarities, protein family signatures, structure prediction, and a protein function association network, has assigned putative functions for at least 58% of the 4242 predicted proteins, a far larger number than is usually achieved in most newly sequenced microorganisms. Among these assigned functions were genes encoding six opsins, 19 MCP and/or HAMP domain signal transducers, and an unusually large number of environmental response regulators-nearly five times as many as those encoded in Halobacterium sp. NRC-1--suggesting H. marismortui is significantly more physiologically capable of exploiting diverse environments. In comparing the physiologies of the two halophilic archaea, in addition to the expected extensive similarity, we discovered several differences in their metabolic strategies and physiological responses such as distinct pathways for arginine breakdown in each halophile. Finally, as expected from the larger genome, H. marismortui encodes many more functions and seems to have fewer nutritional requirements for survival than does Halobacterium sp. NRC-1.

Identification and characterization of inosine monophosphate dehydrogenase fromHalobacterium salinarum

Extremely halophilic Archaea, Halobacterium salinarum live in hypersaline habitats and maintain an osmotic balance of their cytoplasm by accumulating high concentrations of salt (mainly KCl). Therefore, their enzymes adapted to high NaCl concentrations offer a multitude of acutal or potential applications such as biocatalysts in the presence of high salt concentrations. In this study, the protein expression profile of H. salinarum cultured under different NaCl concentrations (3.5 M, 4.3 M, and 6.0 M) was investigated using two-dimensional gel electrophoresis (2-DE). As a result of 2-DE, the protein spots concentrated in acidic range at pH 3-10 were separated effectively using pH 3.5-4.5 ultrazoom IPG DryStrips. The proteins which proved to be upregulated or downregulated in 2-DE gel were digested with trypsin and identified with matrix-assisted laser desorption/ionization-time of flight (MALDI-TOF) and electrospray ionization quadrupole (ESI-Q) TOF-mass spectrometry. Most proteins were identified as known annotated proteins based on sequence homology and few as unknown hypothetical proteins. Among proteins identified, an enzyme named inosine monophosphate dehydrogenase (IMPDH) was selected based on the possibility of its industrial application. IMPDH gene (1.6 kb fragment) expected to exist in H. salinarum was amplified by polymerase chain reaction (PCR) and expressed in Escherichia coli strain, BL21 (DE3) using a pGEX-KG vector. Recombinant IMPDH purified from H. salinarum has a higher activity in the presence of salt than in the absence of salt.

Cassette-based presentation of SIV epitopes with recombinant gas vesicles from halophilic archaea

In earlier studies we demonstrated recombinant gas vesicles from Halobacterium sp. NRC-1, expressing a model six amino acid insert, or native vesicles displaying chemically coupled TNP, each were immunogenic, and antigenic. Long-lived responses displaying immunologic memory were elicited without exogenous adjuvant. Here we report the generation and expression of cassettes containing SIV derived DNA. The results indicate a cassette-based display/delivery system derived from recombinant halobacterial gas vesicle genes is highly feasible. Data specifically support four conclusions: (i) Recombinants carrying up to 705 bp of SIV DNA inserted into the gvpC gene form functional gas vesicles; (ii) SIV peptides contained as part of the expressed recombinant, surface exposed GvpC protein are recognized by antibody elicited in monkeys exposed to native SIV in vivo; (iii) in the absence of adjuvant, mice immunized with the recombinant gas vesicle (r-GV) preparations mount a solid, titratable antibody response to the test SIV insert that is long lived and exhibits immunologic memory; (iv) recombinant organelles, created through the generation of cassettes encoding epitopes inserted into the gvpC DNA, can be used to construct a multiepitope display (MED) library, a potentially cost effective vehicle to express and deliver peptides of SIV, HIV or other pathogens.

Comprehensive de novo structure prediction in a systems-biology context for the archaea Halobacterium sp. NRC-1

BACKGROUND

Large fractions of all fully sequenced genomes code for proteins of unknown function. Annotating these proteins of unknown function remains a critical bottleneck for systems biology and is crucial to understanding the biological relevance of genome-wide changes in mRNA and protein expression, protein-protein and protein-DNA interactions. The work reported here demonstrates that de novo structure prediction is now a viable option for providing general function information for many proteins of unknown function.

RESULTS

We have used Rosetta de novo structure prediction to predict three-dimensional structures for 1,185 proteins and protein domains (<150 residues in length) found in Halobacterium NRC-1, a widely studied halophilic archaeon. Predicted structures were searched against the Protein Data Bank to identify fold similarities and extrapolate putative functions. They were analyzed in the context of a predicted association network composed of several sources of functional associations such as: predicted protein interactions, predicted operons, phylogenetic profile similarity and domain fusion. To illustrate this approach, we highlight three cases where our combined procedure has provided novel insights into our understanding of chemotaxis, possible prophage remnants in Halobacterium NRC-1 and archaeal transcriptional regulators.

CONCLUSIONS

Simultaneous analysis of the association network, coordinated mRNA level changes in microarray experiments and genome-wide structure prediction has allowed us to glean significant biological insights into the roles of several Halobacterium NRC-1 proteins of previously unknown function, and significantly reduce the number of proteins encoded in the genome of this haloarchaeon for which no annotation is available.

Global analysis of predicted proteomes: functional adaptation of physical properties

The physical characteristics of proteins are fundamentally important in organismal function. We used the complete predicted proteomes of >100 organisms spanning the three domains of life to investigate the comparative biology and evolution of proteomes. Theoretical 2D gels were constructed with axes of protein mass and charge (pI) and converted to density estimates comparable across all types and sizes of proteome. We asked whether we could detect general patterns of proteome conservation and variation. The overall pattern of theoretical 2D gels was strongly conserved across all life forms. Nevertheless, coevolved replicons from the same organism (different chromosomes or plasmid and host chromosomes) encode proteomes more similar to each other than those from different organisms. Furthermore, there was disparity between the membrane and nonmembrane subproteomes within organisms (proteins of membrane proteomes are on the average more basic and heavier) and their variation across organisms, suggesting that membrane proteomes evolve most rapidly. Experimentally, a significant positive relationship independent of phylogeny was found between the predicted proteome and Biolog profile, a measure associated with the ecological niche. Finally, we show that, for the smallest and most alkaline proteomes, there is a negative relationship between proteome size and basicity. This relationship is not adequately explained by AT bias at the DNA sequence level. Together, these data provide evidence of functional adaptation in the properties of complete proteomes.

Systems level insights into the stress response to UV radiation in the halophilic archaeon Halobacterium NRC-1

We report a remarkably high UV-radiation resistance in the extremely halophilic archaeon Halobacterium NRC-1 withstanding up to 110 J/m2 with no loss of viability. Gene knockout analysis in two putative photolyase-like genes (phr1 and phr2) implicated only phr2 in photoreactivation. The UV-response was further characterized by analyzing simultaneously, along with gene function and protein interactions inferred through comparative genomics approaches, mRNA changes for all 2400 genes during light and dark repair. In addition to photoreactivation, three other putative repair mechanisms were identified including d(CTAG) methylation-directed mismatch repair, four oxidative damage repair enzymes, and two proteases for eliminating damaged proteins. Moreover, a UV-induced down-regulation of many important metabolic functions was observed during light repair and seems to be a phenomenon shared by all three domains of life. The systems analysis has facilitated the assignment of putative functions to 26 of 33 key proteins in the UV response through sequence-based methods and/or similarities of their predicted three-dimensional structures to known structures in the PDB. Finally, the systems analysis has raised, through the integration of experimentally determined and computationally inferred data, many experimentally testable hypotheses that describe the metabolic and regulatory networks of Halobacterium NRC-1.

What we have learned about prokaryotes from structural genomics

Five years ago systematic determination and theoretical analysis of all protein structures encoded in model prokaryotic organisms was proposed as a powerful way to obtain new insights into protein function and the variety of protein folds. What has been the pay-off from studying structures in genomic context? Have we learned anything new about protein structure? Can we now predict protein function better? In this contribution, I summarize the status of large-scale structure determination projects on prokaryotes and provide an overview of the main results obtained from experimental and theoretical studies in this dynamic research field.

Archaeal phylogeny based on proteins of the transcription and translation machineries: tackling the Methanopyrus kandleri paradox

This article presents a phylogenetic analysis of the Archea based on sets of transcription and translation proteins. The phylogenies shed light on the evolutionary position of Methanopyrus kandleri.

Background

Phylogenetic analysis of the Archaea has been mainly established by 16S rRNA sequence comparison. With the accumulation of completely sequenced genomes, it is now possible to test alternative approaches by using large sequence datasets. We analyzed archaeal phylogeny using two concatenated datasets consisting of 14 proteins involved in transcription and 53 ribosomal proteins (3,275 and 6,377 positions, respectively).

Results

Important relationships were confirmed, notably the dichotomy of the archaeal domain as represented by the Crenarchaeota and Euryarchaeota, the sister grouping of Sulfolobales and Aeropyrum pernix, and the monophyly of a large group comprising Thermoplasmatales, Archaeoglobus fulgidus, Methanosarcinales and Halobacteriales, with the latter two orders forming a robust cluster. The main difference concerned the position of Methanopyrus kandleri, which grouped with Methanococcales and Methanobacteriales in the translation tree, whereas it emerged at the base of the euryarchaeotes in the transcription tree. The incongruent placement of M. kandleri is likely to be the result of a reconstruction artifact due to the high evolutionary rates displayed by the components of its transcription apparatus.

Conclusions

We show that two informational systems, transcription and translation, provide a largely congruent signal for archaeal phylogeny. In particular, our analyses support the appearance of methanogenesis after the divergence of the Thermococcales and a late emergence of aerobic respiration from within methanogenic ancestors. We discuss the possible link between the evolutionary acceleration of the transcription machinery in M. kandleri and several unique features of this archaeon, in particular the absence of the elongation transcription factor TFS.

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.

Archeal DNA replication: eukaryal proteins in a bacterial context

Genome sequences of a number of archaea have revealed an apparent paradox in the phylogenies of the bacteria, archaea, and eukarya, as well as an intriguing set of problems to be resolved in the study of DNA replication. The archaea, long thought to be bacteria, are not only different enough to merit their own domain but also appear to be an interesting mosaic of bacterial, eukaryal, and unique features. Most archaeal proteins participating in DNA replication are more similar in sequence to those found in eukarya than to analogous replication proteins in bacteria. However, archaea have only a subset of the eukaryal replication machinery, apparently needing fewer polypeptides and structurally simpler complexes. The archaeal replication apparatus also contains features not found in other organisms owing, in part, to the broad range of environmental conditions, some extreme, in which members of this domain thrive. In this review the current knowledge of the mechanisms governing DNA replication in archaea is summarized and the similarities and differences of those of bacteria and eukarya are highlighted.

Low-pass sequencing for microbial comparative genomics

BACKGROUND

We studied four extremely halophilic archaea by low-pass shotgun sequencing: (1) the metabolically versatile Haloarcula marismortui; (2) the non-pigmented Natrialba asiatica; (3) the psychrophile Halorubrum lacusprofundi and (4) the Dead Sea isolate Halobaculum gomorrense. Approximately one thousand single pass genomic sequences per genome were obtained. The data were analyzed by comparative genomic analyses using the completed Halobacterium sp. NRC-1 genome as a reference. Low-pass shotgun sequencing is a simple, inexpensive, and rapid approach that can readily be performed on any cultured microbe.

RESULTS

As expected, the four archaeal halophiles analyzed exhibit both bacterial and eukaryotic characteristics as well as uniquely archaeal traits. All five halophiles exhibit greater than sixty percent GC content and low isoelectric points (pI) for their predicted proteins. Multiple insertion sequence (IS) elements, often involved in genome rearrangements, were identified in H. lacusprofundi and H. marismortui. The core biological functions that govern cellular and genetic mechanisms of H. sp. NRC-1 appear to be conserved in these four other halophiles. Multiple TATA box binding protein (TBP) and transcription factor IIB (TFB) homologs were identified from most of the four shotgunned halophiles. The reconstructed molecular tree of all five halophiles shows a large divergence between these species, but with the closest relationship being between H. sp. NRC-1 and H. lacusprofundi.

CONCLUSION

Despite the diverse habitats of these species, all five halophiles share (1) high GC content and (2) low protein isoelectric points, which are characteristics associated with environmental exposure to UV radiation and hypersalinity, respectively. Identification of multiple IS elements in the genome of H. lacusprofundi and H. marismortui suggest that genome structure and dynamic genome reorganization might be similar to that previously observed in the IS-element rich genome of H. sp. NRC-1. Identification of multiple TBP and TFB homologs in these four halophiles are consistent with the hypothesis that different types of complex transcriptional regulation may occur through multiple TBP-TFB combinations in response to rapidly changing environmental conditions. Low-pass shotgun sequence analyses of genomes permit extensive and diverse analyses, and should be generally useful for comparative microbial genomics.

An archaeal chromosomal autonomously replicating sequence element from an extreme halophile, Halobacterium sp. strain NRC-1

We report on the identification and first cloning of an autonomously replicating sequence element from the chromosome of an archaeon, the extreme halophile Halobacterium strain NRC-1. The putative replication origin was identified by association with the orc7 gene and replication ability in the host strain, demonstrated by cloning into a nonreplicating plasmid. Deletion analysis showed that sequences located up to 750 bp upstream of the orc7 gene translational start, plus the orc7 gene and 50 bp downstream, are sufficient to endow the plasmid with replication ability, as judged by expression of a plasmid-encoded mevinolin resistance selectable marker and plasmid recovery after transformation. Sequences located proximal to the two other chromosomally carried haloarchaeal orc genes (orc6 and orc8) are not able to promote efficient autonomous replication. Located within the 750-bp region upstream of orc7 is a nearly perfect inverted repeat of 31 bp, which flanks an extremely AT-rich (44%) stretch of 189 bp. The replication ability of the plasmid was lost when one copy of the inverted repeat was deleted. Additionally, the inverted repeat structure near orc7 homologs in the genomic sequences of two other halophiles, Haloarcula marismortui and Haloferax volcanii, is highly conserved. Our results indicate that, in halophilic archaea, a chromosomal origin of replication is physically linked to orc7 homologs and that this element is sufficient to promote autonomous replication. We discuss the finding of a functional haloarchaeal origin in relation to the large number of orc1-cdc6 homologs identified in the genomes of all haloarchaea to date.

Purification and biological characterization of halocin C8, a novel peptide antibiotic from Halobacterium strain AS7092

Halocins are bacteriocin-like proteins or peptides produced by many species of the family Halobacteriaceae. Halocin C8, excreted by the Halobacterium strain AS7092, is a single 6.3-kDa polypeptide with an isoelectric point of 4.4, which is sensitive to proteinase K but not to trypsin. Halocin C8 is quite stable, as it can be desalted, boiled, frozen, subjected to organic solvents, and stored in culture supernatant at 4 degrees C or in dH(2)O at -20 degrees C for more than 1 year without losing activity. The purification of this halocin was achieved by combination of tangential flow filtration (TFF), Sephadex G50 and DEAE-sepharose chromatography. The N-terminal amino acid sequence was also determined by Edman degradation. Halocin C8 appeared to have a very wide activity spectrum, including most haloarchaea and even some haloalkaliphilic rods. When a sensitive strain of Halorubrum saccharovorum was exposed to halocin C8, the treated cells swelled at the initial stage, the cell wall appeared to be nicked and the cytoplasm was then extruded out, and the whole cell was eventually completely lysed. These results indicate that halocin C8 is a novel microhalocin and its primary target might be located in the cell wall of the sensitive cells.

Mechanisms of thermal adaptation revealed from the genomes of the Antarctic Archaea Methanogenium frigidum and Methanococcoides burtonii

We generated draft genome sequences for two cold-adapted Archaea, Methanogenium frigidum and Methanococcoides burtonii, to identify genotypic characteristics that distinguish them from Archaea with a higher optimal growth temperature (OGT). Comparative genomics revealed trends in amino acid and tRNA composition, and structural features of proteins. Proteins from the cold-adapted Archaea are characterized by a higher content of noncharged polar amino acids, particularly Gln and Thr and a lower content of hydrophobic amino acids, particularly Leu. Sequence data from nine methanogen genomes (OGT 15 degrees -98 degrees C) were used to generate 1111 modeled protein structures. Analysis of the models from the cold-adapted Archaea showed a strong tendency in the solvent-accessible area for more Gln, Thr, and hydrophobic residues and fewer charged residues. A cold shock domain (CSD) protein (CspA homolog) was identified in M. frigidum, two hypothetical proteins with CSD-folds in M. burtonii, and a unique winged helix DNA-binding domain protein in M. burtonii. This suggests that these types of nucleic acid binding proteins have a critical role in cold-adapted Archaea. Structural analysis of tRNA sequences from the Archaea indicated that GC content is the major factor influencing tRNA stability in hyperthermophiles, but not in the psychrophiles, mesophiles or moderate thermophiles. Below an OGT of 60 degrees C, the GC content in tRNA was largely unchanged, indicating that any requirement for flexibility of tRNA in psychrophiles is mediated by other means. This is the first time that comparisons have been performed with genome data from Archaea spanning the growth temperature extremes from psychrophiles to hyperthermophiles.

Aminoacylation of an unusual tRNA(Cys) from an extreme halophile

The extreme halophile Halobacterium species NRC-1 overcomes external near-saturating salt concentrations by accumulating intracellular salts comparable to those of the medium. This raises the fundamental question of how halophiles can maintain the specificity of protein-nucleic acid interactions that are particularly sensitive to high salts in mesophiles. Here we address the specificity of the essential aminoacylation reaction of the halophile, by focusing on molecular recognition of tRNA(Cys) by the cognate cysteinyl-tRNA synthetase. Despite the high salt environments of the aminoacylation reaction, and despite an unusual structure of the tRNA with an exceptionally large dihydrouridine loop, we show that aminoacylation of the tRNA proceeds with a catalytic efficiency similar to that of its mesophilic counterparts. This is manifested by an essentially identical K(m) for tRNA to those of the mesophiles, and by recognition of the same nucleotide determinants that are conserved in evolution. Interestingly, aminoacylation of the halophile tRNA(Cys) is more closely related to that of bacteria than eukarya by placing a strong emphasis on features of the tRNA tertiary core. This suggests an adaptation to the highly negatively charged tRNA sugar-phosphate backbone groups that are the key elements of the tertiary core.

Archaea: an archetype for replication initiation studies?

Whereas the process of DNA replication is fundamentally conserved in the three domains of life, the archaeal system is closer to that of eukarya than bacteria. In the time since the complete genome sequences of several members of the archaeal domain became available, there has been a burst of research on archaeal DNA replication. These studies have led to both expected and surprising findings. This review summarizes the search for origins of replication in archaea, and our current knowledge of initiation, the process by which replication origins are recognized, the DNA molecule is unwound and the replicative helicase is loaded onto the DNA in preparation for DNA synthesis. The similarities and differences of the initiation process in archea, bacteria and eukarya are also summarized.

Unique amino acid composition of proteins in halophilic bacteria

The amino acid compositions of proteins from halophilic archaea were compared with those from non-halophilic mesophiles and thermophiles, in terms of the protein surface and interior, on a genome-wide scale. As we previously reported for proteins from thermophiles, a biased amino acid composition also exists in halophiles, in which an abundance of acidic residues was found on the protein surface as compared to the interior. This general feature did not seem to depend on the individual protein structures, but was applicable to all proteins encoded within the entire genome. Unique protein surface compositions are common in both halophiles and thermophiles. Statistical tests have shown that significant surface compositional differences exist among halophiles, non-halophiles, and thermophiles, while the interior composition within each of the three types of organisms does not significantly differ. Although thermophilic proteins have an almost equal abundance of both acidic and basic residues, a large excess of acidic residues in halophilic proteins seems to be compensated by fewer basic residues. Aspartic acid, lysine, asparagine, alanine, and threonine significantly contributed to the compositional differences of halophiles from meso- and thermophiles. Among them, however, only aspartic acid deviated largely from the expected amount estimated from the dinucleotide composition of the genomic DNA sequence of the halophile, which has an extremely high G+C content (68%). Thus, the other residues with large deviations (Lys, Ala, etc.) from their non-halophilic frequencies could have arisen merely as "dragging effects" caused by the compositional shift of the DNA, which would have changed to increase principally the fraction of aspartic acid alone.

Multiple replication origins of the archaeon Halobacterium species NRC-1

The genomic sequence of the halophilic archaeon Halobacterium NRC-1 has been analyzed by the Z curve method. The Z curve is a three-dimensional curve that uniquely represents a given DNA sequence. Based on the known behaviors of the Z curves for the archaea whose replication origins have been identified, the analysis of the Z curve for the genome of Halobacterium NRC-1 strongly suggests that the large genome has two replication origins, oriC1 (921,863-922,014) and oriC2 (1,806,444-1,807,229), which are located at two sharp peaks of the Z curve. These two regions are next to the cdc6 genes and contain multiple copies of stretches of G and C, i.e., ggggtgggg and ccccacccc, which may also be regarded as direct and inverted repeats. Based on the above analysis, a model of replication of Halobacterium NRC-1 with two replication origins and two termini has been proposed. The experimental confirmation of this model would constitute the first example of multiple replication origins of archaea, which will finally provide much insight into the understanding of replication mechanisms of eukaryotic organisms, including human. In addition, the potential multiple replication origins of the archaeon Sulfolobus solfataricus are suggested by the analysis based on the Z curve method.

Simple sequence repeats and compositional bias in the bipartite Ralstonia solanacearum GMI1000 genome

BACKGROUND

Ralstonia solanacearum is an important plant pathogen. The genome of R. solananearum GMI1000 is organised into two replicons (a 3.7-Mb chromosome and a 2.1-Mb megaplasmid) and this bipartite genome structure is characteristic for most R. solanacearum strains. To determine whether the megaplasmid was acquired via recent horizontal gene transfer or is part of an ancestral single chromosome, we compared the abundance, distribution and composition of simple sequence repeats (SSRs) between both replicons and also compared the respective compositional biases.

RESULTS

Our data show that both replicons are very similar in respect to distribution and composition of SSRs and presence of compositional biases. Minor variations in SSR and compositional biases observed may be attributable to minor differences in gene expression and regulation of gene expression or can be attributed to the small sample numbers observed.

CONCLUSIONS

The observed similarities indicate that both replicons have shared a similar evolutionary history and thus suggest that the megaplasmid was not recently acquired from other organisms by lateral gene transfer but is a part of an ancestral R. solanacearum chromosome.

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

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

Structural plasticity of thermophilic serine hydroxymethyltransferases

Serine hydroxymethyltransferase (SHMT) catalyzes the reversible cleavage of serine to form glycine and monocarbonic groups, essential in several biosynthetic pathways. The availability of crystallographic structures of SHMT from mesophilic organisms and information produced by the genomic projects prompted the analysis of the adaptation of SHMT to "extreme" environments, such as high temperatures, by exploitation of structural data from thermophilic organisms. The sequences of 10 thermophilic/hyperthermophilic SHMTs were multiply aligned to 53 mesophilic homologs and analyzed by a comparative approach, examining the amino acid compositions and preferred residue exchanges between mesophiles and extremophiles. The structural basis of the observed exchanges was further investigated through the application of homology modeling to the 10 extremophilic SHMTs. The results of this study indicate that, in SHMT, thermal stability can be achieved mainly through three strategies: (i) increased number of charged residues at the protein surface; (ii) increased hydrophobicity of the protein core; and (iii) substitution of thermolabile residues exposed to the solvent. Additional features of the archaeal SHMTs, for which no structural data are available yet, were also investigated to explain their quaternary assemblage and the interaction with modified folates.

Genomic BLAST: custom-defined virtual databases for complete and unfinished genomes

BLAST (Basic Local Alignment Search Tool) searches against DNA and protein sequence databases have become an indispensable tool for biomedical research. The proliferation of the genome sequencing projects is steadily increasing the fraction of genome-derived sequences in the public databases and their importance as a public resource. We report here the availability of Genomic BLAST, a novel graphical tool for simplifying BLAST searches against complete and unfinished genome sequences. This tool allows the user to compare the query sequence against a virtual database of DNA and/or protein sequences from a selected group of organisms with finished or unfinished genomes. The organisms for such a database can be selected using either a graphic taxonomy-based tree or an alphabetical list of organism-specific sequences. The first option is designed to help explore the evolutionary relationships among organisms within a certain taxonomy group when performing BLAST searches. The use of an alphabetical list allows the user to perform a more elaborate set of selections, assembling any given number of organism-specific databases from unfinished or complete genomes. This tool, available at the NCBI web site http://www.ncbi.nlm.nih.gov/cgi-bin/Entrez/genom_table_cgi, currently provides access to over 170 bacterial and archaeal genomes and over 40 eukaryotic genomes.

Chromosome replication patterns in the hyperthermophilic euryarchaea Archaeoglobus fulgidus and Methanocaldococcus (Methanococcus) jannaschii

We analysed chromosome replication patterns in the two hyperthermophilic euryarchaea Archaeoglobus fulgidus and Methanocaldococcus(Methanococcus) jannaschii by marker frequency analysis (MFA). For A. fulgidus, the central region of the chromosomal physical map displayed a higher relative abundance in gene dosage during exponential growth, with two continuous gradients to a region of lower abundance at the diametrically opposite side of the genome map. This suggests bidirectional replication of the A. fulgidus chromosome from a single origin. The organization of the putative replication origin region relative to the cdc6, mcm and DNA polymerase genes differed from that reported for Pyrococcus species. No single replication origin or termination regions could be identified for M. jannaschii, adding to the list of unusual properties of this organism. The organization of the A. fulgidus cell cycle was characterized by flow cytometry analysis of the samples from which genomic DNA was extracted for MFA. The relative lengths of the cell cycle periods were found to be similar to those of crenarchaea.

The sequence of the major gas vesicle protein, GvpA, influences the width and strength of halobacterial gas vesicles

Transformation experiments with Haloferax volcanii show that the amino acid sequence of the gas vesicle protein GvpA influences the morphology and strength of gas vesicles produced by halophilic archaea. A modified expression vector containing p-gvpA was used to complement a Vac(-) strain of Hfx. volcanii that harboured the entire p-vac region (from Halobacterium salinarum PHH1) except for p-gvpA. Replacement of p-gvpA with mc-gvpA (from Haloferax mediterranei) led to the synthesis of gas vesicles that were narrower and stronger. Other gene replacements (using c-gvpA from Hbt. salinarum or mutated p-gvpA sequences) led to a significant but smaller increase in gas vesicle strength, and less marked effects on gas vesicle morphology.

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.

A draft sequence of the rice genome (Oryza sativa L. ssp. indica)

We have produced a draft sequence of the rice genome for the most widely cultivated subspecies in China, Oryza sativa L. ssp. indica, by whole-genome shotgun sequencing. The genome was 466 megabases in size, with an estimated 46,022 to 55,615 genes. Functional coverage in the assembled sequences was 92.0%. About 42.2% of the genome was in exact 20-nucleotide oligomer repeats, and most of the transposons were in the intergenic regions between genes. Although 80.6% of predicted Arabidopsis thaliana genes had a homolog in rice, only 49.4% of predicted rice genes had a homolog in A. thaliana. The large proportion of rice genes with no recognizable homologs is due to a gradient in the GC content of rice coding sequences.

Extremophiles and the search for extraterrestrial life

Extremophiles thrive in ice, boiling water, acid, the water core of nuclear reactors, salt crystals, and toxic waste and in a range of other extreme habitats that were previously thought to be inhospitable for life. Extremophiles include representatives of all three domains (Bacteria, Archaea, and Eucarya); however, the majority are microorganisms, and a high proportion of these are Archaea. Knowledge of extremophile habitats is expanding the number and types of extraterrestrial locations that may be targeted for exploration. In addition, contemporary biological studies are being fueled by the increasing availability of genome sequences and associated functional studies of extremophiles. This is leading to the identification of new biomarkers, an accurate assessment of cellular evolution, insight into the ability of microorganisms to survive in meteorites and during periods of global extinction, and knowledge of how to process and examine environmental samples to detect viable life forms. This paper evaluates extremophiles and extreme environments in the context of astrobiology and the search for extraterrestrial life.

Genomic perspective on the photobiology of Halobacterium species NRC-1, a phototrophic, phototactic, and UV-tolerant haloarchaeon

Halobacterium species display a variety of responses to light, including phototrophic growth, phototactic behavior, and photoprotective mechanisms. The complete genome sequence of Halobacterium species NRC-1 (Proc Natl Acad Sci USA 97: 12176-12181, 2000), coupled with the availability of a battery of methods for its analysis makes this an ideal model system for studying photobiology among the archaea. Here, we review: (1) the structure of the 2.57 Mbp Halobacterium NRC-1 genome, including a large chromosome, two minichromosomes, and 91 transposable IS elements; (2) the purple membrane regulon, which programs the accumulation of large quantities of the light-driven proton pump, bacteriorhodopsin, and allows for a period of phototrophic growth; (3) components of the sophisticated pathways for color-sensitive phototaxis; (4) the gas vesicle gene cluster, which codes for cell buoyancy organelles; (5) pathways for the production of carotenoid pigments and retinal, (6) processes for the repair of DNA damage; and (7) putative homologs of circadian rhythm regulators. We conclude with a discussion of the power of systems biology for comprehensive understanding of Halobacterium NRC-1 photobiology.