Skip to main content
  • Research article
  • Open access
  • Published:

Diversity and evolution of multiple orc/cdc6-adjacent replication origins in haloarchaea

Abstract

Background

While multiple replication origins have been observed in archaea, considerably less is known about their evolutionary processes. Here, we performed a comparative analysis of the predicted (proved in part) orc/cdc6-associated replication origins in 15 completely sequenced haloarchaeal genomes to investigate the diversity and evolution of replication origins in halophilic Archaea.

Results

Multiple orc/cdc6-associated replication origins were predicted in all of the analyzed haloarchaeal genomes following the identification of putative ORBs (origin recognition boxes) that are associated with orc/cdc6 genes. Five of these predicted replication origins in Haloarcula hispanica were experimentally confirmed via autonomous replication activities. Strikingly, several predicted replication origins in H. hispanica and Haloarcula marismortui are located in the distinct regions of their highly homologous chromosomes, suggesting that these replication origins might have been introduced as parts of new genomic content. A comparison of the origin-associated Orc/Cdc6 homologs and the corresponding predicted ORB elements revealed that the replication origins in a given haloarchaeon are quite diverse, while different haloarchaea can share a few conserved origins. Phylogenetic and genomic context analyses suggested that there is an original replication origin (oriC1) that was inherited from the ancestor of archaea, and several other origins were likely evolved and/or translocated within the haloarchaeal species.

Conclusion

This study provides detailed information about the diversity of multiple orc/cdc6-associated replication origins in haloarchaeal genomes, and provides novel insight into the evolution of multiple replication origins in Archaea.

Background

DNA replication is an essential process and is generally conserved across all three domains of life, making use of two different DNA replication apparatuses (bacterial-type and eukaryotic-type)[1, 2]. DNA replication initiates from a single origin in bacteria, whereas multiple origins are utilized in eukaryotes[3]. The study of replication origins in archaea has been ongoing for more than a decade, and multiple replication origins have been identified in several archaeal species[410]. It is not yet understood, however, why archaea adopt multiple origins to initiate replication of their bacterial-like chromosomes or how these multiple replication origins evolved. Notably, answering these questions may also provide insight into the mechanisms of the more complex replication origins found in eukaryotes.

Our current knowledge of archaeal replication origins comes from a few model strains, including Pyrococcus abyssi (Euryarchaeota)[4, 11, 12], Sulfolobus spp. (Crenarchaeota)[5, 6, 13], and two halophilic archaea (Euryarchaeota), Haloferax volcanii[9] and Halobacterium sp. strain NRC-1[10, 14]. A single origin was mapped near the only orc1 gene in the genome of P. abyssi via a combination of in silico[4] and biochemical[11, 12] analyses. Three origins were mapped in Sulfolobus solfataricus and Sulfolobus acidocaldarius by means of microarray-based marker frequency analysis (MFA)[6], and two origins adjacent to the orc/cdc6 genes were previously identified by two-dimensional gel electrophoresis[5]. On the chromosome of Halobacterium sp. NRC-1, one origin was verified with autonomous replication activity[14], and four were mapped using whole-genome MFA, with three origins located in the vicinity of orc/cdc6 genes[10]. In H. volcanii, five replication origins were identified, two within the chromosome and one each within the three megaplasmids pHV1, pHV3 and pHV4[9].

These experimental data revealed that the basic structure of replication origins is conserved among archaea, normally containing an AT-rich unwinding element and several conserved repeats (Origin Recognition Box, ORB)[9]. The ORB elements were proven to be the recognition sites for the Orc/Cdc6 initiation protein via biochemical[5] and structural approaches[15, 16]. In addition, distinct from the ORBs identified in the oriC1 of S. solfataricus[5], a halophile-specific “G-string” (long G-stretches locating at the end of ORBs) was observed in all origins from H. volcanii[9]. Whereas the Cdc6 and the ORC complex proteins (Orc1-6) act together to recruit the MCM (minichromosome maintenance) complex to an origin of replication in eukaryotes[3, 17], a subset of initiator proteins (Orc/Cdc6), which are related to both Orc1 and Cdc6 of eukaryotes, were adopted by archaea. Therefore, archaeal Orc/Cdc6 is considered to possess both origin recognition and MCM-loading activities[3]. Previous studies in S. solfataricus revealed that origin identity was determined by the specific recognition of Orc/Cdc6 proteins[18]. Interestingly, the multiple origins, especially the ORB sequences and their associated Orc/Cdc6 proteins, are quite diverse in all three experimentally characterized archaea (S. solfataricus, H. volcanii and Halobacterium sp. NRC-1)[5, 9, 10], indicating independent evolutionary history. In particular, an origin comparison between two hyperthermophilic archaeal genera, Aeropyrum and Sulfolobus, suggests that the capture of extrachromosomal elements accounts for replicon evolution[7]. However, as the study of replication origins has been limited to only a few archaeal species, it is still difficult to determine the evolutionary relationship of multiple replication origins within Archaea.

Haloarchaeal genomes are normally composed of multiple replicons (chromosome, minichromosome, and plasmids) with multiple Orc/Cdc6 homologs (usually more than 10 homologs)[1930], indicating that the occurrence of multiple replication origins is widespread in haloarchaea. To date, however, studies have been limited to two model haloarchaea, H. volcanii[9] and Halobacterium sp. NRC-1[10]. There is not enough information to understand the diversity and evolution of multiple replication origins in this distinct group of archaea. In this work, which is based on both previous experimental data and the identification of replication origins in Haloarcula hispanica, a haloarchaeon that was recently sequenced in our laboratory[20], we performed a comparative analysis of predicted orc/cdc6-associated replication origins in 15 completely sequenced haloarchaeal genomes. These comparative analyses indicated that the introduction of novel replication origins usually accompanied the acquisition of new genomic content by insertions into a chromosome or the reconstruction of novel extrachromosomal replicons, which may be linked to an adaptive mechanism of haloarchaea to diverse environments, similar to Salinibacter ruber[31]. In particular, various families of orc/cdc6-associated replication origins were identified, and different evolutionary mechanisms, including ancestral preservation, translocation among haloarchaea and likely differential loss, were proposed to account for the current multiple origins of replication in the haloarchaeal genomes.

Results and discussion

Identification of orc/cdc6-associated replication origins in H. hispanica

Research in archaeal model strains indicates that most replication origins share conserved characteristics, such as the presence of inverted ORB elements and being located directly adjacent to orc/cdc6 genes[5, 9, 10]. In addition, a “G-string” at the end of ORB elements was observed in all identified origins from H. volcanii[9]. These common features provided us a reference standard to predict replication origins in H. hispanica. Briefly, only those intergenic regions (IRs) that contain ORB-like elements and are directly adjacent to orc/cdc6 genes were considered to be putative orc/cdc6-associated replication origins. Necessarily, although they were not included in the scope of this study, we do not exclude the possibility of replication origins that are not directly adjacent to orc/cdc6 genes or are without classical ORB-like elements. Replication origins with these characteristics were shown to exist in Sulfolobus spp.[5] and may exist in Halobacterium sp. NRC-1[10]; however, to our knowledge, they constitute only a small proportion of the replication origins in archaea and are not easily predicted with current information.

To identify replication origins in H. hispanica, the IRs around the orc/cdc6 genes were examined for the presence of ORB elements. H. hispanica encodes eleven orc/cdc6 genes, with six copies (cdc6A F) in the main chromosome, four (cdc6G J) in the minichromosome and only one (cdc6K) in the megaplasmid. ORB repeats harboring a G-rich motif were observed adjacent to eight orc/cdc6 genes (Additional file1 and Figure 1A), in agreement with the halophile-specific “G-string” elements found in H. volcanii[9]. However, in contrast to other characterized archaeal origins with at least two ORB repeats flanking an AT-rich unwinding element, only one ORB-like element was observed in each IR flanking the cdc6D gene, which was considered to be a deficient origin (oriC3-cdc6D*) when examined by hand (Figure 1A and Additional file1). Accordingly, seven replication origins were predicted in H. hispanica: two were in the main chromosome (oriC1 cdc6A and oriC2 cdc6E), four were in the minichromosome (oriC4 cdc6G, oriC5 cdc6H, oriC6 cdc6I and oriC7 cdc6J), and one was in the megaplasmid (oriP cdc6K) (Figure 1A and B).

Figure 1
figure 1

Bioinformatic and genetic identification of replication origins in H. hispanica. A. Seven replication origins, oriC1-cdc6A and oriC2-cdc6E in the main chromosome; oriC4-cdc6G, oriC5-cdc6H, oriC6-cdc6I and oriC7-cdc6J in the minichromosome; and oriP-cdc6K in the megaplasmid, were predicted by searching ORB motifs (indicated with small triangles) in the IRs located directly adjacent to orc/cdc6 genes (indicated with red arrowheads) using MEME software. Logo representations of ORB elements are presented on the right, and the spaces represent sequences that are not conserved. oriC3*: predicted deficient origin adjacent to cdc6D gene. B. Replication assay for plasmids containing the origins predicted in A. (Up) Southern blot analysis with a bla gene probe: lane T contains crude DNA extracted from the H. hispanica transformants, and lane P represents the purified plasmid as an input control; (down) summaries of the identification of origins in H. hispanica and the five origins with ARS activity (oriC 1, 2, 6, 7, P) are indicated with filled ovals and are bolded in A.

To confirm these putative replication origins, we performed a genetic assay to test their autonomous replication activities. As a control, we also examined whether oriC3-cdc6D* and IRs around cdc6B, cdc6C and cdc6F, where no ORBs were detected, could engage in origin activities. DNA fragments, including the orc/cdc6 genes plus their flanking IRs, were cloned into a nonreplicating plasmid, pBI101[32, 33], to assay for the presence of an autonomously replicating sequence (ARS) (Figure 1, Additional file2). Of the eleven orc/cdc6 genes with adjacent IRs, oriC1 cdc6A and oriC2 cdc6E in the main chromosome, oriC6 cdc6I and oriC7 cdc6J in the minichromosome and oriP-cdc6K in the megaplasmid were able to confer replication ability to the non-replicating plasmid (Figure 1B, Additional file2), which was indicative of the ARS activities of these origins. As expected, no replicating ability was observed for plasmids constructed with oriC3-cdc6D* or with the fragments containing cdc6B, cdc6C and cdc6F (Additional file2). Although the remaining two predicted replication origins, oriC4-cdc6G and oriC5-cdc6H, shared a conserved structure with characteristic archaeal origin (Additional file1), they could not drive the autonomous replication activities under our experimental conditions, which is reminiscent of the dormant origins found in eukaryotes[34]. Dormant replication origins are normally inactive, but they can be activated for cellular response to replicative stress[35, 36]. In the future, it would be interesting to further analyze the utilization of these likely dormant replication origins in H. hispanica.

Most orc/cdc6 genes are predicted to associate with replication origins in haloarchaea

To date, the genomes of 15 haloarchaea have been made available through NCBI (before October 2011), and 14 of these 15 genomes include the minichromosomes and/or megaplasmids, which provided us the opportunity to perform a comparative genomic analysis of replication origins in haloarchaea. To focus on the orc/cdc6-associated replication origins, we first conducted an exhaustive search of the orc/cdc6 genes in the 15 sequenced haloarchaeal genomes (Table 1).

Table 1 Predicted origin-associated Orc/Cdc6 homologs in the haloarchaeal genomes

Multiple Orc/Cdc6 homologs are encoded in each of the 15 sequenced haloarchaeal genomes. Based on a previous study[15], origin-associated Orc/Cdc6 proteins contain two important domains, a N-terminal AAA + domain and a C-terminal winged-helix domain, and almost all have a length greater than 300 amino acids. A total of 154 Orc/Cdc6 homologs fulfilling these criteria were collected from the 15 sequenced haloarchaeal genomes (Table 1 and Additional file3), and the IRs flanking these orc/cdc6 genes were collected for motif searching. Interestingly, distinct ORB-like elements harboring G-string were found in the IRs flanking nearly two-thirds (102 of 154) of the orc/cdc6 genes (Table 1 and Additional files3 and4), and the predicted replication origins were rechecked manually to remove deficient origins such as oriC3-cdc6D* in H. hispanica. As expected, multiple replication origins were predicted in all of the analyzed haloarchaeal genomes (Table 1). Haloterrigena turkmenica has the greatest number of predicted origins at 12, and 7 of those origins are located on its chromosome (Table 1). On average, within the haloarchaeal chromosomes, more than half of the orc/cdc6 genes have predicted origins nearby: a maximum of 75% (3 of 4) in Halobacterium spp. and a minimum of 33% (1 of 3) in Halomicrobium mukohataei (Table 1). Compared with the chromosome, the overwhelming majority (greater than 80%) of the orc/cdc6 genes in the extrachromosomal elements (minichromosomes and megaplasmids) are associated with predicted replication origins (Table 1).

As several replication origins have been experimentally mapped in H. hispanica (Figure 1), H. volcanii[9] and Halobacterium sp. NRC-1[10], these experimental data were used to evaluate the efficiency of the origin prediction performed in this study. In H. hispanica, five out of the seven predicted replication origins were confirmed to have ARS activity. For the replication origins in the chromosome of Halobacterium sp. NRC-1, a high consistency between our predicted results and the whole-genome MFA[10], except for one uncertain origin (oriC4), proves the efficiency of the bioinformatic approach in this study. For H. volcanii, in addition to the five replication origins that were previously genetically mapped[9], five additional replication origins were also predicted in this study. As discussed above, these additional predicted origins might be weak or dormant replication origins, which are not easily identified by experimental approaches.

In summary, our bioinformatic approach not only is important for identifying active replication origins in haloarchaea but also provides novel information for predicting likely dormant replication origins, which is also important for the future study of replication regulation and adaptation in archaea.

Diversity of orc/cdc6-associated replication origins in haloarchaea

A recent report suggested that Orc/Cdc6 initiators specifically determine origin discrimination in archaea[18]. To investigate this further, a phylogenetic analysis of ori-associated Orc/Cdc6 proteins in haloarchaea was performed, and the resulting tree showed that Orc/Cdc6 homologs cluster into different families (Figure 2A), which suggested that various orc/cdc6-associated replication origins have been adopted in haloarchaea. Different Orc/Cdc6 families have been suggested in previous work[14, 37]; herein, we focused on the putative origin-associated Orc/Cdc6 homologs with the intention of providing a detailed classification of predicted replication origins. Although setting precise boundaries was difficult, the predicted replication origins could be sorted into distinct families based on a combination of the phylogenetic tree of the Orc/Cdc6 homologs (Figure 2A) and a comparison of ORB sequences (Figure 2B). It is noteworthy that BLAST analyses confirmed that only those Orc/Cdc6 homologs showing high identities (at least 80%) were grouped into the same family in this study. Specifically, the origins adjacent to the specific Orc/Cdc6 conserved among all haloarchaea were named oriC1, as previously reported[10, 37], and two other families with the top two members were selected and named oriCa and oriCb to facilitate additional evolutionary analyses (Figure 2A). This classification of replication origins will become more complete when more haloarchaeal genomes become available and will aid in the understanding of replication origins in novel haloarchaea.

Figure 2
figure 2

Families of orc/cdc6 -associated replication origins in the haloarchaeal genomes. A. Phylogenetic tree of origin-associated Orc/Cdc6 homologs (Hbo: Halogeometricum borinquense, Hhi: Haloarcula hispanica, Hje: Halalkalicoccus jeotgali B3, Hla: Halorubrum lacusprofundi, Hma: Haloarcula marismortui, Hmu: Halomicrobium mukohataei, HR1: Halobacterium salinarum R1, Htu: Haloterrigena turkmenica, Hut: Halorhabdus utahensis, Hvo: Haloferax volcanii DS2, Hwa: Haloquadratum walsbyi, Nma: Natrialba magadii, Nph: Natronomonas pharaonis, NRC-1: Halobacterium sp. NRC-1). * indicates the Orc/Cdc6 proteins encoded on extrachromosomal elements. Orc/Cdc6 proteins from other archaea (APE: Aeropyrum pernix, Pab: Pyrococcus abyssi, Sso: Sulfolobus solfataricus) are highlighted with yellow background. The origin conserved in all genomes was assigned the name oriC1, as in previous reports (in red), and the other two origin clusters with the top two members, excluding oriC1 in this study, were assigned the names oriCa (in blue) and oriCb (in pink). B. Logo representations of the putative ORB elements identified in the predicted replication origins adjacent to the orc/cdc6 genes. The ORB elements are not shown for origin families with only one member.

Notably, Orc/Cdc6 proteins associated with oriC1 not only are conserved in all haloarchaeal genomes but were also identified in other archaeal genomes (Figure 2A). For instance, Cdc6-1 in S. solfataricus was experimentally proven to recognize the ORB elements of oriC1 in Halobacterium sp. NRC-1[5]. Thus, we suggest that this conserved origin (oriC1) might be present in an archaea ancestor and maintained in different lineages during the evolutionary history of Archaea. By contrast, other Orc/Cdc6 homologs from different haloarchaeal genomes could be clustered into several distinct families (Figure 2A), indicating that these origins are shared by different haloarchaea. As different Orc/Cdc6 proteins from the same haloarchaeal genome are normally distributed into different families (Figure 2A), the results suggested that multiple origins within a haloarchaeon were unlikely to arise from internal duplications.

To further characterize these predicted origins in haloarchaea, ORB sequences were extracted from all of the putative origins (Additional file4) for comparison analysis (Figure 2B). The results revealed the linkage-specificity of Orc/Cdc6 homologs and the corresponding ORB elements, (i.e., the predicted ORB sequences proximal to Orc/Cdc6 homologs within the same family are highly similar) (Figure 2B and Additional file5). These observations suggested that Orc/Cdc6 proteins specifically recognize adjacent ORB elements, consistent with a recent report about origin discrimination by Orc/Cdc6 initiators[18].

Notably, the structures of origins from even the same family are also diverse in haloarchaea. Although the predicted replication origins are primarily located directly upstream of orc/cdc6 genes, as previously observed[4, 5, 9, 10], there are a few exceptions (Additional file4) in which ORB elements are located in IRs downstream of the orc/cdc6 gene (e.g., proximal to cdc6n in H. marismortui, orc5 in Halobacterium species, orc8 in H. volcanii and Nmag_3611 in N. magadii) (Additional file4). In addition, in some predicted replication origins, ORB elements are occasionally observed in both IRs flanking the orc/cdc6 gene, such as those flanking cdc6E in H. hispanica, cdc6h, cdc6g and cdc6m in H. marismortui, orc10 and orc8 in Halobacterium species, orc4 in H. volcanii and Hlac_1078, Hlac_2747 and Hlac_2997 in H. lacusprofundi (Additional file4). Interestingly, although a previous genetic experiment found that the predicted origin proximal to NRC-1_orc8 is not able to promote efficient autonomous replication[14], this origin, in combination with the origin proximal to NRC-1_orc10, was experimentally proven to be active in Halobacterium sp. strain NRC-1 by whole-genome MFA[10]. The origin proximal to orc10 in Halobacterium species has a different structure comparing to others in the oriCa family (Additional file4), suggesting that structurally diverse origins function in different haloarchaea.

Surprisingly, two different ORB-like elements were found in the IRs proximal to one orc/cdc6 gene, in either the megaplasmid of H. lacusprofundi (Hla_3512) or the main chromosome of H. walsbyi (Hwa_HQ2959A) (Figure 3A). One ORB-like element, TAACAGCGGAAACAGTGGGGTGGGGGGGT, is shared by these two different origins, while the other shows no similarity (Figure 3B).

Figure 3
figure 3

Two different candidate replication origins are adjacent to one orc/cdc6 gene (Hla_3512 or Hwa_HQ2959A). A. The sequence features of the two putative origins. The orc/cdc6 genes are indicated with black-boxed arrows, and their start site is numbered 1. The adjacent origins (I and II) were mapped with small arrowheads and rectangles indicating the ORB elements and AT-rich regions, respectively. Origin II, in both cases, contains conserved ORB elements, which are highlighted in blue. B. Logo representations of the ORB elements in the four candidate replication origins. The ORB elements (boxed) are highly conserved in Origin II in both cases.

Taken together, the multiple replication origins in haloarchaea are dramatically diverse. In addition to the diversity of the ORB elements and corresponding orc/cdc6 genes, the number and position of ORB elements also contributes to the diversity of the origins in haloarchaea. This diversity may facilitate the differential utilization of multiple replication origins in haloarchaea. Strikingly, the origin proximal to orc10 in Halobacterium species was active in vivo[10], while its conserved origin (oriCa family) proximal to orc2 in H. volcanii was not proven functional[9]. As the two origins exhibit different structures and these two haloarchaea grow in different environmental conditions, these observations may provide novel insight into differential utilization of replication origins in haloarchaea.

Novel replication origins accompany newly acquired genomic content

As described above, the replication origins of two Haloarcula species, H. hispanica and H. marismortui were predicted, and their ARS activities were also examined in H. hispanica (Figure 1). Although their chromosomes show a high degree of conservation (Figure 4B), the two species harbor several different replication origins (Table 1 and Figure 4A). Thus, an in-depth study of these origins would be helpful in understanding the processes involved in the diversity of haloarchaeal replication origins.

Figure 4
figure 4

Comparative analysis of the orc/cdc6 -associated replication origins between the chromosomes of H. hispanica and H. marismortui. A. Distribution of the candidate orc/cdc6-associated replication origins in the chromosomes of H. hispanica (inside) and H. marismortui (outside). G + C content of the chromosome of H. hispanica was plotted, and significant variations in the two divergent regions are indicated with blue arrows. The predicted orc/cdc6-associated replication origins are indicated as ovals on the chromosome circle, and the shared orc/cdc6-associated replication origins in the two Haloarcula species, oriC1 and oriC2, are highlighted as filled ovals. B. Genome alignment of the chromosomes of H. hispanica and H. marismortui. Their shared orc/cdc6-associated replication origins are indicated as in A. Regions A and B represent discrepancies between the two chromosomes, which are exactly in accordance with the positions of their specific orc/cdc6-associated replication origins; oriC3-cdc6D* of H. hispanica and oriC3-cdc6i of H. marismortui are located in region A, and oriC4-cdc6g of H. marismortui is located in region B. The divergent regions and the edges of the similar regions were confirmed by BLASTN alignments of sequences, and shaded regions denote a similarity of over 70%. Linearized scaled bars are provided. C. A schematic representation of the two divergent regions (1 kb scale for Hhis_A, Hmar_A and Hhis_B; 2 kb scale for Hmar_B) between the two chromosomes. The orc/cdc6 genes are indicated. The polysaccharide biosynthesis genes are in yellow, transposase genes in purple, other genes with known functions in pink and hypothetical genes in gray. The species with the closest matches in the BLAST analysis is indicated on top of the gene: M, Methanobacterium; A, other non-halophilic archaea; B, eubacteria (the colors are designed to correspond to the marks in Additional file6). The genes in clusters are also in clusters in other haloarchaea, as indicated at the top of the clusters.

With the exception of the two shared replication origins, oriC1 cdc6A and oriC2-cdc6E in H. hispanica and the corresponding oriC1-cdc6d and oriC2-cdc6h in H. marismortui, there are one or two other predicted origins specific to each strain: oriC3-cdc6D* in H. hispanica, oriC3-cdc6i and oriC4-cdc6g in H. marismortui (Figure 4A and B). The two shared origins, oriC1 and oriC2, were likely present in the ancestor of Haloarcula, and their specific origins, oriC3 cdc6D* in H. hispanica and oriC3 cdc6i and oriC4 cdc6g in H. marismortui, may have been acquired later through translocation processes following the divergence of these species. An alternative hypothesis is that all three species-specific origins were also present in the ancestor of Haloarcula but were lost differently in H. hispanica and H. marismortui. However, these three predicted origins (oriC3 cdc6D*, oriC3 cdc6i and oriC4 cdc6g) are located in two divergent regions (region A and B in Figure 4A and B) with significant G + C content variations (Figure 4A), which is indicative of newly acquired genomic content specific to each of the two strains[38]. Thus, the most likely explanation is that these predicted species-specific origins were newly acquired as a part of new genomic content (i.e., the haloarchaeal genomes might recruit novel replication origins accompanying new genes). This hypothesis is reinforced by the abundance of transposases observed around these specific origins (Figure 4C and Additional file3).

Concentrating on the genes with annotated functions, except for those predicted to be transposases, the majority of genes within the two divergent regions were found to be involved in polysaccharide biosynthesis (Figure 4C). Subsequently, a BLAST analysis against the NCBI non-redundant proteins database was performed on all of the genes in regions A and B in both chromosomes (Figures 4B and C and Additional file6). The genes were conserved across several different organisms (Figure 4C); most were similar to other haloarchaeal homologs, but for several genes, their closest homologs were outside of haloarchaea. The two linked glycosyltransferase genes in region A of H. hispanica were most similar to those found in Methanobacterium (Figure 4C and Additional file6). Several genes in region B of both chromosomes showed the greatest similarity to genes found in bacteria, especially a cluster in region B of H. hispanica (Figure 4C and Additional file6). In addition, those genes found in clusters in the two Haloarcula species were also usually found in clusters in other organisms (Figure 4C), suggesting that these genes were acquired in clusters.

A previous report in Salinibacter ruber suggested that genes with related functions but different origins might have been assembled together and introduced concurrently into the genome of S. rubber[31]. Similarly, our comparative analyses indicated that the convergence of closely related functional genes from different sources is an important way through which new genomic content is acquired in haloarchaea and that foreign replication origins are usually introduced as a component of this new content. We cannot be certain whether the new genomic content (mixture of new genes and foreign replication origins) is introduced with single or multiple transfer(s), as the mechanism is not well understood; however, our analyses strongly suggested that the novel replication origins may be important for the acquisition of new genomic content and that the newly acquired genes from the surroundings may be favorable for the haloarchaeal cells to improve their ability to adapt to changeable environments.

Recruitment of novel replication origins in the reconstruction of the extrachromosomal replicons

The haloarchaeal genomes in this study, except that of Halorhabdus utahensis, generally harbor extrachromosomal replicon(s), ranging in number from one in H. mukohataei and H. walsbyi to eight in H. marismortui (Table 1). In addition, orc/cdc6 genes were found on most of the extrachromosomal elements (Table 1), suggesting that the orc/cdc6-associated replication origins are responsible for replication initiation on most of these replicons. Therefore, an in-depth analysis could further elucidate the evolution of these replication origins.

Compared to H. hispanica, the H. marismortui genome contains a greater number of extrachromosomal replicons, with eight (minichromosome II and 7 megaplasmids, pNG100 to pNG700), while H. hispanica contains only two (minichromosome II and megaplasmid pHH400). Among these minireplicons, only megaplasmids pHH400 and pNG700 are collinear (Figure 5), suggesting that they may have been present in a common ancestor of the two Haloarcula species. The lengths of the minichromosomes of H. marismortui and H. hispanica are 288 kb and 488 kb, respectively. They share homology over approximately 100 kb, with a few inversions and gaps (Figure 5), indicating that this region was likely rearranged in the two Haloarcula species and thus that the two minichromosomes are only distantly related. In addition, the megaplasmids from pNG100 to pNG600 are unique to H. marismortui. However, pairs of orthologous to the minichromosome of H. hispanica are observed, especially in pNG500, with orthologs as large as 30 kb (Figure 5). Together with the abundant ISH (i nsertion s equence from H alobacteriaceae) elements encoded in these replicons, our data imply that the extrachromosomal replicons were significantly rearranged after the divergence of the two species and that new DNA contents were acquired from surrounding organisms. These results are also reminiscent of previous reports on the evolution of the large dynamic replicons found in Halobacterium spp.[22, 39].

Figure 5
figure 5

Comparative genomic analysis of the extrachromosomal replicons of H. hispanica and H. marismortui. The orc/cdc6 genes (those from H. hispanica and H. marismortui are highlighted with a purple asterisk and a dark green round dot, respectively) that are associated with candidate replication origins are indicated, and the shared origins associated with cdc6G/cdc6a, cdc6K/cdc6k of the two strains are highlighted in bold. The homologous regions are boxed, and the lines in the box represent the regions that are continuous in H. marismortui.

To understand the different composition of the extrachromosomal elements in the two Haloarcula species, the orc/cdc6-associated replication origins in these minireplicons were also examined. In H. hispanica, four predicted orc/cdc6-associated replication origins are distributed in the minichromosome, and one is present in the megaplasmid pHH400. The two origins (oriC6 cdc6I and oriC7 cdc6J) in the minichromosome and the one (oriP-cdc6K) in pHH400 were confirmed by ARS activity (Figure 1 and5). In H. marismortui, the predicted orc/cdc6-associated replication origins are distributed among the extrachromosomal replicons as follows: two in the minichromosome, one in pNG700, one in pNG600, two in pNG500 and one in pNG100 (Figure 5). No orc/cdc6 genes are encoded by either pNG400 or pNG200, and no candidate replication origin was identified adjacent to the orc/cdc6 gene in pNG300, indicating that other types of replication origins are involved in the initiation of replication in these replicons. This concept is reinforced by the identification of rep genes in these replicons (Table 1)[40]. Among these replication origins, only two are shared by the two Haloarcula species, oriP-cdc6K in pHH400 and the origin (proximal to cdc6k) in pNG700, as well as the origin proximal to cdc6G and cdc6a in the minichromosomes of H. hispanica and H. marismortui, respectively (Figure 5). In contrast to the high conservation found in the megaplasmids pHH400 and pNG700, the region around cdc6G and cdc6a shows no collinearity (Figure 5), strongly suggesting that this origin might not have been present in their ancestor and instead was employed by H. hispanica and H. marismortui after their divergence. Surprisingly, a specific origin (oriC7-cdc6J) in the minichromosome of H. hispanica, which proved functional (Figure 1), was located in the region with high orthology to H. marismortui (Figure 5). This observation suggested that this replication origin was recruited into this region in H. hispanica or was lost in H. marismortui during rearrangement of minichromosomes in the two Haloarcula species. Similarly, the specific origins in pNG600, pNG500 and pNG100 and the rep-associated origins in pNG400, pNG300 and pNG200 were all likely recruited to accomplish the construction of these replicons in H. marismortui.

Multiple evolutionary mechanisms account for multiple orc/cdc6-associated origins in haloarchaea

Our above analysis clearly indicated that the replication origins in haloarchaea are quite diverse and that different haloarchaea can share a few different origins. Although we cannot exclude the possibility that origin loss contributes to mosaic replication origins in haloarchaea, it is unlikely that all of the origins currently shared by different haloarchaea were present in the ancestor of each genus of Halobacteriaceae as oriC1. Archaea species often harbor mobile elements within their genome, which are mobilized via integrases[41] or transposases encoded by insertion sequence (IS) elements[42]. Our comparative analyses of the genomic context of the replication origins in the two Haloarcula species demonstrated the presence of mobile elements near their specific origins (Figure 4). These indicators of translocation processes were further analyzed in the genomes proximal to the origins in other haloarchaea. Forty-two of 102 potential replication origins have integrases or transposases nearby (Table 2 and Additional file3), which might contribute to accelerate the translocation of these origins. In haloarchaeal chromosomes, the ratios of later-acquired origins are comparatively low, with a maximum of 50% for H. marismortui, H. utahensis and H. walsbyi and none in H. borinquense, H. jeotgali B3, H. mukohataei, H. xanaduensis and N. pharaonis (Table 2). By comparison, these later-acquired replication origins are widespread in extrachromosomal elements. For example, they account for 80% (4 of 5), 83% (5 of 6) and 87.5% (7 of 8) of the replication origins in the extrachromosomal elements of H. salinarum R1, H. volcanii DS2 and H. lacusprofundi, respectively (Table 2). These observations suggest that a portion of the replication origins in haloarchaea, especially those in extrachromosomal elements, were introduced through recent translocation processes.

Table 2 Predicted later-acquired replication origins in the haloarchaeal genomes

Contrary to the complete conservation of the replication origin oriC1[10, 37], the other origins are distributed almost randomly among haloarchaea, within both the chromosome and extrachromosomal elements (Figure 6A). To better understand the evolutionary history of these replication origins in haloarchaea, two distinct origin families (oriCa and oriCb, Figure 2), with the top two members excluding oriC1 in this study, were selected for further comparative analyses. Interestingly, while the genes around oriC1 are highly syntenic[10, 37], the genomic context around oriCa reveals no similarity among the different haloarchaeal genomes (Figure 6B and Additional file7). These observations indicated that the origins belonging to this family had different evolutionary processes from those in the oriC1 family. Furthermore, transposases were observed near the origin in five out of eight genomes Halobacterium species (HR1_orc10 and NRC-1_orc10), H. marismortui (Hma_cdc6o), H. jeotgali (Hje_17938), and H. lacusprofundi (Hla_2958)] (Figure 6B). These results suggested that these replication origins were likely mobilized via transposases, implying that association with transposases might result in an acceleration of translocation rates of oriCa among haloarchaea. Additionally, this acceleration may account for the random distribution of this origin family among different haloarchaea.

Figure 6
figure 6

Genome context analyses of the origins in the oriCa and oriCb families. A. Phylogenetic tree based on the 16S rRNA genes (Methanocaldococcus jannaschii was added as an outgroup) and the distribution of the oriCa and oriCb families in the haloarchaeal genomes. + indicates the presence of this family of origins in the chromosome, and (+) indicates its presence in an extrachromosomal element. B and C. The regions around oriCa (B) and oriCb (C). Shaded regions denote similarity greater than 70% by BLASTN analyses. The orc/cdc6 genes are highlighted in red; genes annotated with transposase (tnp) are highlighted in purple; rRNAs, including 16S rRNA, 23S rRNA and 5S rRNA, are highlighted in yellow; tRNAs are highlighted in green, with T, C and A representing the Thr-tRNA, Cys-tRNA and Ala-tRNA genes, respectively. oriCb origins of replication in H. volcanii, H. borinquense and H. lacusprofundi are indicated by teal rectangles.

When focusing on the origins of the oriCb family, two different types were observed. One type consists of origins proximal to Hvo_orc5 in H. volcanii, Hbo_14700 in H. borinquense and Hla_1524 in H. lacusprofundi. Some homologs were detected adjacent to this type of origin in all of the three genomes, showing that this type of origin directly links to a syntenic rRNA region containing the 16S rRNA, Ala-tRNA, 23S rRNA and 5S rRNA genes (Figure 6C and Additional file7). In conjunction with the clustering in the 16S rRNA tree, this finding strongly suggests that these three haloarchaea shared their closest ancestor and that oriCb might be present in their ancestor. Notably, the two functional replication origins in the chromosome of H. volcanii were observed close to the two rRNA operons[9]. This closeness might benefit the preservation of origins over evolutionary time. Another type of origin in the oriCb family, including the origins proximal to Htu_5222 in H. turkmenica, Hje_08365 in H. jeotgali, Nma_3611 in N. magadii and Hxa_0635 in H. xanaduensis, was observed; this origin showed no similarity with respect to the order of the genes flanking the origin in different genomes (Figure 6C and Additional file7). This finding implies a totally distinct evolutionary process. As three out of the four members of this type of origin were identified in extrachromosomal elements, it is plausible to propose that these origins were recruited for the construction of novel extrachromosomal replicons independently. Gene pools in environments were proposed to explain the adaption of prokaryotes under changeable environments[31, 43]. Similarly, the diversity of replication origins can be thought of as a pool of origins in environments that can be recruited for the construction of novel replicons. This hypothesis sheds light on not only the random distribution of conserved origins in different haloarchaea but also the presence of extremely variable extrachromosomal replicons in haloarchaea.

Conclusion

In this study, orc/cdc6-associated replication origins were predicted in 15 sequenced haloarchaeal genomes through Orc/Cdc6 protein analyses and adjacent ORB searching. Multiple replication origins were found in all of the analyzed genomes, and nearly two-thirds of the orc/cdc6 genes were found to be associated with the predicted replication origins. We also experimentally investigated the predicted replication origins in H. hispanica and demonstrated that 5 out of 7 predicted origins possess ARS activity and that the remaining 2 putative replication origins appear to be dormant in experimental conditions. In conjunction with ORB comparisons and phylogenetic analysis of the Orc/Cdc6 homologs, various families of these predicted replication origins were revealed in haloarchaea. The diversity of multiple replication origins in haloarchaea was mainly driven by the diversity of Orc/Cdc6 proteins that specifically associate with distinct ORB elements. Interestingly, origins within the same family may have different functions among the various haloarchaea, e.g., although belonging to the oriCa family, the active origin in Halobacterium sp. NRC-1 (proximal to orc10)[10] was not proven functional in H. volcanii (proximal to orc2)[9]. These observations suggested differential origin utilization under different replicative conditions and demonstrated the advantage of our bioinformatic approaches in the identification of dormant or weak replication origins in haloarchaea.

Phylogenetic analysis of Orc/Cdc6 proteins suggested that multiple replication origins in haloarchaeal genomes can be categorized into at least two types: oriC1, which is present in an ancestor of archaea, and the other origins, which are likely specific in haloarchaea. We also revealed that transposases or integrases flank more than 40% of predicted replication origins; this flanking is indicative of the translocation of a portion of the replication origins among haloarchaea. In conjunction with comparative analyses of two families of replication origins (oriCa and oriCb), we suggested that different evolutionary mechanisms account for the diversity of replication origins in haloarchaea: preservation from ancestors (e.g., oriC1 was maintained from the original ancestor of archaea, and one type of origin in oriCb was maintained from the closest ancestor of H. volcanii, H. borinquense and H. lacusprofundi), differential loss, and translocation among haloarchaea. In particular, a comparative genomic analysis of two Haloarcula species revealed that species-specific origins in the main chromosome were introduced along with new genes, whereas in the extrachromosomal replicons, the recruitment of novel replication origins usually accompanied the construction and/or rearrangement of minireplicons. The concept of an “origins pool” was proposed, and the introduction of novel origins in conjunction with the acquisition of new genomic content may be linked to the mechanisms involved in the adaptation of haloarchaeal cells to changeable environments. Taken together, our analyses of the diversity and evolution of the potential replication origins in haloarchaea may open avenues to understanding the significance of the multiple replication origins in the domain of Archaea.

Methods

Strains, plasmids and culturing

Escherichia coli were grown in Luria-Bertani medium at 37 °C, and 100 μg/mL of ampicillin was added when required. H. hispanica was cultivated at 37 °C in nutrient-rich medium AS-168 (per liter: 5.0 g Bacto Casamino Acids, 5.0 g yeast extract, 1.0 g sodium glutamate, 3.0 g trisodium citrate, 200 g NaCl, 20 g MgSO4 · 7H2O, 2.0 g KCl, traces of FeSO4 · 4H2O and MnCl2 · 4H2O, pH 7.2), and 3 μg/ml of mevinolin was added when required[44]. Plasmid pBI101[32, 33] was used for the investigation of the autonomous replication ability of the predicted origins. These plasmids were normally constructed in E. coli and then introduced into H. hispanica by a polyethylene glycol-mediated transformation method[45, 46].

Autonomous replication ability assay

Each PCR fragment, including the intergenic sequences and orc/cdc6, was amplified (see in Additional file4 for primers) from wild type H. hispanica genomic DNA and was cloned into the nonreplicating plasmid pBI101[32, 33]. After sequencing, the plasmids were then introduced into H. hispanica or the corresponding origin-deletion strains (unpublished data) using a polyethylene glycol-mediated transformation method[45, 46], and the mevinolin-resistant transformants were selected on AS-168 plates with 3 μg/mL of mevinolin. Plasmid recovery in H. hispanica transformants indicated the autonomous replication ability of the corresponding origins, which was verified by Southern blot analysis[9, 14]. Briefly, the transformant on the plate was transferred into 200 μL of double-distilled H2O and 100 μL of phenol-chloroform and vortexed briefly. The supernatant (crude DNA) was collected for Southern blot analysis.

Identification of Orc/Cdc6 homologs in the haloarchaeal genomes

Fifteen haloarchaeal genomes were available through NCBI, including the H. hispanica genome sequenced by our laboratory[20]. When searching the Orc/Cdc6 homologs in these genomes, a BlastP search (BLOSUM62 matrix; 1 × 10-6 as an e-value cutoff) was performed against all haloarchaeal genomes (http://www.ncbi.nlm.nih.gov/sutils/genom_table.cgi) using the Orc/Cdc6 sequences from H. hispanica as seeds[47]. To focus on origin-associated Orc/Cdc6 homologs, PSI-BLAST (BLOSUM62 matrix, -e = 0.005) was also performed[47] using the profile from multiple alignments of experimentally functional Orc/Cdc6 homologs as a query. The results are summarized in Additional file3.

Prediction of orc/cdc6-associated replication origins

The IRs flanking orc/cdc6 genes were collected, and a motif (predicted ORB) search was performed using MEME software (motif size: 20–40; ZOOPS model)[48]. Consensus repeats in the IRs were confirmed using DNAMAN software (for windows, version 2.6)[49], and those harboring “G-string”[9] were considered as candidate ORB elements. The IRs were verified by hand, and only those contained inverted ORB repeats and were structurally similar to characterized archaeal replication origins were considered to be candidate orc/cdc6-associated replication origins. The results are summarized in Additional file4. Logo representation of ORB elements was performed using the program WebLogo (http://weblogo.berkeley.edu).

Phylogenetic analysis

16S rRNAs were collected from the 15 haloarchaeal genomes to estimate the evolutionary distance between them. The 16S rRNA sequence nearest the haloarchaeal-conserved replication origin (oriC1) was selected when there was more than one rRNA operon in the genome. Multiple alignments of the 16S rRNA sequences were performed using Clustal[50] implemented in MEGA[51]. A phylogenetic tree was constructed using neighbor-joining method[52] and maximum composite likelihood model implemented in MEGA, and 1000 bootstrap replicates were carried out. The Orc/Cdc6 homologs that were predicted to be associated with replication origins were collected from each of the 15 haloarchaeal genomes. The Orc/Cdc6 proteins, experimentally proven functional in their ability to recognize replication origins in other archaea (Pyrococcus abyssi[4], Sulfolobus solfataricus[5, 6], Aeropyrum pernix[7, 8]), were also included in this phylogenetic analysis. Multiple alignments of Orc/Cdc6 homologs were generated using Clustal (substitution matrix = BLOSUM; gap-opening penalty =10; gap-extension penalty = 0.1), and the result was adjusted manually to remove columns with many gaps. For maximum likelihood (ML) phylogeny, we used PHYML v3.0 with an LG substitution model and 100 nonparametric bootstrap replicates[53]. The data used to build the trees were deposited in TreeBASE (http://purl.org/phylo/treebase/phylows/study/TB2:S12601).

Comparative genomics and gene analysis

Whole genome alignments were performed using mummer and mummerplot algorithms in MUMmer[54] with default parameters. The GC plot was drawn using DNAplotter (window size: 50000; step size: 1000)[55]. Genome context analysis of the regions flanking the orc/cdc6-associated replication origins was performed using the NCBI Genome Workbench and scrutinized manually. Gene analysis was carried out using BlastP against the NCBI non-redundant proteins database (http://blast.ncbi.nlm.nih.gov/).

References

  1. Kelman LM, Kelman Z: Archaea: an archetype for replication initiation studies?. Mol Microbiol. 2003, 48 (3): 605-615. 10.1046/j.1365-2958.2003.03369.x.

    Article  CAS  PubMed  Google Scholar 

  2. Giraldo R: Common domains in the initiators of DNA replication in Bacteria, Archaea and Eukarya: combined structural, functional and phylogenetic perspectives. FEMS Microbiol Rev. 2003, 26 (5): 533-554. 10.1111/j.1574-6976.2003.tb00629.x.

    Article  CAS  PubMed  Google Scholar 

  3. Robinson NP, Bell SD: Origins of DNA replication in the three domains of life. FEBS J. 2005, 272 (15): 3757-3766. 10.1111/j.1742-4658.2005.04768.x.

    Article  CAS  PubMed  Google Scholar 

  4. Myllykallio H, Lopez P, Lopez-Garcia P, Heilig R, Saurin W, Zivanovic Y, Philippe H, Forterre P: Bacterial mode of replication with eukaryotic-like machinery in a hyperthermophilic archaeon. Science. 2000, 288 (5474): 2212-2215. 10.1126/science.288.5474.2212.

    Article  CAS  PubMed  Google Scholar 

  5. Robinson NP, Dionne I, Lundgren M, Marsh VL, Bernander R, Bell SD: Identification of two origins of replication in the single chromosome of the archaeon Sulfolobus solfataricus. Cell. 2004, 116 (1): 25-38. 10.1016/S0092-8674(03)01034-1.

    Article  CAS  PubMed  Google Scholar 

  6. Lundgren M, Andersson A, Chen L, Nilsson P, Bernander R: Three replication origins in Sulfolobus species: synchronous initiation of chromosome replication and asynchronous termination. Proc Natl Acad Sci USA. 2004, 101 (18): 7046-7051. 10.1073/pnas.0400656101.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  7. Robinson NP, Bell SD: Extrachromosomal element capture and the evolution of multiple replication origins in archaeal chromosomes. Proc Natl Acad Sci USA. 2007, 104 (14): 5806-5811. 10.1073/pnas.0700206104.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  8. Grainge I, Gaudier M, Schuwirth BS, Westcott SL, Sandall J, Atanassova N, Wigley DB: Biochemical analysis of a DNA replication origin in the archaeon Aeropyrum pernix. J Mol Biol. 2006, 363 (2): 355-369. 10.1016/j.jmb.2006.07.076.

    Article  CAS  PubMed  Google Scholar 

  9. Norais C, Hawkins M, Hartman AL, Eisen JA, Myllykallio H, Allers T: Genetic and physical mapping of DNA replication origins in Haloferax volcanii. PLoS Genet. 2007, 3 (5): e77-10.1371/journal.pgen.0030077.

    Article  PubMed Central  PubMed  Google Scholar 

  10. Coker JA, DasSarma P, Capes M, Wallace T, McGarrity K, Gessler R, Liu J, Xiang H, Tatusov R, Berquist BR, et al: Multiple replication origins of Halobacterium sp. strain NRC-1: properties of the conserved orc7-dependent oriC1. J Bacteriol. 2009, 191 (16): 5253-5261. 10.1128/JB.00210-09.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  11. Matsunaga F, Forterre P, Ishino Y, Myllykallio H: In vivo interactions of archaeal Cdc6/Orc1 and minichromosome maintenance proteins with the replication origin. Proc Natl Acad Sci USA. 2001, 98 (20): 11152-11157. 10.1073/pnas.191387498.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  12. Matsunaga F, Norais C, Forterre P, Myllykallio H: Identification of short 'eukaryotic' Okazaki fragments synthesized from a prokaryotic replication origin. EMBO Rep. 2003, 4 (2): 154-158. 10.1038/sj.embor.embor732.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  13. Duggin IG, McCallum SA, Bell SD: Chromosome replication dynamics in the archaeon Sulfolobus acidocaldarius. Proc Natl Acad Sci USA. 2008, 105 (43): 16737-16742. 10.1073/pnas.0806414105.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  14. Berquist BR, DasSarma S: An archaeal chromosomal autonomously replicating sequence element from an extreme halophile, Halobacterium sp. strain NRC-1. J Bacteriol. 2003, 185 (20): 5959-5966. 10.1128/JB.185.20.5959-5966.2003.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  15. Gaudier M, Schuwirth BS, Westcott SL, Wigley DB: Structural basis of DNA replication origin recognition by an ORC protein. Science. 2007, 317 (5842): 1213-1216. 10.1126/science.1143664.

    Article  CAS  PubMed  Google Scholar 

  16. Dueber EL, Corn JE, Bell SD, Berger JM: Replication origin recognition and deformation by a heterodimeric archaeal Orc1 complex. Science. 2007, 317 (5842): 1210-1213. 10.1126/science.1143690.

    Article  CAS  PubMed  Google Scholar 

  17. Bell SP, Dutta A: DNA replication in eukaryotic cells. Annu Rev Biochem. 2002, 71: 333-374. 10.1146/annurev.biochem.71.110601.135425.

    Article  CAS  PubMed  Google Scholar 

  18. Dueber EC, Costa A, Corn JE, Bell SD, Berger JM: Molecular determinants of origin discrimination by Orc1 initiators in archaea. Nucleic Acids Res. 2011, 39 (9): 3621-3631. 10.1093/nar/gkq1308.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  19. Roh SW, Nam YD, Nam SH, Choi SH, Park HS, Bae JW: Complete genome sequence of Halalkalicoccus jeotgali B3(T), an extremely halophilic archaeon. J Bacteriol. 2010, 192 (17): 4528-4529. 10.1128/JB.00663-10.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  20. Liu H, Wu Z, Li M, Zhang F, Zheng H, Han J, Liu J, Zhou J, Wang S, Xiang H: Complete genome sequence of Haloarcula hispanica, a model haloarchaeon for studying genetics, metabolism, and virus-host interaction. J Bacteriol. 2011, 193 (21): 6086-6087. 10.1128/JB.05953-11.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  21. Baliga NS, Bonneau R, Facciotti MT, Pan M, Glusman G, Deutsch EW, Shannon P, Chiu Y, Weng RS, Gan RR, et al: Genome sequence of Haloarcula marismortui: a halophilic archaeon from the Dead Sea. Genome Res. 2004, 14 (11): 2221-2234. 10.1101/gr.2700304.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  22. Pfeiffer F, Schuster SC, Broicher A, Falb M, Palm P, Rodewald K, Ruepp A, Soppa J, Tittor J, Oesterhelt D: Evolution in the laboratory: the genome of Halobacterium salinarum strain R1 compared to that of strain NRC-1. Genomics. 2008, 91 (4): 335-346. 10.1016/j.ygeno.2008.01.001.

    Article  CAS  PubMed  Google Scholar 

  23. Ng WV, Kennedy SP, Mahairas GG, Berquist B, Pan M, Shukla HD, Lasky SR, Baliga NS, Thorsson V, Sbrogna J, et al: Genome sequence of Halobacterium species NRC-1. Proc Natl Acad Sci USA. 2000, 97 (22): 12176-12181. 10.1073/pnas.190337797.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  24. Hartman AL, Norais C, Badger JH, Delmas S, Haldenby S, Madupu R, Robinson J, Khouri H, Ren Q, Lowe TM, et al: The complete genome sequence of Haloferax volcanii DS2, a model archaeon. PLoS One. 2010, 5 (3): e9605-10.1371/journal.pone.0009605.

    Article  PubMed Central  PubMed  Google Scholar 

  25. Malfatti S, Tindall BJ, Schneider S, Fahnrich R, Lapidus A, Labuttii K, Copeland A, Glavina Del Rio T, Nolan M, Chen F, et al: Complete genome sequence of Halogeometricum borinquense type strain (PR3). Stand Genomic Sci. 2009, 1 (2): 150-159.

    PubMed Central  PubMed  Google Scholar 

  26. Tindall BJ, Schneider S, Lapidus A, Copeland A, Glavina Del Rio T, Nolan M, Lucas S, Chen F, Tice H, Cheng JF, et al: Complete genome sequence of Halomicrobium mukohataei type strain (arg-2). Stand Genomic Sci. 2009, 1 (3): 270-277. 10.4056/sigs.42644.

    Article  PubMed Central  PubMed  Google Scholar 

  27. Bolhuis H, Palm P, Wende A, Falb M, Rampp M, Rodriguez-Valera F, Pfeiffer F, Oesterhelt D: The genome of the square archaeon Haloquadratum walsbyi : life at the limits of water activity. BMC Genomics. 2006, 7: 169-10.1186/1471-2164-7-169.

    Article  PubMed Central  PubMed  Google Scholar 

  28. Anderson I, Tindall BJ, Pomrenke H, Goker M, Lapidus A, Nolan M, Copeland A, Glavina Del Rio T, Chen F, Tice H, et al: Complete genome sequence of Halorhabdus utahensis type strain (AX-2). Stand Genomic Sci. 2009, 1 (3): 218-225. 10.4056/sigs.31864.

    Article  PubMed Central  PubMed  Google Scholar 

  29. Saunders E, Tindall BJ, Fahnrich R, Lapidus A, Copeland A, Del Rio TG, Lucas S, Chen F, Tice H, Cheng JF, et al: Complete genome sequence of Haloterrigena turkmenica type strain (4k). Stand Genomic Sci. 2010, 2 (1): 107-116. 10.4056/sigs.681272.

    Article  PubMed Central  PubMed  Google Scholar 

  30. Falb M, Pfeiffer F, Palm P, Rodewald K, Hickmann V, Tittor J, Oesterhelt D: Living with two extremes: conclusions from the genome sequence of Natronomonas pharaonis. Genome Res. 2005, 15 (10): 1336-1343. 10.1101/gr.3952905.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  31. Mongodin EF, Nelson KE, Daugherty S, Deboy RT, Wister J, Khouri H, Weidman J, Walsh DA, Papke RT, Sanchez Perez G, et al: The genome of Salinibacter ruber: convergence and gene exchange among hyperhalophilic bacteria and archaea. Proc Natl Acad Sci USA. 2005, 102 (50): 18147-18152. 10.1073/pnas.0509073102.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  32. Zhou L, Zhou M, Sun C, Xiang H, Tan H: Genetic analysis of a novel plasmid pZMX101 from Halorubrum saccharovorum: determination of the minimal replicon and comparison with the related haloarchaeal plasmid pSCM201. FEMS Microbiol Lett. 2007, 270 (1): 104-108. 10.1111/j.1574-6968.2007.00656.x.

    Article  CAS  PubMed  Google Scholar 

  33. Sun C, Zhou M, Li Y, Xiang H: Molecular characterization of the minimal replicon and the unidirectional theta replication of pSCM201 in extremely halophilic archaea. J Bacteriol. 2006, 188 (23): 8136-8144. 10.1128/JB.00988-06.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  34. Blow JJ, Ge XQ, Jackson DA: How dormant origins promote complete genome replication. Trends Biochem Sci. 2011, 36 (8): 405-414. 10.1016/j.tibs.2011.05.002.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  35. Woodward AM, Gohler T, Luciani MG, Oehlmann M, Ge X, Gartner A, Jackson DA, Blow JJ: Excess Mcm2-7 license dormant origins of replication that can be used under conditions of replicative stress. J Cell Biol. 2006, 173 (5): 673-683. 10.1083/jcb.200602108.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  36. Ge XQ, Jackson DA, Blow JJ: Dormant origins licensed by excess Mcm2-7 are required for human cells to survive replicative stress. Genes Dev. 2007, 21 (24): 3331-3341. 10.1101/gad.457807.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  37. Capes MD, Coker JA, Gessler R, Grinblat-Huse V, DasSarma SL, Jacob CG, Kim JM, DasSarma P, DasSarma S: The information transfer system of halophilic archaea. Plasmid. 2011, 65 (2): 77-101. 10.1016/j.plasmid.2010.11.005.

    Article  CAS  PubMed  Google Scholar 

  38. Mann S, Li J, Chen YP: Insights into bacterial genome composition through variable target GC content profiling. J Comput Biol. 2010, 17 (1): 79-96. 10.1089/cmb.2009.0058.

    Article  CAS  PubMed  Google Scholar 

  39. Ng WV, Ciufo SA, Smith TM, Bumgarner RE, Baskin D, Faust J, Hall B, Loretz C, Seto J, Slagel J, et al: Snapshot of a large dynamic replicon in a halophilic archaeon: megaplasmid or minichromosome?. Genome Res. 1998, 8 (11): 1131-1141.

    CAS  PubMed  Google Scholar 

  40. Ng WL, DasSarma S: Minimal replication origin of the 200-kilobase Halobacterium plasmid pNRC100. J Bacteriol. 1993, 175 (15): 4584-4596.

    PubMed Central  CAS  PubMed  Google Scholar 

  41. She Q, Brugger K, Chen L: Archaeal integrative genetic elements and their impact on genome evolution. Res Microbiol. 2002, 153 (6): 325-332. 10.1016/S0923-2508(02)01331-1.

    Article  CAS  PubMed  Google Scholar 

  42. Brugger K, Redder P, She Q, Confalonieri F, Zivanovic Y, Garrett RA: Mobile elements in archaeal genomes. FEMS Microbiol Lett. 2002, 206 (2): 131-141.

    Article  CAS  PubMed  Google Scholar 

  43. Rodriguez-Valera F: Approaches to prokaryotic biodiversity: a population genetics perspective. Environ Microbiol. 2002, 4 (11): 628-633. 10.1046/j.1462-2920.2002.00354.x.

    Article  CAS  PubMed  Google Scholar 

  44. Han J, Lu Q, Zhou L, Liu H, Xiang H: Identification of the polyhydroxyalkanoate (PHA)-specific acetoacetyl coenzyme A reductase among multiple FabG paralogs in Haloarcula hispanica and reconstruction of the PHA biosynthetic pathway in Haloferax volcanii. Appl Environ Microbiol. 2009, 75 (19): 6168-6175. 10.1128/AEM.00938-09.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  45. Cline SW, Lam WL, Charlebois RL, Schalkwyk LC, Doolittle WF: Transformation methods for halophilic archaebacteria. Can J Microbiol. 1989, 35 (1): 148-152. 10.1139/m89-022.

    Article  CAS  PubMed  Google Scholar 

  46. Liu H, Han J, Liu X, Zhou J, Xiang H: Development of pyrF-based gene knockout systems for genome-wide manipulation of the archaea Haloferax mediterranei and Haloarcula hispanica. J Genet Genomics. 2011, 38 (6): 261-269. 10.1016/j.jgg.2011.05.003.

    Article  CAS  PubMed  Google Scholar 

  47. Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997, 25 (17): 3389-3402. 10.1093/nar/25.17.3389.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  48. Bailey TL, Williams N, Misleh C, Li WW: MEME: discovering and analyzing DNA and protein sequence motifs. Nucleic Acids Res. 2006, 34 (Web Server issue): W369-W373.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  49. Woffelman C: DNAMAN for windows, version 2.6. Lynno Biosoft, Institute of Molecular Plant Sciences, Leiden University, the Netherlands, 1994-1996.

  50. Higgins DG, Sharp PM: CLUSTAL: a package for performing multiple sequence alignment on a microcomputer. Gene. 1988, 73 (1): 237-244. 10.1016/0378-1119(88)90330-7.

    Article  CAS  PubMed  Google Scholar 

  51. Tamura K, Dudley J, Nei M, Kumar S: MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0. Mol Biol Evol. 2007, 24 (8): 1596-1599. 10.1093/molbev/msm092.

    Article  CAS  PubMed  Google Scholar 

  52. Saitou N, Nei M: The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol. 1987, 4 (4): 406-425.

    CAS  PubMed  Google Scholar 

  53. Guindon S, Dufayard JF, Lefort V, Anisimova M, Hordijk W, Gascuel O: New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0. Syst Biol. 2010, 59 (3): 307-321. 10.1093/sysbio/syq010.

    Article  CAS  PubMed  Google Scholar 

  54. Delcher AL, Salzberg SL, Phillippy AM: Using MUMmer to identify similar regions in large sequence sets. Curr Protoc Bioinformatics. 2003, Chapter 10:10.3.1-10.3.18

    Google Scholar 

  55. Carver T, Thomson M, Bleasby A, Beeriman M, Parkhill J: DNAplotter: circular and linear interactive genome visualization. Bioinformatics. 2009, 25 (1): 119-120. 10.1093/bioinformatics/btn578.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

This work was partially supported by grants from the National Natural Science Foundation of China (30621005, 30830004, 30925001, 31100893) and the Chinese Academy of Sciences (KSCX2-EW-G-2-4).

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Hua Xiang.

Additional information

Competing interests

The authors declare no competing interests.

Authors’ contributions

ZW conducted the experiment, data analysis and drafted the manuscript. HL, JL and XL participated in the data collection and analysis. HX conceived of and coordinated the research and finalized the manuscript. All of the authors read and approved the final manuscript.

Electronic supplementary material

12864_2012_4475_MOESM1_ESM.tiff

Additional file 1: Physical mapping of eight predicted replication origins in H. hispanica . ORBs found in the IRs are indicated with numbered arrowheads, and the sequences are listed below. (TIFF 2 MB)

12864_2012_4475_MOESM2_ESM.tiff

Additional file 2: Screening of origin activity in H. hispanica . A. Schematic of the ARS assay. Δ: Corresponding origin (or cdc6 plus intergenic region)-deletion H. hispanica strains (unpublished data) were used for transformation to avoid plasmid integration. * For the two origins, oriC1-cdc6A and oriP-cdc6K, which cannot be knocked out from the chromosome and megaplasmid, respectively, the wide-type strains were used for transformation and Southern blot was performed to confirm ARS activity (Figure 1). B. ARS assay plates for eleven candidates. Colonies in plates of AS-168 (Mev) were observed after 7 days at 37 °C, and the minus signs (−) represent no visible colonies (no ARS activity). (TIFF 4 MB)

12864_2012_4475_MOESM3_ESM.doc

Additional file 3: Orc/Cdc6 homologues encoded in the haloarchaeal genomes. The complete set of Orc/Cdc6 homologues identified in the 15 sequenced haloarchaeal genomes. (DOC 341 KB)

12864_2012_4475_MOESM4_ESM.doc

Additional file 4: Predicted orc/cdc6 -associated replication origins in the haloarchaeal genomes. Prediction of ORB-containing replication origins directly adjacent to orc/cdc6 genes. The ORB elements are highlighted in red or blue colors. (DOC 250 KB)

12864_2012_4475_MOESM5_ESM.tiff

Additional file 5: Alignments of ORB elements in origin families of oriC1 , oriCa and oriCb . A, B and C respectively represent ORB elements found at origins belonging to origin families of oriC1, oriCa and oriCb, and conserved sequences are highlighted with shaded rectangles. (TIFF 1 MB)

12864_2012_4475_MOESM6_ESM.doc

Additional file 6: List of genes in the divergent regions between H. hispanica and H. marismortui. Blast analysis of the genes in the divergent regions between H. hispanica and H. marismortui. The genes whose closest relative is outside haloarchaea are highlighted in different colors. (DOC 212 KB)

12864_2012_4475_MOESM7_ESM.tiff

Additional file 7: Sequence similarity of regions around the oriCa and oriCb origins of replication in different haloarchaeal genomes. BLASTN analysis of the regions around the oriCa (A) and oriCb (B) origins of replication in different haloarchaeal genomes, and gray shading represents sequence similarity greater than 70%. (TIFF 322 KB)

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Wu, Z., Liu, H., Liu, J. et al. Diversity and evolution of multiple orc/cdc6-adjacent replication origins in haloarchaea. BMC Genomics 13, 478 (2012). https://doi.org/10.1186/1471-2164-13-478

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2164-13-478

Keywords