Identification of replication origins in the genome of the methanogenic archaeon, Methanocaldococcus jannaschii

Diversity of DNA Replication in the Archaea

DNA replication is arguably the most fundamental biological process. On account of their shared evolutionary ancestry, the replication machinery found in archaea is similar to that found in eukaryotes. DNA replication is initiated at origins and is highly conserved in eukaryotes, but our limited understanding of archaea has uncovered a wide diversity of replication initiation mechanisms. Archaeal origins are sequence-based, as in bacteria, but are bound by initiator proteins that share homology with the eukaryotic origin recognition complex subunit Orc1 and helicase loader Cdc6). Unlike bacteria, archaea may have multiple origins per chromosome and multiple Orc1/Cdc6 initiator proteins. There is no consensus on how these archaeal origins are recognised—some are bound by a single Orc1/Cdc6 protein while others require a multi- Orc1/Cdc6 complex. Many archaeal genomes consist of multiple parts—the main chromosome plus several megaplasmids—and in polyploid species these parts are present in multiple copies. This poses a challenge to the regulation of DNA replication. However, one archaeal species (Haloferax volcanii) can survive without replication origins; instead, it uses homologous recombination as an alternative mechanism of initiation. This diversity in DNA replication initiation is all the more remarkable for having been discovered in only three groups of archaea where in vivo studies are possible.

Intergenomic evolution and metabolic cross-talk between rumen and thermophilic autotrophic methanogenic archaea

Methanobrevibacter ruminantium M1 (MRU) is a rumen methanogenic archaean that can be able to utilize formate and CO₂/H₂ as growth substrates. Extensive analysis on the evolutionary genomic contexts considered herein to unravel its intergenomic relationship and metabolic adjustment acquired from the genomic content of Methanothermobacter thermautotrophicus ΔH. We demonstrated its intergenomic distance, genome function, synteny homologs and gene families, origin of replication, and methanogenesis to reveal the evolutionary relationships between Methanobrevibacter and Methanothermobacter. Comparison of the phylogenetic and metabolic markers was suggested for its archaeal metabolic core lineage that might have evolved from Methanothermobacter. Orthologous genes involved in its hydrogenotrophic methanogenesis might be acquired from intergenomic ancestry of Methanothermobacter via Methanobacterium formicicum. Formate dehydrogenase (fdhAB) coding gene cluster and carbon monoxide dehydrogenase (cooF) coding gene might have evolved from duplication events within Methanobrevibacter-Methanothermobacter lineage, and fdhCD gene cluster acquired from bacterial origins. Genome-wide metabolic survey found the existence of four novel pathways viz. l-tyrosine catabolism, mevalonate pathway II, acyl-carrier protein metabolism II and glutathione redox reactions II in MRU. Finding of these pathways suggested that MRU has shown a metabolic potential to tolerate molecular oxygen, antimicrobial metabolite biosynthesis and atypical lipid composition in cell wall, which was acquainted by metabolic cross-talk with mammalian bacterial origins. We conclude that coevolution of genomic contents between Methanobrevibacter and Methanothermobacter provides a clue to understand the metabolic adaptation of MRU in the rumen at different environmental niches.

Global phylogenomic analysis disentangles the complex evolutionary history of DNA replication in archaea

The archaeal machinery responsible for DNA replication is largely homologous to that of eukaryotes and is clearly distinct from its bacterial counterpart. Moreover, it shows high diversity in the various archaeal lineages, including different sets of components, heterogeneous taxonomic distribution, and a large number of additional copies that are sometimes highly divergent. This has made the evolutionary history of this cellular system particularly challenging to dissect. Here, we have carried out an exhaustive identification of homologs of all major replication components in over 140 complete archaeal genomes. Phylogenomic analysis allowed assigning them to either a conserved and probably essential core of replication components that were mainly vertically inherited, or to a variable and highly divergent shell of extra copies that have likely arisen from integrative elements. This suggests that replication proteins are frequently exchanged between extrachromosomal elements and cellular genomes. Our study allowed clarifying the history that shaped this key cellular process (ancestral components, horizontal gene transfers, and gene losses), providing important evolutionary and functional information. Finally, our precise identification of core components permitted to show that the phylogenetic signal carried by DNA replication is highly consistent with that harbored by two other key informational machineries (translation and transcription), strengthening the existence of a robust organismal tree for the 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.

A Brief Review: The Z-curve Theory and its Application in Genome Analysis

In theoretical physics, there exist two basic mathematical approaches, algebraic and geometrical methods, which, in most cases, are complementary. In the area of genome sequence analysis, however, algebraic approaches have been widely used, while geometrical approaches have been less explored for a long time. The Z-curve theory is a geometrical approach to genome analysis. The Z-curve is a three-dimensional curve that represents a given DNA sequence in the sense that each can be uniquely reconstructed given the other. The Z-curve, therefore, contains all the information that the corresponding DNA sequence carries. The analysis of a DNA sequence can then be performed through studying the corresponding Z-curve. The Z-curve method has found applications in a wide range of areas in the past two decades, including the identifications of protein-coding genes, replication origins, horizontally-transferred genomic islands, promoters, translational start sides and isochores, as well as studies on phylogenetics, genome visualization and comparative genomics. Here, we review the progress of Z-curve studies from aspects of both theory and applications in genome analysis.

Mapping of active replication origins in vivo in thaum- and euryarchaeal replicons

We report mapping of active replication origins in thaum- and euryarchaeal replicons using high-throughput sequencing-based marker frequency analysis. The chromosome of the thaumarchaeon Nitrosopumilus maritimus is shown to contain a single origin of replication, whereas the main chromosome in the halophilic euryarchaea Haloferax mediterranei and Haloferax volcanii each contains two origins. All replication origins specified bidirectional replication, and the two origins in the halophiles were initiated in synchrony. The pHM500 plasmid of H. mediterranei is shown to contain a single origin, and the copy numbers of five plasmid replicons in the two halophiles were inferred to be close to that of the main chromosome. Origin recognition boxes (ORBs) that provide binding sites for Orc1/Cdc6 replication initiator proteins are identified at all chromosomal origins, as well as in a range of additional thaumarchaeal species. An annotation update is provided for all three species.

DoriC 5.0: an updated database of oriC regions in both bacterial and archaeal genomes

Replication of chromosomes is one of the central events in the cell cycle. Chromosome replication begins at specific sites, called origins of replication (oriCs), for all three domains of life. However, the origins of replication still remain unknown in a considerably large number of bacterial and archaeal genomes completely sequenced so far. The availability of increasing complete bacterial and archaeal genomes has created challenges and opportunities for identification of their oriCs in silico, as well as in vivo. Based on the Z-curve theory, we have developed a web-based system Ori-Finder to predict oriCs in bacterial genomes with high accuracy and reliability by taking advantage of comparative genomics, and the predicted oriC regions have been organized into an online database DoriC, which is publicly available at http://tubic.tju.edu.cn/doric/ since 2007. Five years after we constructed DoriC, the database has significant advances over the number of bacterial genomes, increasing about 4-fold. Additionally, oriC regions in archaeal genomes identified by in vivo experiments, as well as in silico analyses, have also been added to the database. Consequently, the latest release of DoriC contains oriCs for >1500 bacterial genomes and 81 archaeal genomes, respectively.

Nucleotide correlation based measure for identifying origin of replication in genomic sequences

Computational prediction of the origin of replication is a challenging problem and of immense interest to biologists. Several methods have been proposed for identifying the replicon site for various classes of organisms. However, these methods have limited applicability since the replication mechanism is different in different organisms. We propose a correlation measure and show that it is correctly able to predict the origin of replication in most of the bacterial genomes. When applied to Methanocaldococcus jannaschii, Plasmodium falciparum apicoplast and Nicotiana tabacum plastid, this correlation based method is able to correctly predict the origin of replication whereas the generally used GC skew measure fails. Thus, this correlation based measure is a novel and promising tool for predicting the origin of replication in a wide class of organisms. This could have important implications in not only gaining a deeper understanding of the replication machinery in higher organisms, but also for drug discovery.

Evolutionary rates and gene dispensability associate with replication timing in the archaeon Sulfolobus islandicus

In bacterial chromosomes, the position of a gene relative to the single origin of replication generally reflects its replication timing, how often it is expressed, and consequently, its rate of evolution. However, because some archaeal genomes contain multiple origins of replication, bias in gene dosage caused by delayed replication should be minimized and hence the substitution rate of genes should associate less with chromosome position. To test this hypothesis, six archaeal genomes from the genus Sulfolobus containing three origins of replication were selected, conserved orthologs were identified, and the evolutionary rates (dN and dS) of these orthologs were quantified. Ortholog families were grouped by their consensus position and designated by their proximity to one of the three origins (O1, O2, O3). Conserved orthologs were concentrated near the origins and most variation in genome content occurred distant from the origins. Linear regressions of both synonymous and nonsynonymous substitution rates on distance from replication origins were significantly positive, the rates being greatest in the region furthest from any of the origins and slowest among genes near the origins. Genes near O1 also evolved faster than those near O2 and O3, which suggest that this origin may fire later in the cell cycle. Increased evolutionary rates and gene dispensability are strongly associated with reduced gene expression caused in part by reduced gene dosage during the cell cycle. Therefore, in this genus of Archaea as well as in many Bacteria, evolutionary rates and variation in genome content associate with replication timing.

The evolutionary history of archaeal MCM helicases: a case study of vertical evolution combined with hitchhiking of mobile genetic elements

Genes encoding DNA replication proteins have been frequently exchanged between cells and mobile elements, such as viruses or plasmids. This raises potential problems to reconstruct their history. Here, we combine phylogenetic and genomic context analyses to study the evolution of the replicative minichromosome maintenance (MCM) helicases in Archaea. Several archaeal genomes encode more than one copy of the mcm gene. Genome context analysis reveals that most of these additional copies are encoded within mobile elements. Exhaustive analysis of these elements reveals diverse groups of integrated archaeal plasmids or viruses, including several head-and-tail proviruses. Some MCMs encoded by mobile elements are structurally distinct from their cellular counterparts, with one case of novel domain organization. Both genome context and phylogenetic analysis indicate that MCM encoded by mobile elements were recruited from cellular genomes. An accelerated evolution and a dramatic expansion of methanococcal MCMs suggest a host-to-virus-to-host transfer loop, possibly triggered by the loss of the archaeal initiator protein Cdc6 in Methanococcales. Surprisingly, despite extensive transfer of mcm genes between viruses, plasmids, and cells, the topology of the MCM tree is strikingly congruent with the consensus archaeal phylogeny, indicating that mobile elements encoding mcm have coevolved with their hosts and that DNA replication proteins can be also useful to reconstruct the history of the archaeal domain.

Localized melting of duplex DNA by Cdc6/Orc1 at the DNA replication origin in the hyperthermophilic archaeon Pyrococcus furiosus

The initiation step is a key process to regulate the frequency of DNA replication. Although recent studies in Archaea defined the origin of DNA replication (oriC) and the Cdc6/Orc1 homolog as an origin recognition protein, the location and mechanism of duplex opening have remained unclear. We have found that Cdc6/Orc1 binds to oriC and unwinds duplex DNA in the hyperthermophilic archaeon Pyrococcus furiosus, by means of a P1 endonuclease assay. A primer extension analysis further revealed that this localized unwinding occurs in the oriC region at a specific site, which is 12-bp long and rich in adenine and thymine. This site is different from the predicted duplex unwinding element (DUE) that we reported previously. We also discovered that Cdc6/Orc1 induces topological changes in supercoiled oriC DNA, and that this process is dependent on the AAA+ domain. These results indicate that topological alterations of oriC DNA by Cdc6/Orc1 introduce a single-stranded region at the 12-mer site, that could possibly serve as an entry point for Mcm helicase.

The relaxing ori-ter balance of Mycoplasma genomes

Mycoplasma are wall-less bacteria with small genomes, which are thought to have resulted from massive genome reductive processes, during which the ori-ter balance may be disrupted. For technical difficulties, ori and ter have been located only in a few Mycoplasma strains. Using the Z curve method, we were able to locate turning points on the Mycoplasma genomes, with the minimum and maximum points co-locating with ori or ter in the reference genomes. Assuming Z curve correctly located ori and ter, we calculated the distances from ori to ter in both directions on the circular genome and calculated the ori-ter balance status. The Mycoplasma genomes were not balanced, possibly as a result of close association of Mycoplasma with hosts, where there would be no other microbes for Mycoplasma to compete with for nutrients, so fastest possible growth related to balanced genomes might not be needed by Mycoplasma, leading to a relaxing ori-ter balance.

Genomic context analysis in Archaea suggests previously unrecognized links between DNA replication and translation

Specific functional interactions of proteins involved in DNA replication and/or DNA repair or transcription might occur in Archaea, suggesting a previously unrecognized regulatory network coupling DNA replication and translation, which might also exist in Eukarya.

Background

Comparative analysis of genomes is valuable to explore evolution of genomes, deduce gene functions, or predict functional linking between proteins. Here, we have systematically analyzed the genomic environment of all known DNA replication genes in 27 archaeal genomes to infer new connections for DNA replication proteins from conserved genomic associations.

Results

Two distinct sets of DNA replication genes frequently co-localize in archaeal genomes: the first includes the genes for PCNA, the small subunit of the DNA primase (PriS), and Gins15; the second comprises the genes for MCM and Gins23. Other genomic associations of genes encoding proteins involved in informational processes that may be functionally relevant at the cellular level have also been noted; in particular, the association between the genes for PCNA, transcription factor S, and NudF. Surprisingly, a conserved cluster of genes coding for proteins involved in translation or ribosome biogenesis (S27E, L44E, aIF-2 alpha, Nop10) is almost systematically contiguous to the group of genes coding for PCNA, PriS, and Gins15. The functional relevance of this cluster encoding proteins conserved in Archaea and Eukarya is strongly supported by statistical analysis. Interestingly, the gene encoding the S27E protein, also known as metallopanstimulin 1 (MPS-1) in human, is overexpressed in multiple cancer cell lines.

Conclusion

Our genome context analysis suggests specific functional interactions for proteins involved in DNA replication between each other or with proteins involved in DNA repair or transcription. Furthermore, it suggests a previously unrecognized regulatory network coupling DNA replication and translation in Archaea that may also exist in Eukarya.

An easier method for locating replication origins based on real-time polymerase chain reaction

Locating replication origins (ROs) is an old problem to which diverse answers have been previously given. However, a better and more convenient method still remains desirable because previously reported methods are time-consuming. In this study, we report an easier method for locating ROs based on real-time polymerase chain reaction. The feasibility of this method was demonstrated using Escherichia coli K-12 and Halobacterium salinarium NRC-1 as models. Our results show that this method is easy, fast, and accurate in locating ROs.

A conserved mechanism for replication origin recognition and binding in archaea

To date, methanogens are the only group within the archaea where firing DNA replication origins have not been demonstrated in vivo. In the present study we show that a previously identified cluster of ORB (origin recognition box) sequences do indeed function as an origin of replication in vivo in the archaeon Methanothermobacter thermautotrophicus. Although the consensus sequence of ORBs in M. thermautotrophicus is somewhat conserved when compared with ORB sequences in other archaea, the Cdc6-1 protein from M. thermautotrophicus (termed MthCdc6-1) displays sequence-specific binding that is selective for the MthORB sequence and does not recognize ORBs from other archaeal species. Stabilization of in vitro MthORB DNA binding by MthCdc6-1 requires additional conserved sequences 3' to those originally described for M. thermautotrophicus. By testing synthetic sequences bearing mutations in the MthORB consensus sequence, we show that Cdc6/ORB binding is critically dependent on the presence of an invariant guanine found in all archaeal ORB sequences. Mutation of a universally conserved arginine residue in the recognition helix of the winged helix domain of archaeal Cdc6-1 shows that specific origin sequence recognition is dependent on the interaction of this arginine residue with the invariant guanine. Recognition of a mutated origin sequence can be achieved by mutation of the conserved arginine residue to a lysine or glutamine residue. Thus despite a number of differences in protein and DNA sequences between species, the mechanism of origin recognition and binding appears to be conserved throughout the archaea.

The impact of comparative genomics on infectious disease research

The past decade has witnessed a revolution in infectious disease research, fuelled by the accumulation of a huge amount of DNA sequence data. The avalanche of genome sequence information has largely promoted the development of comparative genomics, which exploits available genome sequences to perform either inter- or intra-species comparisons of bacterial genome contents, or performs comparisons between the human genome and those of other organisms. This review aims to summarize how comparative genomics is being extensively used in infectious disease research, such as in the studies to identify virulence determinants, antimicrobial drug targets, vaccine candidates and new markers for diagnostics. These applications hold considerable promise for alleviating the burden of infectious diseases in the coming years.

Archaeal cell cycle progress

The discovery of multiple chromosome replication origins in Sulfolobus species has added yet another eukaryotic trait to the archaea, and brought new levels of complexity to the cell cycle in terms of initiation of chromosome replication, replication termination and chromosome decatenation. Conserved repeated DNA elements--origin recognition boxes--have been identified in the origins of replication, and shown to bind the Orc1/Cdc6 proteins involved in cell cycle control. The origin recognition boxes aid in the identification and characterization of new origins, and their conservation suggests that most archaea have a similar replication initiation mechanism. Cell-cycle-dependent variation in Orc1/Cdc6 levels has been demonstrated, reminiscent of variations in cyclin levels during the eukaryotic cell cycle. Information about archaeal chromosome segregation is also accumulating, including the identification of a protein that binds to short regularly spaced repeats that might constitute centromere-like elements. In addition, studies of cell-cycle-specific gene expression have potential to reveal, in the near future, missing components in crenarchaeal chromosome replication, genome segregation and cell division. Together with an increased number of physiological and cytological investigations of the overall organization of the cell cycle, rapid progress of the archaeal cell cycle field is evident, and archaea, in particular Sulfolobus species, are emerging as simple and powerful models for the eukaryotic cell cycle.

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.