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

Quantitative analysis of replication-related mutation and selection pressures in bacterial chromosomes and plasmids using generalised GC skew index

BACKGROUND

Due to their bi-directional replication machinery starting from a single finite origin, bacterial genomes show characteristic nucleotide compositional bias between the two replichores, which can be visualised through GC skew or (C-G)/(C+G). Although this polarisation is used for computational prediction of replication origins in many bacterial genomes, the degree of GC skew visibility varies widely among different species, necessitating a quantitative measurement of GC skew strength in order to provide confidence measures for GC skew-based predictions of replication origins.

RESULTS

Here we discuss a quantitative index for the measurement of GC skew strength, named the generalised GC skew index (gGCSI), which is applicable to genomes of any length, including bacterial chromosomes and plasmids. We demonstrate that gGCSI is independent of the window size and can thus be used to compare genomes with different sizes, such as bacterial chromosomes and plasmids. It can suggest the existence of different replication mechanisms in archaea and of rolling-circle replication in plasmids. Correlation of gGCSI values between plasmids and their corresponding host chromosomes suggests that within the same strain, these replicons have reproduced using the same replication machinery and thus exhibit similar strengths of replication strand skew.

CONCLUSIONS

gGCSI can be applied to genomes of any length and thus allows comparative study of replication-related mutation and selection pressures in genomes of different lengths such as bacterial chromosomes and plasmids. Using gGCSI, we showed that replication-related mutation or selection pressure is similar for replicons with similar machinery.

Halophilic beta-lactamase as a new solubility- and folding-enhancing tag protein: production of native human interleukin 1alpha and human neutrophil alpha-defensin

The amino acid composition of halophilic enzymes is characterized by an abundant content of acidic amino acid, which confers to the halophilic enzymes extensive negative charges at neutral pH and high aqueous solubility. This negative charge prevents protein aggregation when denatured and thereby leads to highly efficient protein refolding. Beta-lactamase from periplasmic space of moderate halophile (BLA), a typical halophilic enzyme, can be readily expressed as a native, active form in Escherichia coli cytoplasm. Similar to other halophilic enzymes, BLA is soluble upon denaturation by heat or urea treatments and, hence, can be efficiently refolded. Such high solubility and refolding efficiency make BLA a potential fusion partner for expression of aggregation-prone heterologous proteins to be expressed in E. coli. Here, we succeeded in the soluble expression of several "difficult-to-express" proteins as a BLA fusion protein and verified biological activities of human interleukin 1alpha and human neutrophil alpha-defensin, HNP-1.

Getting a better picture of microbial evolution en route to a network of genomes

Most current thinking about evolution is couched in the concept of trees. The notion of a tree with recursively bifurcating branches representing recurrent divergence events is a plausible metaphor to describe the evolution of multicellular organisms like vertebrates or land plants. But if we try to force the tree metaphor onto the whole of the evolutionary process, things go badly awry, because the more closely we inspect microbial genomes through the looking glass of gene and genome sequence comparisons, the smaller the amount of the data that fits the concept of a bifurcating tree becomes. That is mainly because among microbes, endosymbiosis and lateral gene transfer are important, two mechanisms of natural variation that differ from the kind of natural variation that Darwin had in mind. For such reasons, when it comes to discussing the relationships among all living things, that is, including the microbes and all of their genes rather than just one or a select few, many biologists are now beginning to talk about networks rather than trees in the context of evolutionary relationships among microbial chromosomes. But talk is not enough. If we were to actually construct networks instead of trees to describe the evolutionary process, what would they look like? Here we consider endosymbiosis and an example of a network of genomes involving 181 sequenced prokaryotes and how that squares off with some ideas about early cell evolution.

Multiple replication origins of Halobacterium sp. strain NRC-1: properties of the conserved orc7-dependent oriC1

The eukaryote-like DNA replication system of the model haloarchaeon Halobacterium NRC-1 is encoded within a circular chromosome and two large megaplasmids or minichromosomes, pNRC100 and pNRC200. We previously showed by genetic analysis that 2 (orc2 and orc10) of the 10 genes coding for Orc-Cdc6 replication initiator proteins were essential, while a third (orc7), located near a highly conserved autonomously replicating sequence, oriC1, was nonessential for cell viability. Here we used whole-genome marker frequency analysis (MFA) and found multiple peaks, indicative of multiple replication origins. The largest chromosomal peaks were located proximal to orc7 (oriC1) and orc10 (oriC2), and the largest peaks on the extrachromosomal elements were near orc9 (oriP1) in both pNRC100 and -200 and near orc4 (oriP2) in pNRC200. MFA of deletion strains containing different combinations of chromosomal orc genes showed that replication initiation at oriC1 requires orc7 but not orc6 and orc8. The initiation sites at oriC1 were determined by replication initiation point analysis and found to map divergently within and near an AT-rich element flanked by likely Orc binding sites. The oriC1 region, Orc binding sites, and orc7 gene orthologs were conserved in all sequenced haloarchaea. Serial deletion of orc genes resulted in the construction of a minimal strain containing not only orc2 and orc10 but also orc9. Our results suggest that replication in this model system is intriguing and more complex than previously thought. We discuss these results from the perspective of the replication strategy and evolution of haloarchaeal genomes.

A critique of prokaryotic species concepts

Developments in DNA sequencing and population genetics analyses have revealed unanticipated complexity in prokaryotic variation. The observation that genetic traits are horizontally inherited at unexpected rates within and between closely related asexually reproducing lineages impacts our comprehension of prokaryotic evolution and ecology. As a result, the concepts that point to species as being discrete clusters or monophyletic lineages are at odds with most of the data, suggesting that taxon circumscription can only proceed by informed compromise, pragmatism, and subjectivity.

Stability and solubility of proteins from extremophiles

Charges are important for hyperthermophile protein structure and function. However, the number of charges and their predicted contributions to folded state stability are not correlated, implying that more charge does not imply greater stability. The charge properties that distinguish hyperthermophile proteins also differentiate psychrophile proteins from mesophile proteins, but in the opposite direction and to a smaller extent. We conclude that charge number relates to solubility, whereas protein stability is determined by charge location. Most other structural properties are poorly separated over the ambient temperature range, apart from the burial of certain amino acids. Of particular interest are large non-polar sidechains that tend to increased exposure in proteins evolved to function at higher temperatures. Looking at tryptophan in more detail, this increase is often located close to the termini of secondary structure elements, and is discussed in terms of a novel potential role in protein thermostabilisation.

Environmental proteomics: applications of proteome profiling in environmental microbiology and biotechnology

In this review, we present the use of proteomics to advance knowledge in the field of environmental biotechnology, including studies of bacterial physiology, metabolism and ecology. Bacteria are widely applied in environmental biotechnology for their ability to catalyze dehalogenation, methanogenesis, denitrification and sulfate reduction, among others. Their tolerance to radiation and toxic compounds is also of importance. Proteomics has an important role in helping uncover the pathways behind these cellular processes. Environmental samples are often highly complex, which makes proteome studies in this field especially challenging. Some of these challenges are the lack of genome sequences for the vast majority of environmental bacteria, difficulties in isolating bacteria and proteins from certain environments, and the presence of complex microbial communities. Despite these challenges, proteomics offers a unique dynamic view into cellular function. We present examples of environmental proteomics of model organisms, and then discuss metaproteomics (microbial community proteomics), which has the potential to provide insights into the function of a community without isolating organisms. Finally, the environmental proteomics literature is summarized as it pertains to the specific application areas of wastewater treatment, metabolic engineering, microbial ecology and environmental stress responses.

The fate of new bacterial genes

Bacteria experience a continual influx of novel genetic material from a wide range of sources and yet their genomes remain relatively small. This aspect of bacterial evolution indicates that most newly arriving sequences are rapidly eliminated; however, numerous new genes persist, as evident from the presence of unique genes in almost all bacterial genomes. This review summarizes the methods for identifying new genes in bacterial genomes and examines the features that promote the retention and elimination of these evolutionary novelties.

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.

Halobacterium salinarum NRC-1 PeptideAtlas: toward strategies for targeted proteomics and improved proteome coverage

The relatively small numbers of proteins and fewer possible post-translational modifications in microbes provide a unique opportunity to comprehensively characterize their dynamic proteomes. We have constructed a PeptideAtlas (PA) covering 62.7% of the predicted proteome of the extremely halophilic archaeon Halobacterium salinarum NRC-1 by compiling approximately 636 000 tandem mass spectra from 497 mass spectrometry runs in 88 experiments. Analysis of the PA with respect to biophysical properties of constituent peptides, functional properties of parent proteins of detected peptides, and performance of different mass spectrometry approaches has highlighted plausible strategies for improving proteome coverage and selecting signature peptides for targeted proteomics. Notably, discovery of a significant correlation between absolute abundances of mRNAs and proteins has helped identify low abundance of proteins as the major limitation in peptide detection. Furthermore, we have discovered that iTRAQ labeling for quantitative proteomic analysis introduces a significant bias in peptide detection by mass spectrometry. Therefore, despite identifying at least one proteotypic peptide for almost all proteins in the PA, a context-dependent selection of proteotypic peptides appears to be the most effective approach for targeted proteomics.

Millimeter-scale genetic gradients and community-level molecular convergence in a hypersaline microbial mat

To investigate the extent of genetic stratification in structured microbial communities, we compared the metagenomes of 10 successive layers of a phylogenetically complex hypersaline mat from Guerrero Negro, Mexico. We found pronounced millimeter-scale genetic gradients that were consistent with the physicochemical profile of the mat. Despite these gradients, all layers displayed near-identical and acid-shifted isoelectric point profiles due to a molecular convergence of amino-acid usage, indicating that hypersalinity enforces an overriding selective pressure on the mat community.

Genomics and functional genomics with haloarchaea

The first haloarchaeal genome was published in 2000 and today five genome sequences are available. Transcriptome and proteome analyses have been established for two and three haloarchaeal species, respectively, and more than 20 studies using these functional genomic approaches have been published in the last two years. These studies gave global overviews of metabolic regulation (aerobic and anaerobic respiration, phototrophy, carbon source usage), stress response (UV, X-rays, transition metals, osmotic and temperature stress), cell cycle-dependent transcript level regulation, and transcript half-lives. The only translatome analysis available for any prokaryotic species revealed that 10 and 20% of all transcripts are translationally regulated in Haloferax volcanii and Halobacterium salinarum, respectively. Very effective methods for the construction of in frame deletion mutants have been established recently for haloarchaea and are intensively used to unravel the biological roles of genes in this group. Bioinformatic analyses include both cross-genome comparisons as well as integration of genomic data with experimental results. The first systems biology approaches have been performed that used experimental data to construct predictive models of gene expression and metabolism, respectively. In this contribution the current status of genomics, functional genomics, and molecular genetics of haloarchaea is summarized and selected examples are discussed.

Signature genes as a phylogenomic tool

Gene content has been shown to contain a strong phylogenetic signal, yet its usage for phylogenetic questions is hampered by horizontal gene transfer and parallel gene loss and until now required completely sequenced genomes. Here, we introduce an approach that allows the phylogenetic signal in gene content to be applied to any set of sequences, using signature genes for phylogenetic classification. The hundreds of publicly available genomes allow us to identify signature genes at various taxonomic depths, and we show how the presence of signature genes in an unspecified sample can be used to characterize its taxonomic composition. We identify 8,362 signature genes specific for 112 prokaryotic taxa. We show that these signature genes can be used to address phylogenetic questions on the basis of gene content in cases where classic gene content or sequence analyses provide an ambiguous answer, such as for Nanoarchaeum equitans, and even in cases where complete genomes are not available, such as for metagenomics data. Cross-validation experiments leaving out up to 30% of the species show that ∼92% of the signature genes correctly place the species in a related clade. Analyses of metagenomics data sets with the signature gene approach are in good agreement with the previously reported species distributions based on phylogenetic analysis of marker genes. Summarizing, signature genes can complement traditional sequence-based methods in addressing taxonomic questions.

Microbial life at high salt concentrations: phylogenetic and metabolic diversity

Halophiles are found in all three domains of life. Within the Bacteria we know halophiles within the phyla Cyanobacteria, Proteobacteria, Firmicutes, Actinobacteria, Spirochaetes, and Bacteroidetes. Within the Archaea the most salt-requiring microorganisms are found in the class Halobacteria. Halobacterium and most of its relatives require over 100–150 g/l salt for growth and structural stability. Also within the order Methanococci we encounter halophilic species. Halophiles and non-halophilic relatives are often found together in the phylogenetic tree, and many genera, families and orders have representatives with greatly different salt requirement and tolerance. A few phylogenetically coherent groups consist of halophiles only: the order Halobacteriales, family Halobacteriaceae (Euryarchaeota) and the anaerobic fermentative bacteria of the order Halanaerobiales (Firmicutes). The family Halomonadaceae (Gammaproteobacteria) almost exclusively contains halophiles. Halophilic microorganisms use two strategies to balance their cytoplasm osmotically with their medium. The first involves accumulation of molar concentrations of KCl. This strategy requires adaptation of the intracellular enzymatic machinery, as proteins should maintain their proper conformation and activity at near-saturating salt concentrations. The proteome of such organisms is highly acidic, and most proteins denature when suspended in low salt. Such microorganisms generally cannot survive in low salt media. The second strategy is to exclude salt from the cytoplasm and to synthesize and/or accumulate organic 'compatible' solutes that do not interfere with enzymatic activity. Few adaptations of the cells' proteome are needed, and organisms using the 'organic-solutes-in strategy' often adapt to a surprisingly broad salt concentration range. Most halophilic Bacteria, but also the halophilic methanogenic Archaea use such organic solutes. A variety of such solutes are known, including glycine betaine, ectoine and other amino acid derivatives, sugars and sugar alcohols. The 'high-salt-in strategy' is not limited to the Halobacteriaceae. The Halanaerobiales (Firmicutes) also accumulate salt rather than organic solutes. A third, phylogenetically unrelated organism accumulates KCl: the red extremely halophilic Salinibacter (Bacteroidetes), recently isolated from saltern crystallizer brines. Analysis of its genome showed many points of resemblance with the Halobacteriaceae, probably resulting from extensive horizontal gene transfer. The case of Salinibacter shows that more unusual types of halophiles may be waiting to be discovered.

Molecular signature of hypersaline adaptation: insights from genome and proteome composition of halophilic prokaryotes

BACKGROUND

Halophilic prokaryotes are adapted to thrive in extreme conditions of salinity. Identification and analysis of distinct macromolecular characteristics of halophiles provide insight into the factors responsible for their adaptation to high-salt environments. The current report presents an extensive and systematic comparative analysis of genome and proteome composition of halophilic and non-halophilic microorganisms, with a view to identify such macromolecular signatures of haloadaptation.

RESULTS

Comparative analysis of the genomes and proteomes of halophiles and non-halophiles reveals some common trends in halophiles that transcend the boundary of phylogenetic relationship and the genomic GC-content of the species. At the protein level, halophilic species are characterized by low hydrophobicity, over-representation of acidic residues, especially Asp, under-representation of Cys, lower propensities for helix formation and higher propensities for coil structure. At the DNA level, the dinucleotide abundance profiles of halophilic genomes bear some common characteristics, which are quite distinct from those of non-halophiles, and hence may be regarded as specific genomic signatures for salt-adaptation. The synonymous codon usage in halophiles also exhibits similar patterns regardless of their long-term evolutionary history.

CONCLUSION

The generality of molecular signatures for environmental adaptation of extreme salt-loving organisms, demonstrated in the present study, advocates the convergent evolution of halophilic species towards specific genome and amino acid composition, irrespective of their varying GC-bias and widely disparate taxonomic positions. The adapted features of halophiles seem to be related to physical principles governing DNA and protein stability, in response to the extreme environmental conditions under which they thrive.

Cloning, Purification, and Partial Characterization of the Halobacterium sp. NRC-1 Minichromosome Maintenance (MCM) Helicase

The MCM gene from the archaeon Halobacterium, with and without its intein, was cloned into an Escherichia coli expression vector, overexpressed and the protein was purified and antibodies were generated. The antibodies were used to demonstrate that in vivo only the processed enzyme, without the intein, could be detected.

A site- and time-heterogeneous model of amino acid replacement

We combined the category (CAT) mixture model (Lartillot N, Philippe H. 2004) and the nonstationary break point (BP) model (Blanquart S, Lartillot N. 2006) into a new model, CAT-BP, accounting for variations of the evolutionary process both along the sequence and across lineages. As in CAT, the model implements a mixture of distinct Markovian processes of substitution distributed among sites, thus accommodating site-specific selective constraints induced by protein structure and function. Furthermore, as in BP, these processes are nonstationary, and their equilibrium frequencies are allowed to change along lineages in a correlated way, through discrete shifts in global amino acid composition distributed along the phylogenetic tree. We implemented the CAT-BP model in a Bayesian Markov Chain Monte Carlo framework and compared its predictions with those of 3 simpler models, BP, CAT, and the site- and time-homogeneous general time-reversible (GTR) model, on a concatenation of 4 mitochondrial proteins of 20 arthropod species. In contrast to GTR, BP, and CAT, which all display a phylogenetic reconstruction artifact positioning the bees Apis mellifera and Melipona bicolor among chelicerates, the CAT-BP model is able to recover the monophyly of insects. Using posterior predictive tests, we further show that the CAT-BP combination yields better anticipations of site- and taxon-specific amino acid frequencies and that it better accounts for the homoplasies that are responsible for the artifact. Altogether, our results show that the joint modeling of heterogeneities across sites and along time results in a synergistic improvement of the phylogenetic inference, indicating that it is essential to disentangle the combined effects of both sources of heterogeneity, in order to overcome systematic errors in protein phylogenetic analyses.

Extremely Radiation-Resistant Mutants of a Halophilic Archaeon with Increased Single-Stranded DNA-Binding Protein (RPA) Gene Expression

Extremely halophilic archaea are highly resistant to multiple stressors, including radiation, desiccation and salinity. To study the basis of stress resistance and determine the maximum tolerance to ionizing radiation, we exposed cultures of the model halophile Halobacterium sp. NRC-1 to four cycles of irradiation with high doses of 18-20 MeV electrons. Two independently obtained mutants displayed an LD(50) > 11 kGy, which is higher than the LD(50) of the extremely radiation-resistant bacterium Deinococcus radiodurans. Whole-genome transcriptome analysis comparing the mutants to the parental wild-type strain revealed up-regulation of an operon containing two single-stranded DNA-binding protein (RPA) genes, VNG2160 (rfa3) and VNG2162, and a third gene of unknown function, VNG2163. The putative transcription start site for the rfa3 operon was mapped approximately 40 bp upstream of the ATG start codon, and a classical TATA-box motif was found centered about 25 bp further upstream. We propose that RPA facilitates DNA repair machinery and/or protects repair intermediates to maximize the ionizing radiation resistance of this archaeon.

Large scale comparative codon-pair context analysis unveils general rules that fine-tune evolution of mRNA primary structure

Background

Codon usage and codon-pair context are important gene primary structure features that influence mRNA decoding fidelity. In order to identify general rules that shape codon-pair context and minimize mRNA decoding error, we have carried out a large scale comparative codon-pair context analysis of 119 fully sequenced genomes.

Methodologies/Principal Findings

We have developed mathematical and software tools for large scale comparative codon-pair context analysis. These methodologies unveiled general and species specific codon-pair context rules that govern evolution of mRNAs in the 3 domains of life. We show that evolution of bacterial and archeal mRNA primary structure is mainly dependent on constraints imposed by the translational machinery, while in eukaryotes DNA methylation and tri-nucleotide repeats impose strong biases on codon-pair context.

Conclusions

The data highlight fundamental differences between prokaryotic and eukaryotic mRNA decoding rules, which are partially independent of codon usage.

Structure and function of cold shock proteins in archaea

Archaea are abundant and drive critical microbial processes in the Earth's cold biosphere. Despite this, not enough is known about the molecular mechanisms of cold adaptation and no biochemical studies have been performed on stenopsychrophilic archaea (e.g., Methanogenium frigidum). This study examined the structural and functional properties of cold shock proteins (Csps) from archaea, including biochemical analysis of the Csp from M. frigidum. csp genes are present in most bacteria and some eucarya but absent from most archaeal genome sequences, most notably, those of all archaeal thermophiles and hyperthermophiles. In bacteria, Csps are small, nucleic acid binding proteins involved in a variety of cellular processes, such as transcription. In this study, archaeal Csp function was assessed by examining the ability of csp genes from psychrophilic and mesophilic Euryarchaeota and Crenarchaeota to complement a cold-sensitive growth defect in Escherichia coli. In addition, an archaeal gene with a cold shock domain (CSD) fold but little sequence identity to Csps was also examined. Genes encoding Csps or a CSD structural analog from three psychrophilic archaea rescued the E. coli growth defect. The three proteins were predicted to have a higher content of solvent-exposed basic residues than the noncomplementing proteins, and the basic residues were located on the nucleic acid binding surface, similar to their arrangement in E. coli CspA. The M. frigidum Csp was purified and found to be a single-domain protein that folds by a reversible two-state mechanism and to exhibit a low conformational stability typical of cold-adapted proteins. Moreover, M. frigidum Csp was characterized as binding E. coli single-stranded RNA, consistent with its ability to complement function in E. coli. The studies show that some Csp and CSD fold proteins have retained sufficient similarity throughout evolution in the Archaea to be able to function effectively in the Bacteria and that the function of the archaeal proteins relates to cold adaptation. The initial biochemical analysis of M. frigidum Csp has developed a platform for further characterization and demonstrates the potential for expanding molecular studies of proteins from this important archaeal stenopsychrophile.

Transcriptional profiling of the model Archaeon Halobacterium sp. NRC-1: responses to changes in salinity and temperature

Background

The model halophile Halobacterium sp. NRC-1 was among the first Archaea to be completely sequenced and many post-genomic tools, including whole genome DNA microarrays are now being applied to its analysis. This extremophile displays tolerance to multiple stresses, including high salinity, extreme (non-mesophilic) temperatures, lack of oxygen, and ultraviolet and ionizing radiation.

Results

In order to study the response of Halobacterium sp. NRC-1 to two common stressors, salinity and temperature, we used whole genome DNA microarrays to assay for changes in gene expression under differential growth conditions. Cultures grown aerobically in rich medium at 42°C were compared to cultures grown at elevated or reduced temperature and high or low salinity. The results obtained were analyzed using a custom database and microarray analysis tools. Growth under salt stress conditions resulted in the modulation of genes coding for many ion transporters, including potassium, phosphate, and iron transporters, as well as some peptide transporters and stress proteins. Growth at cold temperature altered the expression of genes involved in lipid metabolism, buoyant gas vesicles, and cold shock proteins. Heat shock showed induction of several known chaperone genes. The results showed that Halobacterium sp. NRC-1 cells are highly responsive to environmental changes at the level of gene expression.

Conclusion

Transcriptional profiling showed that Halobacterium sp. NRC-1 is highly responsive to its environment and provided insights into some of the specific responses at the level of gene expression. Responses to changes in salt conditions appear to be designed to minimize the loss of essential ionic species and abate possible toxic effects of others, while exposure to temperature extremes elicit responses to promote protein folding and limit factors responsible for growth inhibition. This work lays the foundation for further bioinformatic and genetic studies which will lead to a more comprehensive understanding of the biology of a model halophilic Archaeon.

Global analysis of mRNA decay in Halobacterium salinarum NRC-1 at single-gene resolution using DNA microarrays

RNA degradation is an important factor in the regulation of gene expression. It allows organisms to quickly respond to changing environmental conditions by adapting the expression of individual genes. The stability of individual mRNAs within an organism varies considerably, contributing to differential amounts of proteins expressed. In this study we used DNA microarrays to analyze mRNA degradation in exponentially growing cultures of the extremely halophilic euryarchaeon Halobacterium salinarum NRC-1 on a global level. We determined mRNA half-lives for 1,717 open reading frames, 620 of which are part of known or predicted operons. Under the tested conditions transcript stabilities ranged from 5 min to more than 18 min, with 79% of the evaluated mRNAs showing half-lives between 8 and 12 min. The overall mean half-life was 10 min, which is considerably longer than the ones found in the other prokaryotes investigated thus far. As previously observed in Escherichia coli and Saccharomyces cerevisiae, we could not detect a significant correlation between transcript length and transcript stability, but there was a relationship between gene function and transcript stability. Genes that are known or predicted to be transcribed in operons exhibited similar mRNA half-lives. These results provide initial insights into mRNA turnover in a euryarchaeon. Moreover, our model organism, H. salinarum NRC-1, is one of just two archaea sequenced to date that are missing the core subunits of the archaeal exosome. This complex orthologous to the RNA degrading exosome of eukarya is found in all other archaeal genomes sequenced thus far.

Genomic survey of sequence features for ultraviolet tolerance in haloarchaea (family Halobacteriaceae)

We have investigated the strategy of Halobacterium sp. NRC-1 and other members of the family Halobacteriaceae to survive ultraviolet (UV) irradiation, based on an integrated analysis of various genomic and proteomic features such as dinucleotide composition and distribution of tetranucleotides in the genome and amino acid composition of the proteins. The low dipyrimidine content may help Halobacterium reduce formation of photoproducts in its genome. The usage of residues susceptible to reactive oxygen species attack is reduced significantly in Halobacterium, which helps the organism to minimize protein damage. We then correlated the expression of the zim gene with the genomic structure to reexamine the importance of the putative mismatch repair pathway proposed previously. Our results showed that Halobacterium sp. NRC-1 and other haloarchaea (Haloarcula marismortui, Haloquadratum walsbyi) have optimized their genomic and proteomic structures to reduce damage induced by UV irradiation, often present at high levels in habitats where these organisms thrive.

The relationships between the isoelectric point and: length of proteins, taxonomy and ecology of organisms

Background

The distribution of isoelectric point (pI) of proteins in a proteome is universal for all organisms. It is bimodal dividing the proteome into two sets of acidic and basic proteins. Different species however have different abundance of acidic and basic proteins that may be correlated with taxonomy, subcellular localization, ecological niche of organisms and proteome size.

Results

We have analysed 1784 proteomes encoded by chromosomes of Archaea, Bacteria, Eukaryota, and also mitochondria, plastids, prokaryotic plasmids, phages and viruses. We have found significant correlation in more than 95% of proteomes between the protein length and pI in proteomes – positive for acidic proteins and negative for the basic ones. Plastids, viruses and plasmids encode more basic proteomes while chromosomes of Archaea, Bacteria, Eukaryota, mitochondria and phages more acidic ones. Mitochondrial proteomes of Viridiplantae, Protista and Fungi are more basic than Metazoa. It results from the presence of basic proteins in the former proteomes and their absence from the latter ones and is related with reduction of metazoan genomes. Significant correlation was found between the pI bias of proteomes encoded by prokaryotic chromosomes and proteomes encoded by plasmids but there is no correlation between eukaryotic nuclear-coded proteomes and proteomes encoded by organelles. Detailed analyses of prokaryotic proteomes showed significant relationships between pI distribution and habitat, relation to the host cell and salinity of the environment, but no significant correlation with oxygen and temperature requirements. The salinity is positively correlated with acidicity of proteomes. Host-associated organisms and especially intracellular species have more basic proteomes than free-living ones. The higher rate of mutations accumulation in the intracellular parasites and endosymbionts is responsible for the basicity of their tiny proteomes that explains the observed positive correlation between the decrease of genome size and the increase of basicity of proteomes. The results indicate that even conserved proteins subjected to strong selectional constraints follow the global trend in the pI distribution.

Conclusion

The distribution of pI of proteins in proteomes shows clear relationships with length of proteins, subcellular localization, taxonomy and ecology of organisms. The distribution is also strongly affected by mutational pressure especially in intracellular organisms.

Essential and non-essential DNA replication genes in the model halophilic Archaeon, Halobacterium sp. NRC-1

Background

Information transfer systems in Archaea, including many components of the DNA replication machinery, are similar to those found in eukaryotes. Functional assignments of archaeal DNA replication genes have been primarily based upon sequence homology and biochemical studies of replisome components, but few genetic studies have been conducted thus far. We have developed a tractable genetic system for knockout analysis of genes in the model halophilic archaeon, Halobacterium sp. NRC-1, and used it to determine which DNA replication genes are essential.

Results

Using a directed in-frame gene knockout method in Halobacterium sp. NRC-1, we examined nineteen genes predicted to be involved in DNA replication. Preliminary bioinformatic analysis of the large haloarchaeal Orc/Cdc6 family, related to eukaryotic Orc1 and Cdc6, showed five distinct clades of Orc/Cdc6 proteins conserved in all sequenced haloarchaea. Of ten orc/cdc6 genes in Halobacterium sp. NRC-1, only two were found to be essential, orc10, on the large chromosome, and orc2, on the minichromosome, pNRC200. Of the three replicative-type DNA polymerase genes, two were essential: the chromosomally encoded B family, polB1, and the chromosomally encoded euryarchaeal-specific D family, polD1/D2 (formerly called polA1/polA2 in the Halobacterium sp. NRC-1 genome sequence). The pNRC200-encoded B family polymerase, polB2, was non-essential. Accessory genes for DNA replication initiation and elongation factors, including the putative replicative helicase, mcm, the eukaryotic-type DNA primase, pri1/pri2, the DNA polymerase sliding clamp, pcn, and the flap endonuclease, rad2, were all essential. Targeted genes were classified as non-essential if knockouts were obtained and essential based on statistical analysis and/or by demonstrating the inability to isolate chromosomal knockouts except in the presence of a complementing plasmid copy of the gene.

Conclusion

The results showed that ten out of nineteen eukaryotic-type DNA replication genes are essential for Halobacterium sp. NRC-1, consistent with their requirement for DNA replication. The essential genes code for two of ten Orc/Cdc6 proteins, two out of three DNA polymerases, the MCM helicase, two DNA primase subunits, the DNA polymerase sliding clamp, and the flap endonuclease.

Sequence comparison and environmental adaptation of a bacterial endonuclease

The periplasmic/extracellular bacterial enzyme endonuclease I was chosen as a model system to identify features that might be responsible for temperature- and salt adaptation. A statistical study of amino acid sequence properties belonging to endonuclease I enzymes from three mesophilic habitats (non-marine, brackish water and marine), and three marine temperature groups (psychrophile, intermediate and mesophile) has been conducted. Ten new endonuclease I genes have been sequenced in order to increase the sample size. A bioinformatical method of property dependent statistical analysis of alignments has been applied. To our knowledge this is the first time these methods have been used in order to investigate environmental adaptation of enzymes. Adaptation to low temperature seems to involve increased surface isoelectric point and hydrophobicity in contrast to salt adaptation in which the isoelectric point and hydrophobicity at the surface decreases. Redistribution of charge and hydrophobicity might be the most important signature for cold adaptation and salt adaptation of this enzyme class. The results indicate that general trends of adaptation are possible to elucidate from the amino acid sequences. Also in this paper a new scale of stratified B-factors, derived from the Protein Data Bank, is presented.

Genomewide and biochemical analyses of DNA-binding activity of Cdc6/Orc1 and Mcm proteins in Pyrococcus sp

The origin of DNA replication (oriC) of the hyperthermophilic archaeon Pyrococcus abyssi contains multiple ORB and mini-ORB repeats that show sequence similarities to other archaeal ORB (origin recognition box). We report here that the binding of Cdc6/Orc1 to a 5 kb region containing oriC in vivo was highly specific both in exponential and stationary phases, by means of chromatin immunoprecipitation coupled with hybridization on a whole genome microarray (ChIP-chip). The oriC region is practically the sole binding site for the Cdc6/Orc1, thereby distinguishing oriC in the 1.8 M bp genome. We found that the 5 kb region contains a previously unnoticed cluster of ORB and mini-ORB repeats in the gene encoding the small subunit (dp1) for DNA polymerase II (PolD). ChIP and the gel retardation analyses further revealed that Cdc6/Orc1 specifically binds both of the ORB clusters in oriC and dp1. The organization of the ORB clusters in the dp1 and oriC is conserved during evolution in the order Thermococcales, suggesting a role in the initiation of DNA replication. Our ChIP-chip analysis also revealed that Mcm alters the binding specificity to the oriC region according to the growth phase, consistent with its role as a licensing factor.

Large-scale identification of N-terminal peptides in the halophilic archaea Halobacterium salinarum and Natronomonas pharaonis

Characterization of protein N-terminal peptides supports the quality assessment of data derived from genomic sequences (e.g., the correct assignment of start codons) and hints to in vivo N-terminal modifications such as N-terminal acetylation and removal of the initiator methionine. The current work represents the first large-scale identification of N-terminal peptides from prokaryotes, of the two halophilic euryarchaeota Halobacterium salinarum and Natronomonas pharaonis. Two methods were used that specifically allow the characterization of protein N-terminal peptides: combined fractional diagonal chromatography (COFRADIC) and strong cation exchange chromatography (SCX), both known to enrich for N-terminally blocked peptides. In addition to these specific methods, N-terminal peptide identifications were extracted from our previous genome-wide proteomic data. Combining all data, 606 N-terminal peptides from Hbt. salinarum and 328 from Nmn. pharaonis were reliably identified. These results constitute the largest available dataset holding identified and characterized protein N-termini for prokaryotes (archaea and bacteria). They allowed the validation/improvement of start codon assignments as automatic gene finders tend to misassign start codons for GC-rich genomes. In addition, the dataset allowed unravelling N-terminal protein maturation in archaea, showing that 60% of the proteins undergo methionine cleavage and that-in contrast to current knowledge-Nalpha-acetylation is common in the archaeal domain of life with 13-18% of the proteins being Nalpha-acetylated. The protein sets described in this paper are available by FTP and might be used as reference sets to test the performance of new gene finders.

Genetic and physical mapping of DNA replication origins in Haloferax volcanii

The halophilic archaeon Haloferax volcanii has a multireplicon genome, consisting of a main chromosome, three secondary chromosomes, and a plasmid. Genes for the initiator protein Cdc6/Orc1, which are commonly located adjacent to archaeal origins of DNA replication, are found on all replicons except plasmid pHV2. However, prediction of DNA replication origins in H. volcanii is complicated by the fact that this species has no less than 14 cdc6/orc1 genes. We have used a combination of genetic, biochemical, and bioinformatic approaches to map DNA replication origins in H. volcanii. Five autonomously replicating sequences were found adjacent to cdc6/orc1 genes and replication initiation point mapping was used to confirm that these sequences function as bidirectional DNA replication origins in vivo. Pulsed field gel analyses revealed that cdc6/orc1-associated replication origins are distributed not only on the main chromosome (2.9 Mb) but also on pHV1 (86 kb), pHV3 (442 kb), and pHV4 (690 kb) replicons. Gene inactivation studies indicate that linkage of the initiator gene to the origin is not required for replication initiation, and genetic tests with autonomously replicating plasmids suggest that the origin located on pHV1 and pHV4 may be dominant to the principal chromosomal origin. The replication origins we have identified appear to show a functional hierarchy or differential usage, which might reflect the different replication requirements of their respective chromosomes. We propose that duplication of H. volcanii replication origins was a prerequisite for the multireplicon structure of this genome, and that this might provide a means for chromosome-specific replication control under certain growth conditions. Our observations also suggest that H. volcanii is an ideal organism for studying how replication of four replicons is regulated in the context of the archaeal cell cycle.

Phylogenomic analysis of proteins that are distinctive of Archaea and its main subgroups and the origin of methanogenesis

BACKGROUND

The Archaea are highly diverse in terms of their physiology, metabolism and ecology. Presently, very few molecular characteristics are known that are uniquely shared by either all archaea or the different main groups within archaea. The evolutionary relationships among different groups within the Euryarchaeota branch are also not clearly understood.

RESULTS

We have carried out comprehensive analyses on each open reading frame (ORFs) in the genomes of 11 archaea (3 Crenarchaeota--Aeropyrum pernix, Pyrobaculum aerophilum and Sulfolobus acidocaldarius; 8 Euryarchaeota--Pyrococcus abyssi, Methanococcus maripaludis, Methanopyrus kandleri, Methanococcoides burtonii, Halobacterium sp. NCR-1, Haloquadratum walsbyi, Thermoplasma acidophilum and Picrophilus torridus) to search for proteins that are unique to either all Archaea or for its main subgroups. These studies have identified 1448 proteins or ORFs that are distinctive characteristics of Archaea and its various subgroups and whose homologues are not found in other organisms. Six of these proteins are unique to all Archaea, 10 others are only missing in Nanoarchaeum equitans and a large number of other proteins are specific for various main groups within the Archaea (e.g. Crenarchaeota, Euryarchaeota, Sulfolobales and Desulfurococcales, Halobacteriales, Thermococci, Thermoplasmata, all methanogenic archaea or particular groups of methanogens). Of particular importance is the observation that 31 proteins are uniquely present in virtually all methanogens (including M. kandleri) and 10 additional proteins are only found in different methanogens as well as A. fulgidus. In contrast, no protein was exclusively shared by various methanogen and any of the Halobacteriales or Thermoplasmatales. These results strongly indicate that all methanogenic archaea form a monophyletic group exclusive of other archaea and that this lineage likely evolved from Archaeoglobus. In addition, 15 proteins that are uniquely shared by M. kandleri and Methanobacteriales suggest a close evolutionary relationship between them. In contrast to the phylogenomics studies, a monophyletic grouping of archaea is not supported by phylogenetic analyses based on protein sequences.

CONCLUSION

The identified archaea-specific proteins provide novel molecular markers or signature proteins that are distinctive characteristics of Archaea and all of its major subgroups. The species distributions of these proteins provide novel insights into the evolutionary relationships among different groups within Archaea, particularly regarding the origin of methanogenesis. Most of these proteins are of unknown function and further studies should lead to discovery of novel biochemical and physiological characteristics that are unique to either all archaea or its different subgroups.

Amino acid contribution to protein solubility: Asp, Glu, and Ser contribute more favorably than the other hydrophilic amino acids in RNase Sa

Poor protein solubility is a common problem in high-resolution structural studies, formulation of protein pharmaceuticals, and biochemical characterization of proteins. One popular strategy to improve protein solubility is to use site-directed mutagenesis to make hydrophobic to hydrophilic mutations on the protein surface. However, a systematic investigation of the relative contributions of all 20 amino acids to protein solubility has not been done. Here, 20 variants at the completely solvent-exposed position 76 of ribonuclease (RNase) Sa are made to compare the contributions of each amino acid. Stability measurements were also made for these variants, which occur at the i+1 position of a type II beta-turn. Solubility measurements in ammonium sulfate solutions were made at high positive net charge, low net charge, and high negative net charge. Surprisingly, there was a wide range of contributions to protein solubility even among the hydrophilic amino acids. The results suggest that aspartic acid, glutamic acid, and serine contribute significantly more favorably than the other hydrophilic amino acids especially at high net charge. Therefore, to increase protein solubility, asparagine, glutamine, or threonine should be replaced with aspartic acid, glutamic acid or serine.

The role of laterally transferred genes in adaptive evolution

Background

Bacterial genomes develop new mechanisms to tide them over the imposing conditions they encounter during the course of their evolution. Acquisition of new genes by lateral gene transfer may be one of the dominant ways of adaptation in bacterial genome evolution. Lateral gene transfer provides the bacterial genome with a new set of genes that help it to explore and adapt to new ecological niches.

Methods

A maximum likelihood analysis was done on the five sequenced corynebacterial genomes to model the rates of gene insertions/deletions at various depths of the phylogeny.

Results

The study shows that most of the laterally acquired genes are transient and the inferred rates of gene movement are higher on the external branches of the phylogeny and decrease as the phylogenetic depth increases. The newly acquired genes are under relaxed selection and evolve faster than their older counterparts. Analysis of some of the functionally characterised LGTs in each species has indicated that they may have a possible adaptive role.

Conclusion

The five Corynebacterial genomes sequenced to date have evolved by acquiring between 8 – 14% of their genomes by LGT and some of these genes may have a role in adaptation.

Cold-adapted enzymes

By far the largest proportion of the Earth's biosphere is comprised of organisms that thrive in cold environments (psychrophiles). Their ability to proliferate in the cold is predicated on a capacity to synthesize cold-adapted enzymes. These enzymes have evolved a range of structural features that confer a high level of flexibility compared to thermostable homologs. High flexibility, particularly around the active site, is translated into low-activation enthalpy, low-substrate affinity, and high specific activity at low temperatures. High flexibility is also accompanied by a trade-off in stability, resulting in heat lability and, in the few cases studied, cold lability. This review addresses the structure, function, and stability of cold-adapted enzymes, highlighting the challenges for immediate and future consideration. Because of the unique properties of cold-adapted enzymes, they are not only an important focus in extremophile biology, but also represent a valuable model for fundamental research into protein folding and catalysis.

Proteome analysis of Halobacterium salinarum and characterization of proteins related to the degradation of isopropyl alcohol

We reported in a previous study that proteomic approach, coupled with genomic techniques, could be used to screen and develop multiple candidates for halophilic enzymes from Halobacterium salinarum. In order to evaluate the biodegradation of isopropyl alcohol (IPA) by H. salinarum, the amounts of residual IPA and acetone generated in the growth media were determined using a gas chromatography-flame ionization detector (GC-FID). The protein expression profiles of cells which had been cultured with IPA were obtained with the two-dimensional gel electrophoresis. Proteins evidencing different expression levels in the presence of 0.5% IPA were identified by electrospray ionization-quadruple-time of flight (ESI-Q-TOF) mass spectrometry. We found 12 proteins which were down-regulated, and another 12 proteins which were up-regulated, in the presence of 0.5% IPA and we further identified 17 proteins among them using ESI-TOF MS/MS. Among these identified proteins, we selected glyceraldehyde 3-phosphate dehydrogenase (GAPDH) for further characterization as a halophilic enzyme. We have demonstrated for the first time that H. salinarum possesses the ability to degrade IPA and GAPDH was both stable and active at high salt concentrations, with maximum activity occurring at 1 M NaCl, although the optimal salt concentration with regard to the growth of H. salinarum is 4.3 M.

Molecular characterization of the minimal replicon and the unidirectional theta replication of pSCM201 in extremely halophilic archaea

A 3,463-bp plasmid, pSCM201, was isolated from a halophilic archaeon, Haloarcula sp. strain AS7094. The minimal replicon that is essential and sufficient for autonomous replication and stable maintenance in Haloarcula hispanica was determined by deletion analysis of the plasmid. This minimal replicon ( approximately 1.8 kb) consisted of only two functionally related segments: (i) a putative origin (ori201) containing an AT-rich region and sets of repeats and (ii) an adjacent gene encoding a putative replication initiation protein (Rep201). Electron microscopic observation and Southern blotting analysis demonstrated that pSCM201 replicates via a theta mechanism. Precise mapping of the putative origin suggested that the replication initiated from a fixed site close to the AT-rich region and proceeded unidirectionally toward the downstream rep201 gene, which was further confirmed by electron microscopic analysis of the ClaI-digested replication intermediates. To our knowledge, this is the first unidirectional theta replication plasmid experimentally identified in the domain of archaea. It provides a novel plasmid system to conduct research on archaeal DNA replication.

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.

Extremely halophilic archaea and the issue of long-term microbial survival

Halophilic archaebacteria (haloarchaea) thrive in environments with salt concentrations approaching saturation, such as natural brines, the Dead Sea, alkaline salt lakes and marine solar salterns; they have also been isolated from rock salt of great geological age (195-250 million years). An overview of their taxonomy, including novel isolates from rock salt, is presented here; in addition, some of their unique characteristics and physiological adaptations to environments of low water activity are reviewed. The issue of extreme long-term microbial survival is considered and its implications for the search for extraterrestrial life. The development of detection methods for subterranean haloarchaea, which might also be applicable to samples from future missions to space, is presented.

PROMPT: a protein mapping and comparison tool

Background

Comparison of large protein datasets has become a standard task in bioinformatics. Typically researchers wish to know whether one group of proteins is significantly enriched in certain annotation attributes or sequence properties compared to another group, and whether this enrichment is statistically significant. In order to conduct such comparisons it is often required to integrate molecular sequence data and experimental information from disparate incompatible sources. While many specialized programs exist for comparisons of this kind in individual problem domains, such as expression data analysis, no generic software solution capable of addressing a wide spectrum of routine tasks in comparative proteomics is currently available.

Results

PROMPT is a comprehensive bioinformatics software environment which enables the user to compare arbitrary protein sequence sets, revealing statistically significant differences in their annotation features. It allows automatic retrieval and integration of data from a multitude of molecular biological databases as well as from a custom XML format. Similarity-based mapping of sequence IDs makes it possible to link experimental information obtained from different sources despite discrepancies in gene identifiers and minor sequence variation. PROMPT provides a full set of statistical procedures to address the following four use cases: i) comparison of the frequencies of categorical annotations between two sets, ii) enrichment of nominal features in one set with respect to another one, iii) comparison of numeric distributions, and iv) correlation of numeric variables. Analysis results can be visualized in the form of plots and spreadsheets and exported in various formats, including Microsoft Excel.

Conclusion

PROMPT is a versatile, platform-independent, easily expandable, stand-alone application designed to be a practical workhorse in analysing and mining protein sequences and associated annotation. The availability of the Java Application Programming Interface and scripting capabilities on one hand, and the intuitive Graphical User Interface with context-sensitive help system on the other, make it equally accessible to professional bioinformaticians and biologically-oriented users. PROMPT is freely available for academic users from .

Ten years of bacterial genome sequencing: comparative-genomics-based discoveries

It has been more than 10 years since the first bacterial genome sequence was published. Hundreds of bacterial genome sequences are now available for comparative genomics, and searching a given protein against more than a thousand genomes will soon be possible. The subject of this review will address a relatively straightforward question: "What have we learned from this vast amount of new genomic data?" Perhaps one of the most important lessons has been that genetic diversity, at the level of large-scale variation amongst even genomes of the same species, is far greater than was thought. The classical textbook view of evolution relying on the relatively slow accumulation of mutational events at the level of individual bases scattered throughout the genome has changed. One of the most obvious conclusions from examining the sequences from several hundred bacterial genomes is the enormous amount of diversity--even in different genomes from the same bacterial species. This diversity is generated by a variety of mechanisms, including mobile genetic elements and bacteriophages. An examination of the 20 Escherichia coli genomes sequenced so far dramatically illustrates this, with the genome size ranging from 4.6 to 5.5 Mbp; much of the variation appears to be of phage origin. This review also addresses mobile genetic elements, including pathogenicity islands and the structure of transposable elements. There are at least 20 different methods available to compare bacterial genomes. Metagenomics offers the chance to study genomic sequences found in ecosystems, including genomes of species that are difficult to culture. It has become clear that a genome sequence represents more than just a collection of gene sequences for an organism and that information concerning the environment and growth conditions for the organism are important for interpretation of the genomic data. The newly proposed Minimal Information about a Genome Sequence standard has been developed to obtain this information.

Protein secretion in the Archaea: multiple paths towards a unique cell surface

Archaea are similar to other prokaryotes in most aspects of cell structure but are unique with respect to the lipid composition of the cytoplasmic membrane and the structure of the cell surface. Membranes of archaea are composed of glycerol-ether lipids instead of glycerol-ester lipids and are based on isoprenoid side chains, whereas the cell walls are formed by surface-layer proteins. The unique cell surface of archaea requires distinct solutions to the problem of how proteins cross this barrier to be either secreted into the medium or assembled as appendages at the cell surface.

A systems view of haloarchaeal strategies to withstand stress from transition metals

Given that transition metals are essential cofactors in central biological processes, misallocation of the wrong metal ion to a metalloprotein can have resounding and often detrimental effects on diverse aspects of cellular physiology. Therefore, in an attempt to characterize unique and shared responses to chemically similar metals, we have reconstructed physiological behaviors of Halobacterium NRC-1, an archaeal halophile, in sublethal levels of Mn(II), Fe(II), Co(II), Ni(II), Cu(II), and Zn(II). Over 20% of all genes responded transiently within minutes of exposure to Fe(II), perhaps reflecting immediate large-scale physiological adjustments to maintain homeostasis. At steady state, each transition metal induced growth arrest, attempts to minimize oxidative stress, toxic ion scavenging, increased protein turnover and DNA repair, and modulation of active ion transport. While several of these constitute generalized stress responses, up-regulation of active efflux of Co(II), Ni(II), Cu(II), and Zn(II), down-regulation of Mn(II) uptake and up-regulation of Fe(II) chelation, confer resistance to the respective metals. We have synthesized all of these discoveries into a unified systems-level model to provide an integrated perspective of responses to six transition metals with emphasis on experimentally verified regulatory mechanisms. Finally, through comparisons across global transcriptional responses to different metals, we provide insights into putative in vivo metal selectivity of metalloregulatory proteins and demonstrate that a systems approach can help rapidly unravel novel metabolic potential and regulatory programs of poorly studied organisms.

Subcellular pH and predicted pH-dependent features of proteins

A characteristic of two-dimensional proteomics gels is a general bimodal distribution of isoelectric (pI) values. Discussion of this feature has focussed on the balance of acidic and basic ionisable residues, and potential relationships between pI distributions and organism classification or protein subcellular location. Electrostatics calculations on a set of protein structures with known subcellular location show that predicted folded state pI are similar to those calculated from sequence alone, but adjusted according to a general stabilising effect from interactions between ionisable groups. Bimodal distributions dominate both pI and the predicted pH of maximal stability. However, there are significant differences between these features. The average pH of maximal stability generally follows organelle pH. Average pI values are well removed from organelle pH in most subcellular environments, consistent with the view that proteins have evolved to carry (on average) net charge in a given subcellular location, and relevant to discussion of solubility in crowded environments. Correlation of the predicted pH of maximum stability with subcellular pH suggests an evolutionary pressure to adjust folded state interactions according to environment. Finally, our analysis of ionisable group contributions to stability suggests that Golgi proteins have the largest such term, although this dataset is small.

Proteome Analysis of Halobacterium sp. NRC-1 Facilitated by the Biomodule Analysis Tool BMSorter

To better understand the extremely halophilic archaeon Halobacterium species NRC-1, we analyzed its soluble proteome by two-dimensional liquid chromatography coupled to electrospray ionization tandem mass spectrometry. A total of 888 unique proteins were identified with a ProteinProphet probability (P) between 0.9 and 1.0. To evaluate the biochemical activities of the organism, the proteomic data were subjected to a biological network analysis using our BMSorter software. This allowed us to examine the proteins expressed in different biomodules and study the interactions between pertinent biomodules. Interestingly an integrated analysis of the enzymes in the amino acid metabolism and citrate cycle networks suggested that up to eight amino acids may be converted to oxaloacetate, fumarate, or oxoglutarate in the citrate cycle for energy production. In addition, glutamate and aspartate may be interconverted from other amino acids or synthesized from citrate cycle intermediates to meet the high demand for the acidic amino acids that are required to build the highly acidic proteome of the organism. Thus this study demonstrated that proteome analysis can provide useful information and help systems analyses of organisms.

Evolution of structure and function in the o-succinylbenzoate synthase/N-acylamino acid racemase family of the enolase superfamily

Understanding how proteins evolve to provide both exquisite specificity and proficient activity is a fundamental problem in biology that has implications for protein function prediction and protein engineering. To study this problem, we analyzed the evolution of structure and function in the o-succinylbenzoate synthase/N-acylamino acid racemase (OSBS/NAAAR) family, part of the mechanistically diverse enolase superfamily. Although all characterized members of the family catalyze the OSBS reaction, this family is extraordinarily divergent, with some members sharing <15% identity. In addition, a member of this family, Amycolatopsis OSBS/NAAAR, is promiscuous, catalyzing both dehydration and racemization. Although the OSBS/NAAAR family appears to have a single evolutionary origin, no sequence or structural motifs unique to this family could be identified; all residues conserved in the family are also found in enolase superfamily members that have different functions. Based on their species distribution, several uncharacterized proteins similar to Amycolatopsis OSBS/NAAAR appear to have been transmitted by lateral gene transfer. Like Amycolatopsis OSBS/NAAAR, these might have additional or alternative functions to OSBS because many are from organisms lacking the pathway in which OSBS is an intermediate. In addition to functional differences, the OSBS/NAAAR family exhibits surprising structural variations, including large differences in orientation between the two domains. These results offer several insights into protein evolution. First, orthologous proteins can exhibit significant structural variation, and specificity can be maintained with little conservation of ligand-contacting residues. Second, the discovery of a set of proteins similar to Amycolatopsis OSBS/NAAAR supports the hypothesis that new protein functions evolve through promiscuous intermediates. Finally, a combination of evolutionary, structural, and sequence analyses identified characteristics that might prime proteins, such as Amycolatopsis OSBS/NAAAR, for the evolution of new activities.

Proteomic analysis of acidic chaperones, and stress proteins in extreme halophile Halobacterium NRC-1: a comparative proteomic approach to study heat shock response

BACKGROUND

Halobacterium sp. NRC-1 is an extremely halophilic archaeon and has adapted to optimal growth under conditions of extremely high salinity. Its proteome is highly acidic with a median pI of 4.9, a unique characteristic which helps the organism to adapt high saline environment. In the natural growth environment, Halobacterium NRC-1 encounters a number of stressful conditions including high temperature and intense solar radiation, oxidative and cold stress. Heat shock proteins and chaperones play indispensable roles in an organism's survival under many stress conditions. The aim of this study was to develop an improved method of 2-D gel electrophoresis with enhanced resolution of the acidic proteome, and to identify proteins with diverse cellular functions using in-gel digestion and LC-MS/MS and MALDI-TOF approach.

RESULTS

A modified 2-D gel electrophoretic procedure, employing IPG strips in the range of pH 3-6, enabled improved separation of acidic proteins relative to previous techniques. Combining experimental data from 2-D gel electrophoresis with available genomic information, allowed the identification of at least 30 cellular proteins involved in many cellular functions: stress response and protein folding (CctB, PpiA, DpsA, and MsrA), DNA replication and repair (DNA polymerase A alpha subunit, Orc4/CDC6, and UvrC), transcriptional regulation (Trh5 and ElfA), translation (ribosomal proteins Rps27ae and Rphs6 of the 30 S ribosomal subunit; Rpl31eand Rpl18e of the 50 S ribosomal subunit), transport (YufN), chemotaxis (CheC2), and housekeeping (ThiC, ThiD, FumC, ImD2, GapB, TpiA, and PurE). In addition, four gene products with undetermined function were also identified: Vng1807H, Vng0683C, Vng1300H, and Vng6254. To study the heat shock response of Halobacterium NRC-1, growth conditions for heat shock were determined and the proteomic profiles under normal (42 degrees C), and heat shock (49 degrees C) conditions, were compared. Using a differential proteomic approach in combination with available genomic information, bioinformatic analysis revealed five putative heat shock proteins that were upregulated in cells subjected to heat stress at 49 degrees C, namely DnaJ, GrpE, sHsp-1, Hsp-5 and sHsp-2.

CONCLUSION

The modified 2-D gel electrophoresis markedly enhanced the resolution of the extremely acidic proteome of Halobacterium NRC-1. Constitutive expression of stress proteins and chaperones help the organism to adapt and survive under extreme salinity and other stress conditions. The upregulated expression pattern of putative chaperones DnaJ, GrpE, sHsp-1, Hsp-5 and sHsp-2 under elevated temperature clearly suggests that Halobacterium NRC-1 has a sophisticated defense mechanism to survive in extreme environments.

A plant locus essential for phylloquinone (vitamin K1) biosynthesis originated from a fusion of four eubacterial genes

Phylloquinone is a compound present in all photosynthetic plants serving as cofactor for Photosystem I-mediated electron transport. Newly identified seedling-lethal Arabidopsis thaliana mutants impaired in the biosynthesis of phylloquinone possess reduced Photosystem I activity. The affected gene, called PHYLLO, consists of a fusion of four previously individual eubacterial genes, menF, menD, menC, and menH, required for the biosynthesis of phylloquinone in photosynthetic cyanobacteria and the respiratory menaquinone in eubacteria. The fact that homologous men genes reside as polycistronic units in eubacterial chromosomes and in plastomes of red algae strongly suggests that PHYLLO derived from a plastid operon during endosymbiosis. The principle architecture of the fused PHYLLO locus is conserved in the nuclear genomes of plants, green algae, and the diatom alga Thalassiosira pseudonana. The latter arose from secondary endosymbiosis of a red algae and a eukaryotic host indicating selective driving forces for maintenance and/or independent generation of the composite gene cluster within the nuclear genomes. Besides, individual menF genes, encoding active isochorismate synthases (ICS), have been established followed by splitting of the essential 3' region of the menF module of PHYLLO only in genomes of higher plants. This resulted in inactivation of the ICS activity encoded by PHYLLO and enabled a metabolic branch from the phylloquinone biosynthetic route to independently regulate the synthesis of salicylic acid required for plant defense. Therefore, gene fusion, duplication, and fission events adapted a eubacterial multienzymatic system to the metabolic requirements of plants.

Analysis of protein solvent interactions in glucose dehydrogenase from the extreme halophile Haloferax mediterranei

The structure of glucose dehydrogenase from the extreme halophile Haloferax mediterranei has been solved at 1.6-A resolution under crystallization conditions which closely mimic the "in vivo" intracellular environment. The decoration of the enzyme's surface with acidic residues is only partially neutralized by bound potassium counterions, which also appear to play a role in substrate binding. The surface shows the expected reduction in hydrophobic character, surprisingly not from changes associated with the loss of exposed hydrophobic residues but rather arising from a loss of lysines consistent with the genome wide-reduction of this residue in extreme halophiles. The structure reveals a highly ordered, multilayered solvation shell that can be seen to be organized into one dominant network covering much of the exposed surface accessible area to an extent not seen in almost any other protein structure solved. This finding is consistent with the requirement of the enzyme to form a protective shell in a dehydrating environment.