Single replication origin of the archaeon Methanosarcina mazei revealed by the Z curve method

The genomic sequence of the archaeon Methanosarcina mazei has been analyzed by the Z curve method. The Z curve is a three-dimensional curve that uniquely represents the given DNA sequence. The three-dimensional Z curve and its x and y components for the genome of M. mazei show a sharp peak and relatively broad peak, respectively. The cdc6 gene is located exactly at the position of the sharp peak. Based on the known behavior of the Z curves for the archaea whose replication origins have been identified, we hypothesize that the replication origin and termination sites correspond to the positions of the sharp peak and broad peak, respectively. We have located an intergenic region that is between the cdc6 gene (MM1314) and the gene for an adjacent protein (MM1315), which shows strong characteristics of the known replication origins. This region is highly rich in AT and contains multiple copies of consecutive repeats. Our results strongly suggest that the single replication origin of M. mazei is situated at the intergenic region between the cdc6 gene and the gene for the adjacent protein, from 1,564,657 to 1,566,241 bp of the genome.     

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

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

Identification of replication origins in the genome of the methanogenic archaeon,<Emphasis Type="Italic"> Methanocaldococcus jannaschii</Emphasis>

Methanocaldococcus jannaschii has been notorious as an archaeon in which the replication origins are difficult to identify. Although extensive efforts have been exerted on this issue, the locations of replication origins still remain elusive 7 years after the publication of its complete genome sequence in 1996. Ambiguous results were obtained in identifying the replication origins of M. jannaschii based on all theoretical and experimental approaches. In the genome of M. jannaschii, we found that an ORF (MJ0774), annotated as a hypothetical protein, is a homologue of the Cdc6 protein. The position of the gene is at a global minimum of the x component of the Z curve, i.e., RY disparity curve, which has been used to identify replication origins in other Archaea. In addition, an intergenic region (694,540–695,226 bp) that is between the cdc6 gene and an adjacent ORF shows almost all the characteristics of known replication origins, i.e., it is highly rich in AT composition (80%) and contains multiple copies of repeat elements and AT stretches. Therefore, these lines of evidence strongly suggest that the identified region is a replication origin, which is designated as oriC1. The analysis of the y component of the Z curve, i.e., MK disparity curve, suggests the presence of another replication origin corresponding to one of the peaks in the MK disparity curve at around 1,388 kb of the genome.Communicated by G. Antranikian     

The cobY gene of the archaeon Halobacterium sp. strain NRC-1 is required for de novo cobamide synthesis

Genetic and nutritional analyses of mutants of the extremely halophilic archaeon Halobacterium sp. strain NRC-1 showed that open reading frame (ORF) Vng1581C encodes a protein with nucleoside triphosphate:adenosylcobinamide-phosphate nucleotidyltransferase enzyme activity. This activity was previously associated with the cobY gene of the methanogenic archaeon Methanobacterium thermoautotrophicum strain DeltaH, but no evidence was obtained to demonstrate the direct involvement of this protein in cobamide biosynthesis in archaea. Computer analysis of the Halobacterium sp. strain NRC-1 ORF Vng1581C gene and the cobY gene of M. thermoautotrophicum strain DeltaH showed the primary amino acid sequence of the proteins encoded by these two genes to be 35% identical and 48% similar. A strain of Halobacterium sp. strain NRC-1 carrying a null allele of the cobY gene was auxotrophic for cobinamide-GDP, a known intermediate of the late steps of cobamide biosynthesis. The auxotrophic requirement for cobinamide-GDP was corrected when a wild-type allele of cobY was introduced into the mutant strain, demonstrating that the lack of cobY function was solely responsible for the observed block in cobamide biosynthesis in this archaeon. The data also show that Halobacterium sp. strain NRC-1 possesses a high-affinity transport system for corrinoids and that this archaeon can synthesize cobamides de novo under aerobic growth conditions. To the best of our knowledge this is the first genetic and nutritional analysis of cobalamin biosynthetic mutants in archaea.     

Recent Advances in the Identification of Replication Origins Based on the Z-curve Method

Precise DNA replication is critical for the maintenance of genetic integrity in all organisms. In all three domains of life, DNA replication starts at a specialized locus, termed as the replication origin, oriC or ORI, and its identification is vital to understanding the complex replication process. In bacteria and eukaryotes, replication initiates from single and multiple origins, respectively, while archaea can adopt either of the two modes. The Z-curve method has been successfully used to identify replication origins in genomes of various species, including multiple oriCs in some archaea. Based on the Z-curve method and comparative genomics analysis, we have developed a web-based system, Ori-Finder, for finding oriCs in bacterial genomes with high accuracy. Predicted oriC regions in bacterial genomes are organized into an online database, DoriC. Recently, archaeal oriC regions identified by both in vivo and in silico methods have also been included in the database. Here, we summarize the recent advances of in silico prediction of oriCs in bacterial and archaeal genomes using the Z-curve based method.     

Multiple replication origins with diverse control mechanisms in Haloarcula hispanica

The use of multiple replication origins in archaea is not well understood. In particular, little is known about their specific control mechanisms. Here, we investigated the active replication origins in the three replicons of a halophilic archaeon, Haloarcula hispanica, by extensive gene deletion, DNA mutation and genome-wide marker frequency analyses. We revealed that individual origins are specifically dependent on their co-located cdc6 genes, and a single active origin/cdc6 pairing is essential and sufficient for each replicon. Notably, we demonstrated that the activities of oriC1 and oriC2, the two origins on the main chromosome, are differently controlled. A G-rich inverted repeat located in the internal region between the two inverted origin recognition boxes (ORBs) plays as an enhancer for oriC1, whereas the replication initiation at oriC2 is negatively regulated by an ORB-rich region located downstream of oriC2-cdc6E, likely via Cdc6E-titrating. The oriC2 placed on a plasmid is incompatible with the wild-type (but not the ΔoriC2) host strain, further indicating that strict control of the oriC2 activity is important for the cell. This is the first report revealing diverse control mechanisms of origins in haloarchaea, which has provided novel insights into the use and coordination of multiple replication origins in the domain of Archaea.     

Multiple origins of replication in archaea

Until recently, the only archaeon for which a bona fide origin of replication was reported was Pyrococcus abyssi, where a single origin was identified. Although several in silico analyses have suggested that some archaeal species might contain more than one origin, this has only been demonstrated recently. Two studies have shown that multiple origins of replication function in two archaeal species. One study identified two origins of replication in the archaeon Sulfolobus solfataricus, whereas a second study used a different technique to show that both S. solfataricus and Sulfolobus acidocaldarius have three functional origins. These are the first reports of archaea having multiple origins. This finding has implications for research on the mechanisms of DNA replication and evolution.     

Cloning and characterization of the chromosomal replication origin region of Amycolatopsis mediterranei U32

The chromosomal replication origins (oriC) of gram positive, acid-fast actinomycetes have been investigated in streptomycetes and mycobacteria. A 1339 bp DNA fragment of the putative oriC region from the rifamycin SV producer Amycolatopsis mediterranei U32 was cloned by PCR amplification employing primers designed based on the conserved flanking genes of dnaA and dnaN. The 884 bp sequence of the intergenic region between dnaA and dnaN genes consists of 19 DnaA-boxes and two 13-mer AT-rich sequences, which is similar to the oriC structure of Streptomyces lividans. A mini-chromosome constructed by cloning the putative U32 oriC DNA fragment into an Escherichia coli plasmid was able to replicate autonomously, but was unstable, in A. mediterranei U32 with an estimated copy number of two per cell. Although efficient replication of the mini-chromosome in U32 requires the complete set of DnaA-boxes and AT-rich regions, only one of the AT-rich sequences together with part of the DnaA-boxes is sufficient, suggesting the presence of combinatorial alternatives for a functional oriC region of A. mediterranei U32. Phylogenetic analysis based on definite oriC sequences among eubacteria reflects well the relationship between these species.     

Poorly conserved ORFs in the genome of the archaea Halobacterium sp. NRC-1 correspond to expressed proteins

MOTIVATION: A large fraction of open reading frames (ORFs) identified as 'hypothetical' proteins correspond to either 'conserved hypothetical' proteins, representing sequences homologous to ORFs of unknown function from other organisms, or to hypothetical proteins lacking any significant sequence similarity to other ORFs in the databases. Elucidating the functions and three-dimensional structures of such orphan ORFs, termed ORFans or poorly conserved ORFs (PCOs), is essential for understanding biodiversity. However, it has been claimed that many ORFans may not encode for expressed proteins. RESULTS: A genome-wide experimental study of 'paralogous PCOs' in the halophilic archaea Halobacterium sp. NRC-1 was conducted. Paralogous PCOs are ORFs with at least one homolog in the same organism, but with no clear homologs in other organisms. The results reveal that mRNA is synthesized for a majority of the Halobacterium sp. NRC-1 paralogous PCO families, including those comprising relatively short proteins, strongly suggesting that these Halobacterium sp. NRC-1 paralogous PCOs correspond to true, expressed proteins. Hence, further computational and experimental studies aimed at characterizing PCOs in this and other organisms are merited. Such efforts could shed light on PCOs' functions and origins, thereby serving to elucidate the vast diversity observed in the genetic material.     

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.Archaea are of considerable interest due to their unusual phylogenetic position and the similarity of their information transfer system to that of eukaryotes. In particular, studies of DNA replication in archaea have revealed characteristics of both bacterial and eukaryotic systems (1). While genome sequencing has shown that archaeal and bacterial genomes are composed of a single or few circular chromosomes, comparative genomic studies have found that most components of the archaeal DNA replication machinery, such as the origin recognition proteins, DNA polymerases, helicases, and primases, are similar to eukaryotic proteins. The hybrid nature of archaeal DNA replication systems raises important questions regarding the mechanism by which they select an origin(s) for initiation and coordinate orderly DNA replication and segregation into daughter cells.Our understanding of DNA replication in archaea has thus far been based primarily on bioinformatic studies, with experimental analysis restricted to only a few tractable systems. An initial study of Pyrococcus species using GC (tetramer) skew analysis suggested that they use a single, unique origin of replication in their chromosomes. Subsequent [³H]uracil labeling analysis of Pyrococcus abyssi (21) showed that newly synthesized DNA mapped to the predicted replication origin region, which contained the only orc gene in the genome, a D family DNA polymerase gene, and a DNA sliding clamp loader subunit. In addition, two-dimensional gel analysis of replicating molecules confirmed the location of the DNA replication origin near the orc1 gene of P. abyssi, with predicted origin binding sequences and AT-rich DNA unwinding elements nearby (18). An investigation of DNA replication in Aeropyrum pernix used a combination of biochemical and two-dimensional gel electrophoresis and identified two potential sites of replication initiation, on opposite sides of the circular genome (14, 28). One of these sites (oriC1_Ap) contained four origin recognition boxes and an AT-rich region and was shown to be bound by the ORC1 gene. The other site (oriC2_Ap) contained repeat elements without an intervening AT-rich region and has been shown by two-dimensional gel electrophoresis to contain an active replication origin (28). An examination of replication in two Sulfolobus spp., Sulfolobus solfataricus and Sulfolobus acidocaldarius (16, 30), by use of a combination of bioinformatic and two-dimensional gel analysis and of marker frequency by use of DNA microarrays identified three well-separated replication origins per genome. Only two of the three origins were originally identified, due to their linkage to orc genes and conserved origin binding sequences, while the third was identified by marker frequency analysis (MFA). Using partially synchronized cells of S. acidocaldarius, the origins were shown to initiate DNA replication synchronously, indicating a highly coordinated and regulated process. Biochemical analysis has shown that either two or all three Orc proteins are able to bind to all Sulfolobus origins; however, binding at the third origin is considerably weaker (29). Replication origins were also recently identified in Methanothermobacter thermoautotrophicus (17).Our laboratory has been investigating DNA replication in a halophilic archaeon capable of growth at saturating NaCl concentrations. The model system, Halobacterium sp. strain NRC-1, was one of the earliest archaeal genomes to be sequenced (23) and provided a DNA knockout method, utilizing the selectable and counterselectable ura3 gene, for genetic analysis (25). The NRC-1 genome was found to be organized into a 2-Mbp chromosome and two large and partially redundant extrachromosomal elements, pNRC100 and pNRC200. The genome sequence showed that the orc gene family was highly expanded, with four genes (orc6, -7, -8, and -10) distributed in the chromosome and six genes (orc1, -2, -3, -4, -5, and -9) in pNRC200, one of which (orc9) was also present in pNRC100. Three rep genes thought to be important for replication initiation were present in one (repJ in pNRC100) or both (repH and repI) of the extrachromosomal elements. Regions near two of these genes, orc7 and repH, were shown to harbor autonomous replicating ability and to contain inverted repeat sequences (IRs) and an AT-rich presumptive DNA unwinding region detectable by χ² analysis (3, 22). Additionally, GC/oligomer skew analyses of Halobacterium sp. strain NRC-1 showed multiple inflection points in the chromosome, suggestive of multiple replication origins in this strain (15, 34).Halobacterium sp. strain NRC-1 is the only archaeal system where gene mutation analysis has established which predicted DNA replication genes are essential to cells (2). As expected, two DNA polymerases (one B family and one D family polymerase), the MCM DNA helicase, DNA primase (Pri1/Pri2), the sliding clamp (PCNA), and flap endonuclease (Rad2) were all found to be essential. However, one B family DNA polymerase gene and 8 of the 10 orc and cdc6 genes, including the orc7 gene, were found to be nonessential by deletion analysis. Only the orc2 gene in pNRC200 and the orc10 gene in the chromosome were found to be essential, suggesting a critical role(s) for these genes in DNA replication.In this study, we used a combination of MFA, employing whole-genome DNA microarrays, the ura3-based gene knockout method, and replication initiation point (RIP) analysis to further investigate DNA replication in Halobacterium sp. strain NRC-1. Our results indicate that initiation of DNA replication in NRC-1 is more complex than originally anticipated, with multiple origins likely present on the chromosome and the extrachromosomal elements.     

A new pathway for salvaging the coenzyme B12 precursor cobinamide in archaea requires cobinamide-phosphate synthase (CbiB) enzyme activity

The ability of archaea to salvage cobinamide has been under question because archaeal genomes lack orthologs to the bacterial nucleoside triphosphate:5'-deoxycobinamide kinase enzyme (cobU in Salmonella enterica). The latter activity is required for cobinamide salvaging in bacteria. This paper reports evidence that archaea salvage cobinamide from the environment by using a pathway different from the one used by bacteria. These studies demanded the functional characterization of two genes whose putative function had been annotated based solely on their homology to the bacterial genes encoding adenosylcobyric acid and adenosylcobinamide-phosphate synthases (cbiP and cbiB, respectively) of S. enterica. A cbiP mutant strain of the archaeon Halobacterium sp. strain NRC-1 was auxotrophic for adenosylcobyric acid, a known intermediate of the de novo cobamide biosynthesis pathway, but efficiently salvaged cobinamide from the environment, suggesting the existence of a salvaging pathway in this archaeon. A cbiB mutant strain of Halobacterium was auxotrophic for adenosylcobinamide-GDP, a known de novo intermediate, and did not salvage cobinamide. The results of the nutritional analyses of the cbiP and cbiB mutants suggested that the entry point for cobinamide salvaging is adenosylcobyric acid. The data are consistent with a salvaging pathway for cobinamide in which an amidohydrolase enzyme cleaves off the aminopropanol moiety of adenosylcobinamide to yield adenosylcobyric acid, which is converted by the adenosylcobinamide-phosphate synthase enzyme to adenosylcobinamide-phosphate, a known intermediate of the de novo biosynthetic pathway. The existence of an adenosylcobinamide amidohydrolase enzyme would explain the lack of an adenosylcobinamide kinase in archaea.     

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

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

MutS and MutL Are Dispensable for Maintenance of the Genomic Mutation Rate in the Halophilic Archaeon Halobacterium salinarum NRC-1

Background

The genome of the halophilic archaeon Halobacterium salinarum NRC-1 encodes for homologs of MutS and MutL, which are key proteins of a DNA mismatch repair pathway conserved in Bacteria and Eukarya. Mismatch repair is essential for retaining the fidelity of genetic information and defects in this pathway result in the deleterious accumulation of mutations and in hereditary diseases in humans.

Methodology/Principal Findings

We calculated the spontaneous genomic mutation rate of H. salinarum NRC-1 using fluctuation tests targeting genes of the uracil monophosphate biosynthesis pathway. We found that H. salinarum NRC-1 has a low incidence of mutation suggesting the presence of active mechanisms to control spontaneous mutations during replication. The spectrum of mutational changes found in H. salinarum NRC-1, and in other archaea, appears to be unique to this domain of life and might be a consequence of their adaption to extreme environmental conditions. In-frame targeted gene deletions of H. salinarum NRC-1 mismatch repair genes and phenotypic characterization of the mutants demonstrated that the mutS and mutL genes are not required for maintenance of the observed mutation rate.

Conclusions/Significance

We established that H. salinarum NRC-1 mutS and mutL genes are redundant to an alternative system that limits spontaneous mutation in this organism. This finding leads to the puzzling question of what mechanism is responsible for maintenance of the low genomic mutation rates observed in the Archaea, which for the most part do not have MutS and MutL homologs.     

The Z curve database: a graphic representation of genome sequences

MOTIVATION: Genome projects for many prokaryotic and eukaryotic species have been completed and more new genome projects are being underway currently. The availability of a large number of genomic sequences for researchers creates a need to find graphic tools to study genomes in a perceivable form. The Z curve is one of such tools available for visualizing genomes. The Z curve is a unique three-dimensional curve representation for a given DNA sequence in the sense that each can be uniquely reconstructed given the other. The Z curve database for more than 1000 genomes have been established here. RESULTS: The database contains the Z curves for archaea, bacteria, eukaryota, organelles, phages, plasmids, viroids and viruses, whose genomic sequences are currently available. All the 3-dimensional Z curves and their three component curves are stored in the database. The applications of the Z curve database on comparative genomics, gene prediction, computation of G+C content with a windowless technique, prediction of replication origins and terminations of bacterial and archaeal genomes and study of local deviations from the Chargaff Parity Rule 2 etc. are presented in detail. The Z curve database reported here is a treasure trove in which biologists could find useful biological knowledge.     

Proteomics of Halophilic archaea

Halophilic archaea is a member of the Halobacteriacea family, the only family in the Halobacteriales order. Most Halophilic archaea require 1.5M NaCl both to grow and retain the structural integrity of the cells. The proteins of these organisms have thus been adapted to be active and stable in the hypersaline condition. Consequently, the unique properties of these biocatalysts have resulted in several novel applications in industrial processes. Halophilic archaea are also to be useful for bioremediation of hypersaline environment. Proteome data have expended enormously with the significant advance recently achieved in two-dimensional gel electrophoresis (2-DE) and mass spectrometry (MS). The whole genome sequencing of Halobacterium species NRC-1 was completed and this would also provide tremendous help to analyze the protein mass data from the similar strain Halobacterium salinarum. Proteomics coupled with genomic databases now has become a basic tool to understand or identify the function of genes and proteins. In addition, the bioinformatics approach will facilitate to predict the function of novel proteins of Halophilic archaea. This review will discuss current proteome study of Halophilic archaea and introduce the efficient procedures for screening, predicting, and confirming the function of novel halophilic enzymes.