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

sefOri: selecting the best-engineered sequence features to predict DNA replication origins

Motivation

Cell divisions start from replicating the double-stranded DNA, and the DNA replication process needs to be precisely regulated both spatially and temporally. The DNA is replicated starting from the DNA replication origins. A few successful prediction models were generated based on the assumption that the DNA replication origin regions have sequence level features like physicochemical properties significantly different from the other DNA regions.

Results

This study proposed a feature selection procedure to further refine the classification model of the DNA replication origins. The experimental data demonstrated that as large as 26% improvement in the prediction accuracy may be achieved on the yeast Saccharomyces cerevisiae. Moreover, the prediction accuracies of the DNA replication origins were improved for all the four yeast genomes investigated in this study.

Availability and implementation

The software sefOri version 1.0 was available at http://www.healthinformaticslab.org/supp/resources.php. An online server was also provided for the convenience of the users, and its web link may be found in the above-mentioned web page.

Supplementary information

Supplementary data are available at Bioinformatics online.

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.

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.

GenomeFingerprinter: the genome fingerprint and the universal genome fingerprint analysis for systematic comparative genomics

Background

No attention has been paid on comparing a set of genome sequences crossing genetic components and biological categories with far divergence over large size range. We define it as the systematic comparative genomics and aim to develop the methodology.

Results

First, we create a method, GenomeFingerprinter, to unambiguously produce a set of three-dimensional coordinates from a sequence, followed by one three-dimensional plot and six two-dimensional trajectory projections, to illustrate the genome fingerprint of a given genome sequence. Second, we develop a set of concepts and tools, and thereby establish a method called the universal genome fingerprint analysis (UGFA). Particularly, we define the total genetic component configuration (TGCC) (including chromosome, plasmid, and phage) for describing a strain as a systematic unit, the universal genome fingerprint map (UGFM) of TGCC for differentiating strains as a universal system, and the systematic comparative genomics (SCG) for comparing a set of genomes crossing genetic components and biological categories. Third, we construct a method of quantitative analysis to compare two genomes by using the outcome dataset of genome fingerprint analysis. Specifically, we define the geometric center and its geometric mean for a given genome fingerprint map, followed by the Euclidean distance, the differentiate rate, and the weighted differentiate rate to quantitatively describe the difference between two genomes of comparison. Moreover, we demonstrate the applications through case studies on various genome sequences, giving tremendous insights into the critical issues in microbial genomics and taxonomy.

Conclusions

We have created a method, GenomeFingerprinter, for rapidly computing, geometrically visualizing, intuitively comparing a set of genomes at genome fingerprint level, and hence established a method called the universal genome fingerprint analysis, as well as developed a method of quantitative analysis of the outcome dataset. These have set up the methodology of systematic comparative genomics based on the genome fingerprint analysis.

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.

Mapping of origin of replication in Themococcales

Genome replication is a crucial and essential process for the continuity of life.In all organisms it starts at a specific region of the genome known as origin of replication (Ori) site. The number of Ori sites varies in prokaryotes and eukaryotes. Replication starts at a single Ori site in bacteria, but in eukaryotes multiple Ori sites are used for fast copying across all chromosomes. The situation becomes complex in archaea, where some groups have single and others have multiple origins of replication. Themococcales, are a hyperthermophilic order of archaea. They are anaerobes and heterotrophs-peptide fermenters, sulphate reducers, methanogens being some of the examples of metabolic types. In this paper we have applied a combination of multiple in silico approaches - Z curve, the cell division cycle (cdc6) gene location and location of consensus origin recognition box (ORB) sequences for location of origin of replication in Thermococcus onnurineus, Thermococcus gammatolerans and other Themococcales and compared the results to that of the well-documented case of Pyrococcus abyssi. The motivation behind this study is to find the number of Ori sites based on the data available for members of this order. Results from this in silico analysis show that the Themococcales have a single origin of replication.

Random and site-specific replication termination

Bi-directionality is a common feature observed for genomic replication for all three phylogenetic kingdoms: Eubacteria, Archaea, and Eukaryotes. A consequence of bi-directional replication, where the two replication forks initiated at an origin move away from each other, is that the replication termination will occur at positions away from the origin sequence(s). The replication termination processes are therefore physically and mechanistically dissociated from the replication initiation. The replication machinery is a highly processive complex that in short time copies huge numbers of bases while competing for the DNA substrate with histones, transcription factors, and other DNA-binding proteins. Importantly, the replication machinery generally wins out; meanwhile, when converging forks meet termination occurs, thus preventing over-replication and genetic instability. Very different scenarios for the replication termination processes have been described for the three phylogenetic kingdoms. In eubacterial genomes replication termination is site specific, while in archaea and eukaryotes termination is thought to occur randomly within zones where converging replication forks meet. However, a few site-specific replication barrier elements that mediate replication termination have been described in eukaryotes. This review gives an overview about what is known about replication termination, with a focus on these natural site-specific replication termination sites.

In silico prediction of the origin of replication among bacteria: a case study of Bacteroides thetaiotaomicron

The initiation of chromosomal replication occurs only once during the prokaryote cell cycle. Some origins of replication have been experimentally determined and have led to the development of in silico approaches to find the origin of replication among other prokaryotes. DNA base composition asymmetry is the basis of numerous in silico methods used to detect the origin and terminus of replication in prokaryotes. However, the composition asymmetry does not allow us to locate precisely the positions of the origin and terminus. Since DNA replication is a key step in the cell cycle it is important to determine properly the origin and terminus regions. Therefore, we have reviewed here the methods, tools, and databases for predicting the origins and terminuses of replication, and we have proposed some complementary analyses to reinforce these predictions. These analyses include finding the dnaA gene and its binding sites; making BLAST analyses of the intergenic sequences compared to related species; studying the gene order around the origin sequence; and studying the distribution of the genes encoded in the leading versus the lagging strand.

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.

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.

Separate base usages of genes located on the leading and lagging strands in Chlamydia muridarum revealed by the Z curve method

Background

The nucleotide compositional asymmetry between the leading and lagging strands in bacterial genomes has been the subject of intensive study in the past few years. It is interesting to mention that almost all bacterial genomes exhibit the same kind of base asymmetry. This work aims to investigate the strand biases in Chlamydia muridarum genome and show the potential of the Z curve method for quantitatively differentiating genes on the leading and lagging strands.

Results

The occurrence frequencies of bases of protein-coding genes in C. muridarum genome were analyzed by the Z curve method. It was found that genes located on the two strands of replication have distinct base usages in C. muridarum genome. According to their positions in the 9-D space spanned by the variables u₁ – u₉ of the Z curve method, K-means clustering algorithm can assign about 94% of genes to the correct strands, which is a few percent higher than those correctly classified by K-means based on the RSCU. The base usage and codon usage analyses show that genes on the leading strand have more G than C and more T than A, particularly at the third codon position. For genes on the lagging strand the biases is reverse. The y component of the Z curves for the complete chromosome sequences show that the excess of G over C and T over A are more remarkable in C. muridarum genome than in other bacterial genomes without separating base and/or codon usages. Furthermore, for the genomes of Borrelia burgdorferi, Treponema pallidum, Chlamydia muridarum and Chlamydia trachomatis, in which distinct base and/or codon usages have been observed, closer phylogenetic distance is found compared with other bacterial genomes.

Conclusion

The nature of the strand biases of base composition in C. muridarum is similar to that in most other bacterial genomes. However, the base composition asymmetry between the leading and lagging strands in C. muridarum is more significant than that in other bacteria. It's supposed that the remarkable strand biases of G/C and T/A are responsible for the appearance of separate base or codon usages in C. muridarum. On the other hand, the closer phylogenetic distance among the four bacterial genomes with separate base and/or codon usages is necessary rather than occasional. It's also shown that the Z curve method may be more sensitive than RSCU when being used to quantitatively analyze DNA sequences.

Replication-associated strand asymmetries in vertebrate genomes and implications for replicon size, DNA replication origin, and termination

Strand compositional asymmetry has been observed in prokaryotes and used in predicting prokaryotic DNA replication origins and termini. However, it was not found in eukaryotic genomes by the same methods. We propose that transcription-associated strand asymmetries mask the replication-associated ones. By analyzing the nucleotide composition of intergenic sequences larger than 50 kb by cumulative skew diagrams (CSD), we found replication-associated strand asymmetry in vertebrate genomes. Furthermore, we found that the most common replicon sizes in vertebrates are 50-100 kb, and show evidence that the replication origin and termination regions of vertebrate genomes range from a discrete site to a broad zone.

Origin of replication in circular prokaryotic chromosomes

To predict origins of replication in prokaryotic chromosomes, we analyse the leading and lagging strands of 200 chromosomes for differences in oligomer composition and show that these correlate strongly with taxonomic grouping, lifestyle and molecular details of the replication process. While all bacteria have a preference for Gs over Cs on the leading strand, we discover that the direction of the A/T skew is determined by the polymerase-alpha subunit that replicates the leading strand. The strength of the strand bias varies greatly between both phyla and environments and appears to correlate with growth rate. Finally we observe much greater diversity of skew among archaea than among bacteria. We have developed a program that accurately locates the origins of replication by measuring the differences between leading and lagging strand of all oligonucleotides up to 8 bp in length. The program and results for all publicly available genomes are available from http://www.cbs.dtu.dk/services/GenomeAtlas/suppl/origin.

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, Methanocaldococcus jannaschii

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.

Multiple replication origins of the archaeon Halobacterium species NRC-1

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