Localizing genes on the map of the genome of Haloferax volcanii, one of the Archaea

Generation of comprehensive transposon insertion mutant library for the model archaeon, Haloferax volcanii, and its use for gene discovery

Background

Archaea share fundamental properties with bacteria and eukaryotes. Yet, they also possess unique attributes, which largely remain poorly characterized. Haloferax volcanii is an aerobic, moderately halophilic archaeon that can be grown in defined media. It serves as an excellent archaeal model organism to study the molecular mechanisms of biological processes and cellular responses to changes in the environment. Studies on haloarchaea have been impeded by the lack of efficient genetic screens that would facilitate the identification of protein functions and respective metabolic pathways.

Results

Here, we devised an insertion mutagenesis strategy that combined Mu in vitro DNA transposition and homologous-recombination-based gene targeting in H. volcanii. We generated an insertion mutant library, in which the clones contained a single genomic insertion. From the library, we isolated pigmentation-defective and auxotrophic mutants, and the respective insertions pinpointed a number of genes previously known to be involved in carotenoid and amino acid biosynthesis pathways, thus validating the performance of the methodologies used. We also identified mutants that had a transposon insertion in a gene encoding a protein of unknown or putative function, demonstrating that novel roles for non-annotated genes could be assigned.

Conclusions

We have generated, for the first time, a random genomic insertion mutant library for a halophilic archaeon and used it for efficient gene discovery. The library will facilitate the identification of non-essential genes behind any specific biochemical pathway. It represents a significant step towards achieving a more complete understanding of the unique characteristics of halophilic archaea.

Electronic supplementary material

The online version of this article (doi:10.1186/s12915-014-0103-3) contains supplementary material, which is available to authorized users.

Model organisms for genetics in the domain Archaea: methanogens, halophiles, Thermococcales and Sulfolobales

The tree of life is split into three main branches: eukaryotes, bacteria, and archaea. Our knowledge of eukaryotic and bacteria cell biology has been built on a foundation of studies in model organisms, using the complementary approaches of genetics and biochemistry. Archaea have led to some exciting discoveries in the field of biochemistry, but archaeal genetics has been slow to get off the ground, not least because these organisms inhabit some of the more inhospitable places on earth and are therefore believed to be difficult to culture. In fact, many species can be cultivated with relative ease and there has been tremendous progress in the development of genetic tools for both major archaeal phyla, the Euryarchaeota and the Crenarchaeota. There are several model organisms available for methanogens, halophiles, and thermophiles; in the latter group, there are genetic systems for Sulfolobales and Thermococcales. In this review, we present the advantages and disadvantages of working with each archaeal group, give an overview of their different genetic systems, and direct the neophyte archaeologist to the most appropriate model organism.

Homing endonucleases residing within inteins: evolutionary puzzles awaiting genetic solutions

Inteins are selfish genetic elements that disrupt the sequence of protein-coding genes and are excised post-translationally. Most inteins also contain a HEN (homing endonuclease) domain, which is important for their horizontal transmission. The present review focuses on the evolution of inteins and their nested HENs, and highlights several unsolved questions that could benefit from molecular genetic approaches. Such approaches can be well carried out in halophilic archaea, which are naturally intein-rich and have highly developed genetic tools for their study. In particular, the fitness effects of habouring an intein/HEN can be tested in direct competition assays, providing additional insights that will improve current evolutionary models.

The complete genome sequence of Haloferax volcanii DS2, a model archaeon

BACKGROUND

Haloferax volcanii is an easily culturable moderate halophile that grows on simple defined media, is readily transformable, and has a relatively stable genome. This, in combination with its biochemical and genetic tractability, has made Hfx. volcanii a key model organism, not only for the study of halophilicity, but also for archaeal biology in general.

METHODOLOGY/PRINCIPAL FINDINGS

We report here the sequencing and analysis of the genome of Hfx. volcanii DS2, the type strain of this species. The genome contains a main 2.848 Mb chromosome, three smaller chromosomes pHV1, 3, 4 (85, 438, 636 kb, respectively) and the pHV2 plasmid (6.4 kb).

CONCLUSIONS/SIGNIFICANCE

The completed genome sequence, presented here, provides an invaluable tool for further in vivo and in vitro studies of Hfx. volcanii.

In vivo analysis of an essential archaeal signal recognition particle in its native host

The evolutionarily conserved signal recognition particle (SRP) plays an integral role in Sec-mediated cotranslational protein translocation and membrane protein insertion, as it has been shown to target nascent secretory and membrane proteins to the bacterial and eukaryotic translocation pores. However, little is known about its function in archaea, since characterization of the SRP in this domain of life has thus far been limited to in vitro reconstitution studies of heterologously expressed archaeal SRP components identified by sequence comparisons. In the present study, the genes encoding the SRP54, SRP19, and 7S RNA homologs (hv54h, hv19h, and hv7Sh, respectively) of the genetically and biochemically tractable archaeon Haloferax volcanii were cloned, providing the tools to analyze the SRP in its native host. As part of this analysis, an hv54h knockout strain was created. In vivo characterization of this strain revealed that the archaeal SRP is required for viability, suggesting that cotranslational protein translocation is an essential process in archaea. Furthermore, a method for the purification of this SRP employing nickel chromatography was developed in H. volcanii, allowing the successful copurification of (i) Hv7Sh with a histidine-tagged Hv54h, as well as (ii) Hv54h and Hv7Sh with a histidine-tagged Hv19h. These results provide the first in vivo evidence that these components interact in archaea. Such copurification studies will provide insight into the significance of the similarities and differences of the protein-targeting systems of the three domains of life, thereby increasing knowledge about the recognition of translocated proteins in general.

Gene content and organization of a 281-kbp contig from the genome of the extremely thermophilic archaeon, Sulfolobus solfataricus P2

The sequence of a 281-kbp contig from the crenarchaeote Sulfolobus solfataricus P2 was determined and analysed. Notable features in this region include 29 ribosomal protein genes, 12 tRNA genes (four of which contain archaeal-type introns), operons encoding enzymes of histidine biosynthesis, pyrimidine biosynthesis, and arginine biosynthesis, an ATPase operon, numerous genes for enzymes of lipopolysaccharide biosynthesis, and six insertion sequences. The content and organization of this contig are compared with sequences from crenarchaeotes, euryarchaeotes, bacteria, and eukaryotes.

Halophilic 20S proteasomes of the archaeon Haloferax volcanii: purification, characterization, and gene sequence analysis

A 20S proteasome, composed of alpha(1) and beta subunits arranged in a barrel-shaped structure of four stacked rings, was purified from a halophilic archaeon Haloferax volcanii. The predominant peptide-hydrolyzing activity of the 600-kDa alpha(1)beta-proteasome on synthetic substrates was cleavage carboxyl to hydrophobic residues (chymotrypsin-like [CL] activity) and was optimal at 2 M NaCl, pH 7.7 to 9.5, and 75 degrees C. The alpha(1)beta-proteasome also hydrolyzed insulin B-chain protein. Removal of NaCl inactivated the CL activity of the alpha(1)beta-proteasome and dissociated the complex into monomers. Rapid equilibration of the monomers into buffer containing 2 M NaCl facilitated their reassociation into fully active alpha(1)beta-proteasomes of 600 kDa. However, long-term incubation of the halophilic proteasome in the absence of salt resulted in hydrolysis and irreversible inactivation of the enzyme. Thus, the isolated proteasome has unusual salt requirements which distinguish it from any proteasome which has been described. Comparison of the beta-subunit protein sequence with the sequence deduced from the gene revealed that a 49-residue propeptide is removed to expose a highly conserved N-terminal threonine which is proposed to serve as the catalytic nucleophile and primary proton acceptor during peptide bond hydrolysis. Consistent with this mechanism, the known proteasome inhibitors carbobenzoxyl-leucinyl-leucinyl-leucinal-H (MG132) and N-acetyl-leucinyl-leucinyl-norleucinal (calpain inhibitor I) were found to inhibit the CL activity of the H. volcanii proteasome (K(i) = 0.2 and 8 microM, respectively). In addition to the genes encoding the alpha(1) and beta subunits, a gene encoding a second alpha-type proteasome protein (alpha(2)) was identified. All three genes coding for the proteasome subunits were mapped in the chromosome and found to be unlinked. Modification of the methods used to purify the alpha(1)beta-proteasome resulted in the copurification of the alpha(2) protein with the alpha(1) and beta subunits in nonstoichometric ratios as cylindrical particles of four stacked rings of 600 kDa with CL activity rates similar to the alpha(1)beta-proteasome, suggesting that at least two separate 20S proteasomes are synthesized. This study is the first description of a prokaryote which produces two separate 20S proteasomes and suggests that there may be distinct physiological roles for the two different alpha subunits in this halophilic archaeon.

Bacterial artificial chromosome (BAC) library as a tool for physical mapping of the archaeon Methanosarcina thermophila TM-1

We have used a variety of methods to characterize the genome of the archaeon Methanosarcina thermophila TM-1. Pulsed-field gel analysis indicates a genome size of 2.8 Mb. We have constructed a bacterial artificial chromosome (BAC) library of M. thermophila and have used it to generate physical maps for this organism. The library is made up of 384 clones with an average insert size of 58 kb representing 8.0 genome equivalents. The utility of the library for low-resolution physical mapping was shown by identifying NotI linking clones and using these to order the NotI macrorestriction fragments of M. thermophila into a 2.8 Mb map. Hybridization of nine single copy genes and a 16S rRNA sequence to these macrorestriction fragments forms the basis for the first genetic map in this organism. High-resolution physical maps, consisting of overlapping clones, have been created using HindIII fingerprints of BAC clones. In this way, we identified a minimal path of five clones that span a 270 kb NotI fragment. The ease of manipulating BAC clones makes the BAC system an excellent choice for the construction of low-resolution and high-resolution physical and genetic maps of archaeal genomes. It also provides a substrate for future genome-sequencing efforts.

Gene amplification and genomic plasticity in prokaryotes

Gene amplification is a common feature of the genome of prokaryotic organisms. In this review, we analyze different instances of gene amplification in a variety of prokaryotes, including their mechanisms of generation and biological role. Growing evidence supports the concept that gene amplification be considered not as a mutation but rather as a dynamic genomic state related to the adaptation of bacterial populations to changing environmental conditions or biological interactions. In this context, the potentially amplifiable DNA regions impose a defined dynamic structure on the genome. If such structure has indeed been selected during evolution, it is a particularly challenging hypothesis.

Structure and evolution of NGRRS-1, a complex, repeated element in the genome of Rhizobium sp. strain NGR234

Much of the remarkable ability of Rhizobium sp. strain NGR234 to nodulate at least 110 genera of legumes, as well as the nonlegume Parasponia andersonii, stems from the more than 80 different Nod factors it secretes. Except for nodE, nodG, and nodPQ, which are on the chromosome, most Nod factor biosynthesis genes are dispersed over the 536,165-bp symbiotic plasmid, pNGR234a. Mosaic sequences and insertion sequences (ISs) comprise 18% of pNGR234a. Many of them are clustered, and these IS islands divide the replicon into large blocks of functionally related genes. At 6 kb, NGRRS-1 is a striking example: there is one copy on pNGR234a and three others on the chromosome. DNA sequence comparisons of two NGRRS-1 elements identified three types of IS, NGRIS-2, NGRIS-4, and NGRIS-10. Here we show that all four copies of NGRRS-1 probably originated from transposition of NGRIS-4 into a more ancient IS-like sequence, NGRIS-10. Remarkably, all nine copies of NGRIS-4 have transposed into other ISs. It is unclear whether the accumulation of potentially mutagenic sequences in large clusters is due to the nature of the IS involved or to some selection process. Nevertheless, a direct consequence of the preferential targeting of transposons into such IS islands is to minimize the likelihood of disrupting vital functions.

Evolutionary analysis of the hisCGABdFDEHI gene cluster from the archaeon Sulfolobus solfataricus P2

While sequencing the genome of the archaeon Sulfolobus solfataricus P2, we found an 8,313-bp sequence containing a cluster of nine histidine biosynthesis genes in an order different from that of any known his operon. Results of phylogenetic analysis of the coding regions in the putative operon give conflicting evolutionary histories for individual his genes.

Organizational characteristics and information content of an archaeal genome: 156 kb of sequence from Sulfolobus solfataricus P2

We have initiated a project to sequence the 3 Mbp genome of the thermoacidophilic archaebacterium Sulfolobus solfataricus P2. Cosmids were selected from a provisional set of minimally overlapping clones, subcloned in pUC18, and sequenced using a hybrid (random plus directed) strategy to give two blocks of contiguous unique sequence, respectively, 100,389 and 56,105 bp. These two contigs contain a total of 163 open reading frames (ORFs) in 26-29 putative operons; 56 ORFs could be identified with reasonable certainty. Clusters of ORFs potentially encode proteins of glycogen biosynthesis, oxidative decarboxylation of pyruvate, ATP-dependent transport across membranes, isoprenoid biosynthesis, protein synthesis, and ribosomes. Putative promoters occur upstream of most ORFs. Thirty per cent of the predicted strong and medium-strength promoters can initiate transcription at the start codon or within 10 nucleotides upstream, indicating a process of initial mRNA-ribosome contact unlike that of most eubacterial genes. A novel termination motif is proposed to account for 15 additional terminations. The two contigs differ in densities of ORFs, insertion elements and repeated sequences; together they contain two copies of the previously reported insertion sequence ISC 1217, five additional IS elements representing four novel types, four classes of long non-IS repeated sequences, and numerous short, perfect repeats.

Comparative genomic analysis of the Haloferax volcanii DS2 and Halobacterium salinarium GRB contig maps reveals extensive rearrangement

Anonymous probes from the genome of Halobacterium salinarium GRB and 12 gene probes were hybridized to the cosmid clones representing the chromosome and plasmids of Halobacterium salinarium GRB and Haloferax volcanii DS2. The order of and pairwise distances between 35 loci uniquely cross-hybridizing to both chromosomes were analyzed in a search for conservation. No conservation between the genomes could be detected at the 15-kbp resolution used in this study. We found distinct sets of low-copy-number repeated sequences in the chromosome and plasmids of Halobacterium salinarium GRB, indicating some degree of partitioning between these replicons. We propose alternative courses for the evolution of the haloarchaeal genome: (i) that the majority of genomic differences that exist between genera came about at the inception of this group or (ii) that the differences have accumulated over the lifetime of the lineage. The strengths and limitations of investigating these models through comparative genomic studies are discussed.

Sizing chromosomes and megaplasmids in haloarchaea

PFGE was used for the genomic analysis of different species and strains belonging to four out of the six recognized haloarchaeal (halobacterial) genera. All of them were found to carry one chromosome from 1.8-3 Mb, and usually several, but at least one, large plasmids of approximately 90-680 kb, which were detected in supercoiled and linear forms. From the data gathered, chromosomal size appears to be conserved at genus level, whereas plasmid composition and size seems to be subjected to certain variability.

Differentially transcribed regions of Haloferax volcanii genome depending on the medium salinity

To identify genomic regions involved in osmoregulation in the extremely halophilic archaeon Haloferax volcanii, we used a technique which involves hybridization of cDNAs obtained at different salinities against a cosmid library of the organism. Both low and high salt concentrations trigger differential expression; however, adaptation to low salinities seems to elicit a wider response. The presence of a large domain within the largest of the megaplasmids with a strong response to low salt concentrations is noteworthy.

Evolution of the cell cycle

Cell proliferation involves duplication of all cell constituents and their more-or-less equal segregation to daughter cells. It seems probable that the performance of primitive cell-like structures would have been dogged by poor duplication and segregation fidelity, and by parasitism. This favoured evolution of the genome and with it the distinction between 'genomic' components like chromosomes whose synthesis is periodic and most other 'functional' components whose synthesis is continuous. Eukaryotic cells evolved from bacterial ancestors whose fused genome was replicated from a single origin and whose means of segregating sister chromatids depended on fixing their identity at replication. Evolution of an endo- or cytoskeleton, initially as means of consuming other bacteria, eventually enabled evolution of the mitotic spindle and a new means of segregating sister chromatids whose replication could be initiated from multiple origins. In this primitive eukaryotic cell, S and M phases might have been triggered by activation of a single cyclin-dependent kinase whose destruction along with that of other proteins would have triggered anaphase. Mitotic non-disjunction would have greatly facilitated genomic expansion, now possible due to multiple origins, and thereby accelerated the tempo of evolution when permitted by environmental conditions.

Genomic stability in the archaeae Haloferax volcanii and Haloferax mediterranei

Through hybridization of available probes, we have added nine genes to the macrorestriction map of the Haloferax mediterranei chromosome and five genes to the contig map of Haloferax volcanii. Additionally, we hybridized 17 of the mapped cosmid clones from H. volcanii to the H. mediterranei genome. The resulting 35-point chromosomal comparison revealed only two inversions and a few translocations. Forces known to promote rearrangement, common in the haloarchaea, have been ineffective in changing global gene order throughout the nearly 10(7) years of these species' divergent evolution.

Tempo, mode, the progenote, and the universal root

Early cellular evolution differed in both mode and tempo from the contemporary process. If modern lineages first began to diverge when the phenotype-genotype coupling was still poorly articulated, then we might be able to learn something about the evolution of that coupling through comparing the molecular biologies of living organisms. The issue is whether the last common ancestor of all life, the cenancestor, was a primitive entity, a progenote, with a more rudimentary genetic information-transfer system. Thinking on this issue is still unsettled. Much depends on the placement of the root of the universal tree and on whether or not lateral transfer renders such rooting meaningless.

Evolutionary relationships of bacterial and archaeal glutamine synthetase genes

Glutamine synthetase (GS), an essential enzyme in ammonia assimilation and glutamine biosynthesis, has three distinctive types: GSI, GSII and GSIII. Genes for GSI have been found only in bacteria (eubacteria) and archaea (archaebacteria), while GSII genes only occur in eukaryotes and a few soil-dwelling bacteria. GSIII genes have been found in only a few bacterial species. Recently, it has been suggested that several lateral gene transfers of archaeal GSI genes to bacteria may have occurred. In order to study the evolution of GS, we cloned and sequenced GSI genes from two divergent archaeal species: the extreme thermophile Pyrococcus furiosus and the extreme halophile Haloferax volcanii. Our phylogenetic analysis, which included most available GS sequences, revealed two significant prokaryotic GSI subdivisions: GSI-alpha and GSI-beta. GSI-alpha-genes are found in the thermophilic bacterium, Thermotoga maritima, the low G+C Gram-positive bacteria, and the Euryarchaeota (includes methanogens, halophiles, and some thermophiles). GSI-beta-type genes occur in all other bacteria. GSI-alpha- and GSI-beta-type genes also differ with respect to a specific 25-amino-acid insertion and adenylylation control of GS enzyme activity, both absent in the former but present in the latter. Cyanobacterial genes lack adenylylation regulation of GS and may have secondarily lost it. The GSI gene of Sulfolobus solfataricus, a member of the Crenarchaeota (extreme thermophiles), is exceptional and could not be definitely placed in either subdivision.

Physical map and set of overlapping cosmid clones representing the genome of the archaeon Halobacterium sp. GRB

We have constructed a complete, five-enzyme restriction map of the genome of the archaeon Halobacterium sp. GRB, based on a set of 84 overlapping cosmid clones. Fewer than 30 kbp, in three gaps, remain uncloned. The genome consists of five replicons: a chromosome (2038 kbp) and four plasmids (305, 90, 37, and 1.8 kbp). The genome of Halobacterium sp. GRB is similar in style to other halobacterial genomes by being partitioned among multiple replicons and by being mosaic in terms of nucleotide composition. It is unlike other halobacterial genomes, however, in lacking multicopy families of insertion sequences.

Chromosomal structure of Rhodobacter capsulatus strain SB1003: cosmid encyclopedia and high-resolution physical and genetic map

A combination of cosmid genome walking and pulsed-field gel electrophoresis was used to construct a high-resolution physical and genetic map of the 3.8-megabase (Mb) genome of Rhodobacter capsulatus SB1003. The mapping was done by hybridization of pulsed-field gel blots and by grouping and further mapping of the cosmids and bacteriophages from genomic libraries. Cosmid clones formed two uninterrupted and ordered groups, one corresponding to the chromosome of R. capsulatus, the other to its 134-kb plasmid. Cos site end-labeling and partial EcoRV digestion of cosmids were used to construct a high-resolution EcoRV map of the genome. Overlapping of the cosmids was confirmed by the resemblance of the cosmid restriction maps and by direct end-to-end hybridization with SP6- and T7-specific transcripts. Twenty-three previously cloned genes and eight groups of repeated sequences, revealed in this work, were located in the ordered gene library and mapped with an accuracy of 1-10 kb. Blots of a minimal set of 192 cosmids, covering the chromosome and the plasmid with the known map position of each cosmid, give to R. capsulatus the same advantages that the Kohara phage panel gives to E. coli.

Physical map of the Methanobacterium thermoautotrophicum Marburg chromosome

A physical map of the Methanobacterium thermoautotrophicum Marburg chromosome was constructed by using pulsed-field gel electrophoresis of restriction fragments generated by NotI, PmeI, and NheI. The order of the fragments was deduced from Southern blot hybridization of NotI fragment probes to various restriction digests and from partial digests. The derived map is circular, and the genome size was estimated to be 1,623 kb. Several cloned genes were hybridized to restriction fragments to locate their positions on the map. Genes coding for proteins involved in the methanogenic pathway were located on the same segment of the circular chromosome. In addition, the genomes of a variety of thermophilic Methanobacterium strains were treated with restriction enzymes and analyzed by pulsed-field gel electrophoresis. The sums of the fragment sizes varied from 1,600 to 1,728 kb among the strains, and widely different macrorestriction patterns were observed.