Skip to main content
  • ASM
    • Antimicrobial Agents and Chemotherapy
    • Applied and Environmental Microbiology
    • Clinical Microbiology Reviews
    • Clinical and Vaccine Immunology
    • EcoSal Plus
    • Eukaryotic Cell
    • Infection and Immunity
    • Journal of Bacteriology
    • Journal of Clinical Microbiology
    • Journal of Microbiology & Biology Education
    • Journal of Virology
    • mBio
    • Microbiology and Molecular Biology Reviews
    • Microbiology Resource Announcements
    • Microbiology Spectrum
    • Molecular and Cellular Biology
    • mSphere
    • mSystems
  • Log in
  • My alerts
  • My Cart

Main menu

  • Home
  • Articles
    • Latest Articles
    • COVID-19 Special Collection
    • Archive
    • Minireviews
  • Topics
    • Applied and Environmental Science
    • Clinical Science and Epidemiology
    • Ecological and Evolutionary Science
    • Host-Microbe Biology
    • Molecular Biology and Physiology
    • Therapeutics and Prevention
  • For Authors
    • Submit a Manuscript
    • Scope
    • Editorial Policy
    • Submission, Review, & Publication Processes
    • Organization and Format
    • Errata, Author Corrections, Retractions
    • Illustrations and Tables
    • Nomenclature
    • Abbreviations and Conventions
    • Publication Fees
    • Ethics Resources and Policies
  • About the Journal
    • About mBio
    • Editor in Chief
    • Board of Editors
    • AAM Fellows
    • For Reviewers
    • For the Media
    • For Librarians
    • For Advertisers
    • Alerts
    • RSS
    • FAQ
  • ASM
    • Antimicrobial Agents and Chemotherapy
    • Applied and Environmental Microbiology
    • Clinical Microbiology Reviews
    • Clinical and Vaccine Immunology
    • EcoSal Plus
    • Eukaryotic Cell
    • Infection and Immunity
    • Journal of Bacteriology
    • Journal of Clinical Microbiology
    • Journal of Microbiology & Biology Education
    • Journal of Virology
    • mBio
    • Microbiology and Molecular Biology Reviews
    • Microbiology Resource Announcements
    • Microbiology Spectrum
    • Molecular and Cellular Biology
    • mSphere
    • mSystems

User menu

  • Log in
  • My alerts
  • My Cart

Search

  • Advanced search
mBio
publisher-logosite-logo

Advanced Search

  • Home
  • Articles
    • Latest Articles
    • COVID-19 Special Collection
    • Archive
    • Minireviews
  • Topics
    • Applied and Environmental Science
    • Clinical Science and Epidemiology
    • Ecological and Evolutionary Science
    • Host-Microbe Biology
    • Molecular Biology and Physiology
    • Therapeutics and Prevention
  • For Authors
    • Submit a Manuscript
    • Scope
    • Editorial Policy
    • Submission, Review, & Publication Processes
    • Organization and Format
    • Errata, Author Corrections, Retractions
    • Illustrations and Tables
    • Nomenclature
    • Abbreviations and Conventions
    • Publication Fees
    • Ethics Resources and Policies
  • About the Journal
    • About mBio
    • Editor in Chief
    • Board of Editors
    • AAM Fellows
    • For Reviewers
    • For the Media
    • For Librarians
    • For Advertisers
    • Alerts
    • RSS
    • FAQ
Research Article | Ecological and Evolutionary Science

Evolutionary Trajectory of the Replication Mode of Bacterial Replicons

Bin-Bin Xie, Jin-Cheng Rong, Bai-Lu Tang, Sishuo Wang, Guiming Liu, Qi-Long Qin, Xi-Ying Zhang, Weipeng Zhang, Qunxin She, Yin Chen, Fuchuan Li, Shengying Li, Xiu-Lan Chen, Haiwei Luo, Yu-Zhong Zhang
Vaughn S. Cooper, Editor
Bin-Bin Xie
aState Key Laboratory of Microbial Technology, Marine Biotechnology Research Center, Shandong University, Qingdao, China
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Jin-Cheng Rong
aState Key Laboratory of Microbial Technology, Marine Biotechnology Research Center, Shandong University, Qingdao, China
bCollege of Marine Life Sciences and Frontiers Science Center for Deep Ocean Multispheres and Earth System, Ocean University of China, Qingdao, China
cLaboratory for Marine Biology and Biotechnology, Pilot National Laboratory for Marine Science and Technology (Qingdao), Qingdao, China
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Bai-Lu Tang
aState Key Laboratory of Microbial Technology, Marine Biotechnology Research Center, Shandong University, Qingdao, China
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Sishuo Wang
dSimon F. S. Li Marine Science Laboratory, School of Life Sciences, and State Key Laboratory of Agrobiotechnology, The Chinese University of Hong Kong, Shatin, Hong Kong SAR
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Guiming Liu
eBeijing Agro-Biotechnology Research Center, Beijing Academy of Agriculture and Forestry Sciences, Beijing, China
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Qi-Long Qin
aState Key Laboratory of Microbial Technology, Marine Biotechnology Research Center, Shandong University, Qingdao, China
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Xi-Ying Zhang
aState Key Laboratory of Microbial Technology, Marine Biotechnology Research Center, Shandong University, Qingdao, China
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Weipeng Zhang
bCollege of Marine Life Sciences and Frontiers Science Center for Deep Ocean Multispheres and Earth System, Ocean University of China, Qingdao, China
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Qunxin She
aState Key Laboratory of Microbial Technology, Marine Biotechnology Research Center, Shandong University, Qingdao, China
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Yin Chen
fSchool of Life Sciences, University of Warwick, Coventry, United Kingdom
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Fuchuan Li
aState Key Laboratory of Microbial Technology, Marine Biotechnology Research Center, Shandong University, Qingdao, China
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Shengying Li
aState Key Laboratory of Microbial Technology, Marine Biotechnology Research Center, Shandong University, Qingdao, China
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Xiu-Lan Chen
aState Key Laboratory of Microbial Technology, Marine Biotechnology Research Center, Shandong University, Qingdao, China
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Haiwei Luo
dSimon F. S. Li Marine Science Laboratory, School of Life Sciences, and State Key Laboratory of Agrobiotechnology, The Chinese University of Hong Kong, Shatin, Hong Kong SAR
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Haiwei Luo
Yu-Zhong Zhang
aState Key Laboratory of Microbial Technology, Marine Biotechnology Research Center, Shandong University, Qingdao, China
bCollege of Marine Life Sciences and Frontiers Science Center for Deep Ocean Multispheres and Earth System, Ocean University of China, Qingdao, China
cLaboratory for Marine Biology and Biotechnology, Pilot National Laboratory for Marine Science and Technology (Qingdao), Qingdao, China
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Yu-Zhong Zhang
Vaughn S. Cooper
University of Pittsburgh
Roles: Editor
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
DOI: 10.1128/mBio.02745-20
  • Article
  • Figures & Data
  • Info & Metrics
  • PDF
Loading

ABSTRACT

As typical bacterial replicons, circular chromosomes replicate bidirectionally and circular plasmids replicate either bidirectionally or unidirectionally. Whereas the finding of chromids (plasmid-derived chromosomes) in multiple bacterial lineages provides circumstantial evidence that chromosomes likely evolved from plasmids, all experimentally assayed chromids were shown to use bidirectional replication. Here, we employed a model system, the marine bacterial genus Pseudoalteromonas, members of which consistently carry a chromosome and a chromid. We provide experimental and bioinformatic evidence that while chromids in a few strains replicate bidirectionally, most replicate unidirectionally. This is the first experimental demonstration of the unidirectional replication mode in bacterial chromids. Phylogenomic and comparative genomic analyses showed that the bidirectional replication evolved only once from a unidirectional ancestor and that this transition was associated with insertions of exogenous DNA and relocation of the replication terminus region (ter2) from near the origin site (ori2) to a position roughly opposite it. This process enables a plasmid-derived chromosome to increase its size and expand the bacterium’s metabolic versatility while keeping its replication synchronized with that of the main chromosome. A major implication of our study is that the uni- and bidirectionally replicating chromids may represent two stages on the evolutionary trajectory from unidirectionally replicating plasmids to bidirectionally replicating chromosomes in bacteria. Further bioinformatic analyses predicted unidirectionally replicating chromids in several unrelated bacterial phyla, suggesting that evolution from unidirectionally to bidirectionally replicating replicons occurred multiple times in bacteria.

IMPORTANCE Chromosome replication is an essential process for cell division. The mode of chromosome replication has important impacts on the structure of the chromosome and replication speed. Bidirectional replication is the rule for bacterial chromosomes, and unidirectional replication has been found only in plasmids. To date, no bacterial chromosomes have been experimentally demonstrated to replicate unidirectionally. Here, we showed that the chromids (plasmid-derived chromosomes) in Pseudoalteromonas replicate either uni- or bidirectionally and that a single evolutionary transition from uni- to bidirectionality explains this diversity. These uni- and bidirectionally replicating chromids likely represent two stages during the evolution from a small and unidirectionally replicating plasmid to a large and bidirectionally replicating chromosome. This study provides insights into both the physiology of chromosome replication and the early evolutionary history of bacterial chromosomes.

INTRODUCTION

Prokaryotic DNA replication has been well studied in model bacteria, including the Gram-negative organism Escherichia coli and the Gram-positive organism Bacillus subtilis, both of which have a single circular chromosome that replicates bidirectionally (1–3). In general, this bidirectional replication is initiated at the origin (ori) site, after which two replication forks proceed in opposite directions and ultimately terminate in the terminus (ter) region, located roughly opposite the ori site on the circular chromosome (see Fig. S1a and b in the supplemental material). DnaA, the bacterial replication initiator protein, commences replication by binding to the ori site; the Tus protein terminates replication by binding to ter sites in E. coli (4). Therefore, with the ori and ter regions, the chromosome is divided into two halves, called replichores (5).

FIG S1

The chromosome structure for model bacteria Escherichia coli (a), Bacillus subtilis (b), and Vibrio cholerae (c). Download FIG S1, PDF file, 0.5 MB.
Copyright © 2021 Xie et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

About 10% of bacteria contain more than one chromosome (6, 7). As exemplified in Vibrio cholerae, the causative agent of cholera, the secondary chromosomes have a nucleotide composition close to that of the main chromosome and frequently carry core genes that are found on the (main) chromosome in other species (8). They often use plasmid-type maintenance and replication systems and are therefore referred to as the chromids (9). Similar to the main chromosomes, the chromids of Vibrio cholerae replicate bidirectionally: the replication is initiated from the origin of the chromid (ori2), and two replication forks proceed in the opposite directions (10) (Fig. S1c).

Chromosome replication exerts genome-wide mutational and selective pressures (11–15). The leading and lagging strands of the replication fork replicate differently, resulting in different mutational patterns (16, 17). One such strand-dependent compositional asymmetry is called GC skew, which means higher frequency of guanines (Gs) than cytosines (Cs) on the leading strand (16, 18). For bidirectional replication of typical bacterial circular chromosomes, GC skew shows a bipartite pattern with two reflection points, one located near ori and the other near ter. A study revealed that the secondary chromosome (Chr2) of the marine bacterium Pseudoalteromonas haloplanktis TAC125 does not show an otherwise expected GC skew, suggesting that it may be replicated unidirectionally (19).

Members of the genus Pseudoalteromonas are ubiquitous in a variety of marine habitats (20, 21). These bacteria utilize a wide spectrum of nutrients, including insoluble polysaccharides (22–24), proteins (25, 26), and bacterial cell wall (27). They are also known to produce a large variety of biologically active natural products (e.g., antimicrobial, antifouling, and algicidal substances) (20, 28). Compared to the oligotrophic paradigm, such as members of the SAR11 bacteria, which consistently have streamlined genomes and grow under exceedingly low nutrient conditions (29), Pseudoalteromonas spp. are typical copiotrophs which harbor large and variable genomes, explore nutrient patches, and interact with eukaryotic hosts (23, 30). Here, the Pseudoalteromonas Chr2 was employed as a model system to investigate the versatility and evolution of bacterial chromosome replication by using a combination of genome sequencing, phylogenomic and comparative genomic analyses, and experimental assays. Our study reveals a definitive example of how an evolutionary transition from unidirectional to bidirectional replication allows the increase of the sizes of chromids, thereby illustrating a process through which plasmids evolve into chromosomes.

RESULTS

Survey of replication directions using complete publicly available bacterial genomes.To search for bacterial chromosomes that might replicate unidirectionally, GC skew was calculated for all complete bacterial genomes available from the NCBI RefSeq database. By using the GC skew patterns of E. coli, B. subtilis, and V. cholerae chromosomes as the standard for bidirectional replication and that of the P. haloplanktis TAC125 chromid as the standard for unidirectional replication, replication directions were predicted for 16,103 large replicons (>200 kb) out of the total 26,371 replicons from 13,550 genomes. Among them, 2,304 were unpredictable and 13,722 were predicted to replicate bidirectionally and 77 unidirectionally (Fig. 1a; also, see Table S1). Results showed that the unidirectionally replicating large replicons may exist in different phyla, including Proteobacteria (Gammaproteobacteria and Alphaproteobacteria), Firmicutes, and Bacteroidetes (Table S1). Among these, the genus Pseudoalteromonas (order Alteromonadales of the class Gammaproteobacteria) was the major group that likely harbored a unidirectionally replicating secondary replicon (Fig. 1b; Table S1).

FIG 1
  • Open in new tab
  • Download powerpoint
FIG 1

Survey of replication directions for publicly available complete bacterial genomes. (a) Prediction of replication direction of large replicons (>200 kb) based on GC skew analyses. (b) Predicted unidirectionally replicating replicons. (c) Comparison of the ratio of Chr2 size to Chr1 size for bidirectionally and unidirectionally replicating Chr2 from bacteria with multiple chromosomes and Chr2 from Pseudoalteromonas.

TABLE S1

Predicted unidirectionally replicating replicons (>200 kb) from NCBI RefSeq database based on GC skew. Download Table S1, PDF file, 0.3 MB.
Copyright © 2021 Xie et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

Genome sequencing of Pseudoalteromonas spp.To investigate the unidirectional chromosome replication in Pseudoalteromonas, we collected type strains of the 26 species in this genus, which accounted for more than half of the currently recognized species of the genus. Using a combination of Illumina sequencing and PCR gap-closing, we obtained complete genome sequences for nine strains (Table S2). Similar to P. haloplanktis TAC125, all nine strains with complete genome sequences had two circular chromosomes, 1 and 2. For the 13 of the 17 strains with draft genome sequences, we obtained the complete sequences for the circular Chr2. The remaining four strains were also predicted to have a Chr2, because their genome assemblies contain contigs that showed good colinearity with the complete Chr2. Therefore, our study suggests that all Pseudoalteromonas spp. had two circular chromosomes.

TABLE S2

Detailed information for 26 Pseudoalteromonas type strains sequenced in this study. Download Table S2, PDF file, 0.3 MB.
Copyright © 2021 Xie et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

The replication system of Pseudoalteromonas.Genomic analyses revealed a conserved replication system that contains genes encoding the replication terminus site-binding protein (tus), helicase (dnaB), DNA polymerase III subunits (dnaE, dnaN, dnaQ, dnaX, holA, holB, and holC), gyrase subunit B (gyrB), primase (dnaG), and the single-strand DNA-binding protein (ssb), all of which show high sequence similarities to those from E. coli (Fig. 2; Table S3a). The chromosome segregation system contains the structural maintenance of chromosomes (SMC) protein and the accessory segregation and condensation proteins ScpA and ScpB, all of which had homologs in B. subtilis but not in E. coli (Table S3a). In addition, genes minCDE, encoding septum site-determining proteins, were also annotated. These proteins are generally involved in plasmid segregation in E. coli and B. subtilis (for a review, see reference 4). Thus, it is likely that they are involved in the segregation of Chr2 in Pseudoalteromonas. While most of these genes were located on Chr1 in Pseudoalteromonas, the tus and minCDE genes, which were located on the circular chromosome in E. coli, were invariably on the Chr2 of Pseudoalteromonas (Fig. 2). Localization of the tus gene on the chromid of Pseudoalteromonas was reported previously (31).

FIG 2
  • Open in new tab
  • Download powerpoint
FIG 2

Comparison of the main chromosome and chromid structure for species P. tunicata (a) and P. spongiae (b). The ori site (blue circle), dif site (red rectangle), and genes related to main chromosome/chromid replication and maintenance are indicated on the outer circle. GC skew is shown as the inner circle. Predicted replication directions are shown with blue (clockwise) and green (counterclockwise) arrows.

TABLE S3

(a) Genes involved in chromosome replication and their homologues in model organisms V. cholerae, E. coli, and B. subtilis. (b) Replication origin sites (ori) and termination sites (dif) on the main chromosomes and chromids. Download Table S3, PDF file, 0.5 MB.
Copyright © 2021 Xie et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

Consistent with the fact that Pseudoalteromonas spp. have two chromosomes, two chromosome partitioning systems were identified in their genomes: parA1B1 on Chr1 and parA2B2 on Chr2. The parA2B2 genes were located adjacent to the putative ori2 site in a conserved gene cluster, which also contained the tus gene, the putative replication initiator protein gene repA, the possible DNA helicase gene ywqA, and the DNA recombination protein gene rmuC (Fig. 2). However, the functions of the repA gene and the replication origin ori2 have not been tested experimentally. Here, we conducted a series of truncation experiments to identify the minimal replicon of the Pseudoalteromonas chromid using Pseudoalteromonas spongiae JCM 12884T and Pseudoalteromonas tunicata DSM14096T as representatives of the genus. Our results showed that repA together with the putative ori2 site can initiate and are essential to the replication of a replicon lacking the replication initiator protein, confirming the function of RepA as the replication initiator protein and the predicted ori2 site as the origin of Chr2 (Fig. 2 and Fig. S3). A sequence search against the NCBI nonredundant (nr) protein database with RepA as the query revealed no hits with significant identities in non-Pseudoalteromonas bacteria or plasmids, although it was suggested that it may be related to the initiator protein of the R1 plasmid (19). RepA also showed no detectable sequence similarity to the replication initiator protein DnaA encoded by Chr1. This result indicates that the Chr2 was likely derived from an unknown plasmid and, therefore, that Chr2 is a chromid. In sum, Pseudoalteromonas Chr1 and Chr2 have evolved different strategies to initiate the replication (DnaA+ori1 versus RepA+ori2). As in Pseudoalteromonas, in V. cholerae the initiation of chromid replication is also different from that of the main chromosome, which uses DnaA to initiate the main chromosome replication and RctB to initiate the chromid replication (32).

Unidirectional replication prevails among Pseudoalteromonas chromids.Similar to the chromid of the strain TAC125, most Pseudoalteromonas chromids showed a GC skew different from that of a typical chromosome (Fig. 2a, right, and Fig. S4a), suggesting that unidirectional replication was prevalent among members of the Pseudoalteromonas genus. Another piece of evidence for unidirectional replication was related to the location of the dif site, which participates in chromosome segregation and often locates near the ter region (31, 33, 34). For the main chromosome, which replicated bidirectionally, the dif site (dif1) was located at a position roughly opposite ori (ori1) (Fig. 2 and Fig. S1). However, for most Pseudoalteromonas chromids (Table S3b), dif2 was located adjacent to the ori2 site (only approximately 4 to 27 kb away).

To experimentally determine the replication directions, we performed deep genome sequencing and analyzed the sequencing depth along the replicating main chromosomes and chromids. The underlying principle is that, once chromosome replication is initiated, the chromosomal region that has been replicated (i.e., near the ori site) should have an additional copy compared to the region that has not been replicated (i.e., near the ter site) (10, 35–37). Technically, this will be manifested as a higher coverage (i.e., read counts) for nucleotides near ori than those closer to ter (38). Our data showed that, at the exponential phase, the coverage decreased along Chr1 from ori to dif on both sides of dif (or on both replichores) (Fig. 3a and b), indicating that Chr1 was replicated bidirectionally. However, this symmetric pattern was missing from the P. tunicata chromid, indicating a unidirectional replication of this chromid (Fig. 3a). Similarly, deep sequencing also revealed unidirectional replication in seven other strains tested (Fig. S5a). This is the first experimental demonstration of the unidirectional replication of a bacterial chromosome.

FIG 3
  • Open in new tab
  • Download powerpoint
FIG 3

Analyses of replication direction and synchronization mode by deep genome sequencing. (a and b) Gradient changes of sequencing coverage from the replication origin (ori site) to the replication terminus (represented by dif) along a replicating replicon for P. tunicata DSM 14096T (a) and P. spongiae JCM 12884T (b) at the exponential phase. Data were presented in bins of 1,000 bp and were corrected using data from the stationary phase (see Materials and Methods). Note that the y axis was set to a base 2 logarithmic scale. Arrows indicate the replication direction (blue for clockwise and green for counterclockwise). (c) Sequencing coverage for the replication origin regions (ori1 and ori2) and the terminus regions (ter1 and ter2). Sequencing coverage of ter1 and ter2 was represented by the coverage near the dif1 and dif2 sites, respectively.

Presence of bidirectional replication in Pseudoalteromonas chromids.Interestingly, among the Pseudoalteromonas species sequenced, P. spongiae harbored a chromid that showed a GC skew similar to that of the E. coli, B. subtilis, and V. cholerae chromosomes (Fig. 2b, right, and Fig. S1), suggesting that the P. spongiae chromid may replicate bidirectionally. Consistent with this hypothesis, the dif2 site was located approximately opposite the ori2 site on its chromid (Fig. 2b). Further support for bidirectional replication of the P. spongiae chromid was obtained from the deep-sequencing result, which showed a coverage profile highly similar to that of the bidirectionally replicating Chr1 (Fig. 3b). In addition, we collected experimental and/or bioinformatics evidence for bidirectional replication of the chromids in another two Pseudoalteromonas strains, including the strain SAO4-4 (39) (Fig. S4b) from our culture collection, as confirmed with the deep genome sequencing (Fig. S5b), and P. piratica OCN003T (40) from the NCBI GenBank database (Fig. S4b). In sum, a total of only three strains among the currently available Pseudoalteromonas members carried bidirectionally replicating chromids, indicating that this replication mode appears to be less common than unidirectional replication among the Pseudoalteromonas chromids.

Evolutionary history of the replication mode of Pseudoalteromonas chromids.To infer the evolutionary history of the replication direction of the chromids, we first attempted to construct a reliable species tree delineating the evolutionary history of the available Pseudoalteromonas members using phylogenomic approaches. Ideally, this species tree should be built with sequences from both the chromosome and the chromid of each strain. Unlike the main chromosome tree, where close relatives of Pseudoalteromonas were available to root the tree, however, there were no suitable outgroup sequences for the chromids. As a result, the phylogeny of the chromid was rooted using the less accurate strategy called midpoint rooting (Fig. S2b, d, and f). The identical topologies between the main chromosome tree (Fig. S2a, c, and e) and the chromid tree (Fig. S2b, d, and f) indicate that the main chromosome and the chromid shared the same evolutionary history since their origins. We did notice that one of the chromid trees (Topo 2) (Fig. S2d) differed from all other trees (Topo 1) (Fig. S2). We therefore statistically compared the two topologies (Topo 1 versus Topo 2) with sequences from the main chromosome and showed that Topo 1 was unanimously supported (Fig. S2g). Taken together, these lines of evidence support the idea that the main chromosome tree can represent the species tree (Fig. 4a).

FIG 4
  • Open in new tab
  • Download powerpoint
FIG 4

Phylogenomics analysis of Pseudoalteromonas spp. and comparison of the replicating directions of the chromids of different clades. (a) A maximum-likelihood tree reconstructed based on the concatenated amino acid sequences of shared single-copy genes on the main chromosome. The genus diverged into four major clades: I, II, III, and IV. The replication direction of the chromid was marked for all clades and their ancestors. All branches are supported by the IQ-Tree ultrafast bootstrap percentages of 100%. The bar represents 0.05 substitution per site. See details in Fig. S2a for the complete tree. (b) Colinear analysis of uni- and bidirectionally replicating chromids of Pseudoalteromonas spp. A total of 82 common single-gene families located on the chromids were found among the 43 genomes studied and are indicated using lines between chromids. Genes parA2 and parB2 and the locus ori2 are shown in blue. Genes repA and tus are shown in green. The ter2 region is represented by dif2. The locus dif2 is shown in red. The remaining genes are shown in gray. It is clear that the neighboring region of dif2 is highly conserved. Arrows indicate the replication directions. Asterisks indicate the predicted insertion sequences and integrons, and boxes indicate the predicted prophage-like regions.

FIG S2

Maximum-likelihood trees based on the concatenated amino acid sequences of single-copy genes on the main chromosome (a, c, and e) and chromid (b, d, and f) based on different alignment trimming strategies and the removal of compositionally heterogeneous genes. (g) Topology tests on two different topologies for the chromid trees. Download FIG S2, PDF file, 0.9 MB.
Copyright © 2021 Xie et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

FIG S3

Determination of the minimal replicon for the chromids. Download FIG S3, PDF file, 0.9 MB.
Copyright © 2021 Xie et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

FIG S4

Unidirectional (a) and bidirectional (b) replication of Pseudoalteromonas chromids revealed by GC skew analyses. Download FIG S4, PDF file, 1.1 MB.
Copyright © 2021 Xie et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

FIG S5

(a) Unidirectional replication revealed by deep sequencing at the exponential phase for seven species. (b) Bidirectional replication revealed by deep sequencing at the exponential phase for the strain SAO4-4. Download FIG S5, PDF file, 0.5 MB.
Copyright © 2021 Xie et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

As shown in the species tree, the available Pseudoalteromonas genomes comprised two monophyletic groups, each consisting of two clades (clades I and II versus clades III and IV), with clade I consisting exclusively of the three strains carrying the bidirectionally replicating chromids. With this robust species tree, we applied the parsimonious rule of evolution to infer the evolutionary history of the chromid replication mode. Because a unidirectional replication mode of the chromid at the last common ancestor (LCA) of all Pseudoalteromonas requires only one evolutionary change (from uni- to bidirectional replication), whereas a bidirectional mode at the LCA necessitates a minimum of two independent changes (both from bi- to unidirectional replication) to explain the phylogenetic distribution of the chromid replication mode of the extant species, the former was favored (Fig. 4a). Molecular dating analysis suggests that the LCA of Pseudoalteromonas arose around 502 to 378 million years ago (MYA) (Fig. S6). Clade I diverged from other lineages 469 to 324 MYA, though the extant Pseudoalteromonas members displaying bidirectional replication at the chromid evolved from a more recent ancestor, about 36 to 9 MYA (Fig. S6).

FIG S6

Estimated divergence time of Pseudoalteromonas. Download FIG S6, PDF file, 0.6 MB.
Copyright © 2021 Xie et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

Colinear analysis (Fig. 4b) suggests that different positions of ter2 relative to ori2 in the genus can be explained by gene insertion events between them, resulting in a shift of ter2 away from ori2. Consistently, genome annotation revealed the presence of a number of insertion sequences and integrons in chromids (Fig. 4b, asterisks; Table S4a). Perhaps the most striking feature was that chromids contain long prophage-like regions (57 to 277 kb in total), accounting for 6.5% to 24.1% of the chromid length (Fig. 4b, boxes; Table S4a). Compared to chromids, main chromosomes had much lower ratios of prophage-like regions (1.1% to 6.7%; 36 to 258 kb in length). This result suggests that the insertion of exogenous DNA in chromids occurred at high frequency and that such events may have significantly changed the position of ter in the LCA of clade I. Thus, we identified an example showing how insertions of DNA in a chromid led to a switch from unidirectional to bidirectional replication.

TABLE S4

(a) Integrons, insertion sequences, and prophages predicted on Pseudoalteromonas chromids. (b) Homologs of dif2 sites on chromids (Chr2) of P. spongiae JCM12884T, Pseudoalteromonas sp. SAO4-4, P. piratica OCN003T, and P. prydzensis DSM14232T predicted by BLASTN searches (E value cutoff of 1). Download Table S4, PDF file, 0.7 MB.
Copyright © 2021 Xie et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

Replications of the two chromosomes are coordinated to terminate in synchrony.The appropriate timing of replication is essential for successful cell division (10). Unlike in bacteria with a single chromosome, the presence of multiple chromosomes in a bacterial genome necessitates a cellular mechanism that can coordinate the replication of multiple chromosomes (6). Two conceptual models have been proposed for such coordination: initiation synchrony (control of initiation time) and termination synchrony (control of termination time) (41). Here, we compared P. tunicata DSM14096T with P. spongiae JCM 12884T to investigate the coordination modes used by unidirectionally versus bidirectionally replicating chromids. Again using deep-sequencing data, we analyzed the sequencing coverage for the main chromosome ori and dif sites (ori1 and dif1) and the chromid ori and dif sites (ori2 and dif2). As expected, during the stationary growth phase, a period lacking extensive replication, the sequencing coverage for each ori and dif site was nearly equal in both species (Fig. 3c). At the exponential phase, however, the ori1/ori2 ratios were 1.5 and 1.8 for P. spongiae and P. tunicata, respectively, while their dif1/dif2 ratios each were 1.0 (Fig. 3c), indicating that both Pseudoalteromonas species we examined adopted termination synchrony to coordinate the replication of the multiple chromosomes in their genomes, despite the fact that they differed in the replication mode (i.e., unidirectional versus bidirectional replications) of their chromids.

Bidirectionally replicating chromids are larger than unidirectionally replicating chromids.Yet another striking trend between the bidirectionally and unidirectionally replicating Pseudoalteromonas chromids was that the bidirectionally replicating chromids were larger (1.60 to 1.67 Mb versus 0.64 to 1.39 Mb). The coordination of the replication of multiple chromosomes implies that the time required for the complete replication of the main chromosome sets an upper time limit for the chromid replication (10, 42). Compared with unidirectional replication, bidirectional replication needs only half the replication time and therefore allows the chromid size to increase by 1-fold while keeping the same replication time. This was further supported by higher ratios of chromid size to main-chromosome size for bidirectional replication than for unidirectional replication in both the genus Pseudoalteromonas (50.6% to 50.8% versus 19.3% to 29.8%) and other bacteria (50.4 ± 25.2% versus 18.8 ± 8.7%) (Fig. 1c). It was also noted that the ratios 50.6% to 50.8% are too high to allow the chromids to finish the replication in unidirectional mode within the bidirectional replication time of the main chromosome. In this scenario, using a bidirectional replication system would allow the replication of larger chromids. The increase in the chromid and thus genome size may increase the metabolic versatility, which is particularly beneficial to copiotrophic bacteria like Pseudoalteromonas. Thus, conversion from a unidirectional to a bidirectional replication system may confer a competitive advantage for Pseudoalteromonas to explore multiple niches. Consistently, a survey of the publicly available multipartite genomes revealed that 63% of bidirectionally replicating secondary replicons have a size larger than one half of the corresponding main chromosome (Fig. 1c), though other unknown mechanisms may also play roles in the evolution of replication direction.

DISCUSSION

Bacterial circular chromosomes replicate bidirectionally by the theta mechanism. Circular plasmids replicate by three general mechanisms: theta type, strand displacement, and rolling circle (43). Different machineries are involved in different plasmid replication mechanisms (43). For example, initiation of theta-type plasmid replication generally requires a plasmid-encoded initiator protein. In contrast, initiation of strand displacement replication requires the plasmid-encoded helicase, primase, and initiator protein, and initiation of rolling-circle replication requires the plasmid-encoded protein with DNA strand transferase enzymatic activity. The absence of the above essential proteins for the latter two mechanisms indicated that Pseudoalteromonas chromids replicate by the theta mechanism. Consistent with this, no such proteins in chromids from other bacteria have been reported. Considering the different replication factors involved in different mechanisms, a parsimony assumption is that chromids likely evolved from plasmids that also used the theta-type replication. Unlike in chromosomes, however, theta-type replication in plasmids is mostly unidirectional (43). Therefore, it is possible that bidirectionally replicating chromids may have evolved from unidirectionally replicating plasmids. In this situation, the unidirectionally replicating chromids may represent an intermediate stage in the evolutionary trajectory from a unidirectionally replicating plasmid to a bidirectionally replicating chromosome.

Though chromids have been observed in different bacteria, only bidirectionally replicating chromids have been experimentally verified so far (41). Previous bioinformatic analysis predicted that the chromid of P. haloplanktis TAC125 may replicate unidirectionally (19). Here, the existence of unidirectional replication was experimentally verified in Pseudoalteromonas. Since a bacterial chromosome typically replicates bidirectionally and is larger than a chromid, the question of whether a small and unidirectionally replicating chromid (more plasmid-like) may evolve into a large and bidirectionally replicating one (more chromosome-like) was raised. Deep genome sequencing and phylogenomic analysis revealed that the last common ancestor of Pseudoalteromonas chromids replicated unidirectionally and the unidirectional replication evolved into bidirectional replication only once, occurring at one of the basal clades of the Pseudoalteromonas phylogeny. Therefore, profiling of Pseudoalteromonas from marine environments ultimately revealed a basic biological mechanism that deepens our understanding of the early evolutionary history of bacterial chromosomes.

In E. coli, the replication fork is blocked by complexes formed by binding of the Tus protein to the ter regions (44). The existence of a conserved tus gene within the Pseudoalteromonas genomes suggests that the genus also employs the Tus/ter replication fork trap system. However, with the ter consensus sequence of E. coli as the reference (44), a sequence search failed to identify homologues near the expected ter region on the P. tunicata chromid (i.e., the region with the lowest sequencing coverage as in Fig. 3a), suggesting that the ter sequences of Pseudoalteromonas are either highly divergent from or evolutionarily unrelated to that of E. coli. Two replication forks are formed at ori2 on a bidirectionally replicating chromid. It is unclear whether one or two forks are formed on a unidirectionally replicating chromid. In the former situation, the change of replication initiation machinery occurred during the transition from uni- to bidirectional replication to enable the formation of a second replication fork. In the latter, two replication forks are formed at ori2 on a unidirectionally replicating chromid, and the one proceeding in the counterclockwise direction is immediately blocked at the ter2 region, which is located very close to the ori2 site on the unidirectionally replicating chromid. The molecular mechanism of the initiation of unidirectional replication needs further study.

The chromids contained long prophage-like regions that accounted for up to 6.5% to 24.1% of the chromid length, suggesting that insertion of phage genomes may have significantly changed the chromid structure. It was proposed that the insertion events between the ori2 and dif2 sites in the unidirectionally replicating ancestor moved dif2 away from ori2, eventually resulting in bidirectionally replicating descendants. Consistently, prophages were annotated at or near the dif2 site in a few strains (Fig. 4, boxes). Genetic studies of V. cholerae have identified a number of phages and other integrative mobile elements that hijack tyrosine recombinases XerC and XerD from their host for integration (also referred to as IMEXs) (for a review, see reference 45). In most bacteria, XerC and XerD are used to resolve chromosome dimers by binding to dif (4, 46). A sequence search with the dif2 sequence as the query revealed a few dif2-like sequences within or near the prophage regions (Table S4b). It was noted that the dif2 site was located within a prophage in the bidirectionally replicating chromids but not in the unidirectionally replicating chromids (Fig. 4b; see Table S3b for the positions of the dif2 sites and Table S4a for the positions of the prophages). Prophages were not found between the dif2 and the ori2 sites in the unidirectionally replicating chromids either. This observation was consistent with the hypothesis that a few phages and/or other mobile genetic elements might have specifically integrated at the dif2 site, resulting in the separation of dif2 from ori2 as well as the increase of the chromid size.

The switch from uni- to bidirectional replication of the Pseudoalteromonas chromids represents a simple strategy that bacteria use to add (a large number of) genes to the chromosome while exerting little influence on the replication physiology of host bacteria. This is because replication initiation, elongation, and termination of chromids are not changed in spite of the fact that hundreds of genes are acquired on the two replichores of a chromid. Provided that the cell harboring the new bidirectionally replicating chromid adopted termination synchrony to coordinate replication, each DNA replication complex would have to replicate only half as much DNA in the allotted replication time period. Thus, there would be more than adequate time for the replication of additional genetic material that was inserted into the now bidirectionally replicating chromid, explaining the trend observed among the Pseudoalteromonas genomes that the bidirectionally replicating chromids are larger. While bidirectional replication of the chromid may benefit the organism, the reality is that unidirectional replication is more widespread within the genus. This suggests that other factors may also control the replication mode in this genus, which is a topic of interest for future studies.

The replication of multiple chromosomes within a bacterial cell is well coordinated, and the underlying mechanism has been studied with V. cholerae as a model. The V. cholerae chromid replication is regulated by the initiator RctB encoded by the chromid (Fig. S1c), which binds to iterons (12-mer sites) in ori2 to promote initiation and binds to 39-mer regulatory sites in ori2 to strongly inhibit initiation (47). RctB can also bind to a 150-bp site on the main chromosome, an ori2 replication enhancer, which increases the RctB binding affinity for iterons and decreases affinity for 39-mers in the chromid (48). Replication of this 150-bp site (also named crtS) triggers the initiation of chromid replication (10). In spite of the differences in chromid size and mode of replication across the genus Pseudoalteromonas, the replication of the two chromosomes in all strains studied was precisely coordinated to terminate in synchrony (Fig. 3 and Fig. S5). This observation suggests that the initiation of chromid replication is controlled by the timing of the main chromosome replication in Pseudoalteromonas. Though Pseudoalteromonas and V. cholerae use different initiators for chromid replication, it is possible that Pseudoalteromonas uses a similar molecular mechanism in the communication between chromosomes, i.e., the replication of a specific site on the main chromosome to trigger the initiation of chromid replication. It was not unexpected that search of the main chromosome sequence of Pseudoalteromonas revealed no sites showing similarities to crtS carried in V. cholerae, which is probably a result of the different initiators used by the two organisms. The molecular mechanism underlying the synchrony of the replication of multiple chromosomes in Pseudoalteromonas is worth further study.

The mode of replication is known to strongly impact the organization of bacterial chromosomes (11, 49–51). The replication mode leads to the base composition asymmetries on the chromosome, and one such asymmetry, GC skew, has been used to predict the position of the ori site (52). Here, experimental verification of the existence of unidirectional replication in the Pseudoalteromonas chromid suggests that the GC skew is a valuable predictor of the replication mode of bacterial chromosomes. Other demonstrated replication-dependent structural features include unequal distributions of essential genes between the two replication strands, a propensity for highly expressed genes to be clustered near the ori site, and the accumulation of recombination-associated chi sites in the leading strand (50). The current study presents an example of the opposite situation wherein a structural change shifts the mode of replication used by a chromosome. Furthermore, our study revealed a higher number of genes on the leading strand than the lagging strand for both replichores of the Pseudoalteromonas main chromosomes (leading-lagging ratios of 1.61 ± 0.07 and 1.53 ± 0.07 for the right and left replichores, respectively; n = 9). Similarly, there were more genes on the leading strand for the unidirectionally replicating chromids (1.34 ± 0.10; n = 21) and for the right replichores of the bidirectionally replicating chromids (1.26 ± 0.13; n = 3). However, such unequal distribution of genes on the two strands was not apparent on the left replichore of the bidirectionally replicating chromids (1.04 ± 0.02; n = 3), consistent with the hypothesis that the left replichore was newly generated from exogenous DNA insertion between the ori2 and dif2 sites. Replication also has a global impact on the mutation rate. Generally, there are higher mutation rates in late-replicating regions (14), and the secondary chromosome shows higher mutation rates than the main chromosome (15). Recent studies of genomes possessing one single chromosome and multiple chromosomes revealed that the correlation between the mutation rate and replication timing is not simple but periodic (12, 13, 53–56). The genus Pseudoalteromonas contains both unidirectionally and bidirectionally replicating members, representing a new model to study the mutation rate variation across chromosome.

With the marine bacterial genus Pseudoalteromonas as the model, our study experimentally demonstrated the existence of unidirectional replication in bacterial chromosomes and revealed both uni- and bidirectional replication modes in the homologous chromids from different strains. Our results suggest that the uni- and bidirectionally replicating chromids represent two stages in the evolutionary history from a unidirectionally replicating plasmid to a bidirectionally replicating chromosome. Our bioinformatic analyses of unrelated bacterial lineages predicted that, while most additional large replicons replicate bidirectionally, a few may also replicate unidirectionally, suggesting that the evolution from unidirectionally replicating plasmids to chromosomes (chromids) may have occurred multiple times in bacteria.

MATERIALS AND METHODS

Genome sequencing and assembly.Type strains for 26 Pseudoalteromonas species were purchased from the Leibniz Institute DSMZ-German Collection of Microorganisms and Cell Cultures (strains 14585, 9414, 6820, 17587, 16473, 18437, 26439, 14402, 8810, 6065, 15925, 9166, 14232, 15203, 6057, 8771, 6061, 6842, 14001, 14401, 15557, and 14096) and the Japan Collection of Microorganisms (strains 15903, 21460, 20779, and 12884). The complete list of strains is presented in Table S2.

Genomic DNA was extracted using the PowerMax soil DNA isolation kit (Mo Bio laboratories, Inc., USA) according to the manufacturer’s instructions. Multiple Illumina paired-end DNA libraries of different insert sizes (500 bp, 2 kbp, 5 kbp, and 10 kbp) were prepared for each strain (Table S2). Genome sequencing was performed using the Illumina HiSeq 2000 platform. The read length was 90 bp. Clean reads were obtained by removing those containing ≥36-bp low-quality bases (Phred score, ≤20), those containing ≥9 Ns, and those containing adapter contamination with FASTX-Toolkit version 0.0.13 (http://hannonlab.cshl.edu/fastx_toolkit/). The resultant clean reads were assembled using SOAPdenovo version 1.05 (57), with reads from the 500-bp paired-end library being used to create contigs and reads from 2-kb, 5-kb, and 10-kb libraries used to construct scaffolds (pair_num_cutoff = 5, map_len = 80) (Table S2). A few intrascaffold gaps were closed using reads from the 500-bp paired-end library by GapCloser version 1.1 (http://soap.genomics.org.cn/soapdenovo.html). The remaining intra- and interscaffold gaps were closed using PCR. For the intrascaffold gaps, PCR primer pairs were designed on the flanking contig ends. To design primer pairs for closing the interscaffold gaps, the linkage information between scaffolds was predicted by comparison with genomes of strains TAC125 and SM9913 as well as other species sequenced in this study, using MUMmer version 3.23 (58). Sequences of PCR products were determined by directed sequencing (primer walking). Scaffolds were split into contigs at unclosed gaps. The resultant assemblies were used for further annotation. Detailed information on sequencing and assembly is available in Table S2.

Open reading frames were predicted using Glimmer version 3.02 (59). Gene functions were annotated using both the NCBI nonredundant protein database and the Cluster of Orthologous Groups of proteins (COG) database (60) using BLASTP (E value cutoff, 1e−5; score cutoff, 60; coverage cutoff, 50% for both query and target; and identity cutoff, 35%).

GC skew analysis.The GC skew is calculated as (number of Gs – number of Cs)/(number of Gs + number of Cs), with a sliding window. The cumulative GC skew was calculated by summing up the GC skew values for all windows from the start window to the current window along the replicon. The circular representation of GC skew was produced using DNAPlotter (61).

To explore the replication direction of bacterial chromosomes, all complete bacterial genomes in the NCBI RefSeq database (May 2019 version) were downloaded. Then, GC skew for all large replicons (>200 kb) was calculated using a self-written Perl script, with a window size of 1 kb and a step size of 1 kb. Lastly, the GC skew and cumulative GC skew were manually inspected. The chromosomes of E. coli, B. subtilis, and V. cholerae were used as standards for bidirectional replication. The chromid of P. haloplanktis TAC125 was used as a standard for unidirectional replication.

In silico analysis of ori and dif sites.The ori sites were annotated based on the sequence similarities with the reference sequences, the locus of the neighboring gene dnaA (for ori1) or repA (for ori2), and the internal structures, including DnaA boxes and Dam methylation sites (52). With the ori sites (ori1 and ori2) of TAC125, V. cholerae, and E. coli as references, genome sequences were searched using BLASTN with an E value cutoff of 1e−5, a score cutoff of 60, a coverage cutoff of 50% (for both query and target), and an identity cutoff of 30%. The 28-bp dif sites were searched with dif sites of TAC125, V. cholerae, and E. coli as references using BLASTN with an E value cutoff of 1 and a query coverage cutoff of 40%. The obtained candidate dif sites were further confirmed based on the following criteria: located in noncoding sequences and roughly opposite or near the ori site and possessing a palindromic structure and a highly conserved XerD binding site (62).

Characterization of the minimal replicon.The minimal replicon for the chromid was characterized by testing the ability of the ori2 site and its flanking sequences to initiate the replication of the vector pOriT-4CM in Pseudoalteromonas. pOriT-4CM is highly similar to the shuttle vector pOriT-4Em, which can replicate in both E. coli and Pseudoalteromonas (63). The only difference between the two vectors is that pOriT-4Em contains an erythromycin resistance (Emr) gene, while pOriT-4CM contains a chloramphenicol resistance (Cmr) gene (see Fig. S3a for map of pOriT-4CM). First, the vector pOriT-4CM was digested with BamHI and PstI to delete the sequence for replication initiation in Pseudoalteromonas. Next, a set of sequences including the ori2 site and the flanking sequences from the chromid were amplified by PCR and were then ligated into the linearized pOriT-4CM by using an In-Fusion HD Plus cloning kit (TaKaRa, Japan) to obtain a series of recombinant plasmids (p4CM-pspo-1, p4CM-pspo-2, p4CM-pspo-3, p4CM-pspo-4, p4CM-ptun-1, p4CM-ptun-2, p4CM-ptun-3, and p4CM-ptun-4). These plasmids were transformed into E. coli DH5α cells (Trans, China) and then were isolated and verified by Sanger DNA sequencing. The verified plasmids were then introduced into conjugation donor strain E. coli WM3064. Next, the plasmids were transferred into P. spongiae JCM 12884T or P. tunicata DSM14096T by intergeneric conjugations, as described previously (64). P. spongiae and P. tunicata were grown on marine LB agar medium with 100 μg/ml ampicillin and 35 μg/ml chloramphenicol (Sangon Biotech, China) at 25°C to screen single colonies that contain the recombinant plasmids. Each single colony was then grown in marine LB broth supplemented with 100 μg/ml ampicillin and 35 μg/ml chloramphenicol, and the recombinant plasmids were extracted and verified by restriction enzyme digestion. The bacterial strains and plasmids used in this study are listed in Fig. S3d.

Deep genome sequencing and chromosome coverage analysis.For each strain, synchronized cells were obtained by starving in artificial seawater at 25°C for 24 h as described previously (41). These cells were then cultivated in marine broth 2216 at 25°C and harvested at the exponential phase (1 to 4 h) and the stationary phase (38 to 40 h). Genomic DNA was extracted using the PowerMax soil DNA isolation kit (Mo Bio laboratories, Inc., USA) according to the product instructions and sequenced with an Illumina HiSeq 2500 system (Illumina, USA). A 300-bp insert size library was constructed for each sample. Clean reads were mapped to replicons using Burrows-Wheeler Aligner (BWA) (65). The coverage of every base pair was calculated using the mpileup subprogram in SAMtools v.1.3.1 (66). The coverage data were further grouped in bins of 1,000 bp. For P. tunicata and P. spongiae, the potential bias in the exponential-phase coverage was removed by correction using the stationary-phase coverage. First, based on the stationary-phase data, a correction factor for each bin was calculated as the coverage of the bin divided by the average coverage of the replicon. Next, the exponential-phase coverage of each bin was corrected after division by the corresponding correction factor. Besides, outliers were removed from the coverage profiles by comparing the coverage of each bin with its 50 neighboring bins and removing bins with a coverage difference that exceeded twice the difference between the first quartile and the third quartile of coverages of neighboring bins.

Phylogenomic analysis.Orthologous genes were found by using OrthoFinder v2.2.1 (67). Single-copy genes were selected for the main chromosome and the chromid, respectively. Amino acid sequences were aligned with MAFFT v7.222 with default parameters (68). Alignments were trimmed using trimAl v1.4 with the parameter –st 0.001 (69). We also removed compositionally heterogeneous genes which may affect phylogenetic inference using P4 v1.2.0 (70) with the simulation-based test at a significance level of 0.05. Then, the concatenated alignments were used to reconstruct a maximum-likelihood tree using IQ-TREE v1.6.11 (71), with the model LG+G+F+C60 under the posterior mean site frequency (PMSF) approximation (72). This comprehensive model uses site-specific amino acid profiles based upon the C60 empirical frequency profiles (73). This model has been increasingly recommended for phylogenomic analysis of species that diverged very long ago and can also largely reduce the impact of long branch attraction (72, 74, 75). For the main chromosome tree, Algicola sagamiensis DSM 14643, a close relative of Pseudoalteromonas, was used as the outgroup. For the chromid tree, the tree was rooted at the midpoint, since the close relatives of the genus Pseudoalteromonas have no chromids and therefore cannot be used as the outgroup. Comparisons between different topologies were performed using the seven tests implemented in IQ-Tree (Fig. S2b).

Molecular dating of the divergence time of the Pseudoalteromonas.In order to estimate the time of transition from unidirectional replication to bidirectional replication, molecular dating analysis was performed using MCMCTree (76) with the relaxed-clock model and approximate likelihood calculation (77) on the basis of amino acid alignments of 25 conserved bacterial proteins used in the study (78). Outgroup species were selected based on the study (79) and four calibration points were used, which correspond to the split time of alpha- and gammaproteobacteria, the origin of cyanobacteria, the origin of the Pleurocapsales, and the origin of the Nostocales, respectively (nodes 1 to 4 in Fig. S6). The age of node 1 was set to a range from 2,620 to 2,360 million years ago (MYA) (80), and the age of node 2 was set to range from 3,000 to 2,320 MYA based on geochemical evidence for the rise of oxygen (81, 82). The age ranges of node 3 and node 4 were set to be 1,900 to 1,700 MYA and 1,900 to 1,600 MYA, respectively, based on a previous study (83). Alternatively, the age ranges of node 3 and node 4 were also set to be 2,450 to 1,750 MYA and 2,450 to 2,100 MYA, respectively, based on a previous study (84), and a similar divergence time was predicted for the emergence of bidirectional replication.

Annotation of insertion sequences on chromids.Insertion sequences were annotated using the ISsaga server with default parameters (http://issaga.biotoul.fr/issaga_index.php) (85). Integrons and cassette arrays were annotated using the IntegronFinder server (86). Prophage loci were annotated using phiSpy v2.1 with default parameters (87) and the PHASTER server (88).

Data availability.Genome sequences obtained in this study have been deposited in GenBank/DDBJ under accession numbers CP011011, CP011012, CP011025 to CP011042, CP011924, CP011925, AHBZ00000000, AHCB00000000, AHCD00000000, AHCF00000000, AQGU00000000 to AQGW00000000, AQGY00000000, AQHA00000000 to AQHC00000000, AQHE00000000 to AQHH00000000, and AQHJ00000000. The alignments, phylogenies, and scripts (89) are available at https://figshare.com/s/2145896d75319c896d58.

ACKNOWLEDGMENTS

This work was supported by the National Key R & D Program of China (2016YFA0601303, 2018YFC1406700, 2018YFC1406702), the National Natural Science Foundation of China (31630012, U1706207, 31728001, 91851205, 31770412, and 42006095), Major Scientific and Technological Innovation Project (MSTIP) of Shandong Province (2019JZZY010817), the AoShan Talents Cultivation Program Supported by Pilot National Laboratory for Marine Science and Technology (Qingdao), China (2017ASTCP-OS14), the Taishan Scholars Program of Shandong Province, China (tspd20181203), the Young Scholars Program of Shandong University (2016WLJH41), Youth Interdisciplinary Science and Innovative Research Groups of Shandong University (2020QNQT006), and the Ph.D. Research Start-Up Fund of Liaoning Province (2019-BS-009). H.L. was supported by the Hong Kong Research Grants Council Area of Excellence Scheme (AoE/M-403/16), and S.W. was supported by the impact postdoctoral fellowship scheme (IPDFS) of CUHK.

B.-B.X. and J.-C.R. performed replication direction analyses, genome sequencing, assembly, annotation, and comparative genomic analyses. B.-L.T. performed the minimal replicon assay. S.W. and B.-B.X. performed phylogenomic analyses. S.W. performed the molecular dating analyses. G.L. helped in genome sequencing and assembly. W.Z. helped in data analyses. B.-B.X., J.-C.R., and S.W. wrote the manuscript. H.L., Y.C., X.-L.C., F.L., S.L., and Q.S. did critical revision of the manuscript for important intellectual content. Y.-Z.Z. and B.-B.X. designed the study.

We declare no conflict of interest.

FOOTNOTES

    • Received 5 October 2020
    • Accepted 1 December 2020
    • Published 26 January 2021
  • Copyright © 2021 Xie et al.

This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license.

REFERENCES

  1. 1.↵
    1. Prescott DM,
    2. Kuempel PL
    . 1972. Bidirectional replication of the chromosome in Escherichia coli. Proc Natl Acad Sci U S A 69:2842–2845. doi:10.1073/pnas.69.10.2842.
    OpenUrlAbstract/FREE Full Text
  2. 2.↵
    1. Blattner FR,
    2. Plunkett G, III.,
    3. Bloch CA,
    4. Perna NT,
    5. Burland V,
    6. Riley M,
    7. Collado-Vides J,
    8. Glasner JD,
    9. Rode CK,
    10. Mayhew GF,
    11. Gregor J,
    12. Davis NW,
    13. Kirkpatrick HA,
    14. Goeden MA,
    15. Rose DJ,
    16. Mau B,
    17. Shao Y
    . 1997. The complete genome sequence of Escherichia coli K-12. Science 277:1453–1462. doi:10.1126/science.277.5331.1453.
    OpenUrlAbstract/FREE Full Text
  3. 3.↵
    1. Kunst F,
    2. Ogasawara N,
    3. Moszer I,
    4. Albertini AM,
    5. Alloni G,
    6. Azevedo V,
    7. Bertero MG,
    8. Bessieres P,
    9. Bolotin A,
    10. Borchert S,
    11. Borriss R,
    12. Boursier L,
    13. Brans A,
    14. Braun M,
    15. Brignell SC,
    16. Bron S,
    17. Brouillet S,
    18. Bruschi CV,
    19. Caldwell B,
    20. Capuano V,
    21. Carter NM,
    22. Choi SK,
    23. Cordani JJ,
    24. Connerton IF,
    25. Cummings NJ,
    26. Daniel RA,
    27. Denziot F,
    28. Devine KM,
    29. Dusterhoft A,
    30. Ehrlich SD,
    31. Emmerson PT,
    32. Entian KD,
    33. Errington J,
    34. Fabret C,
    35. Ferrari E,
    36. Foulger D,
    37. Fritz C,
    38. Fujita M,
    39. Fujita Y,
    40. Fuma S,
    41. Galizzi A,
    42. Galleron N,
    43. Ghim SY,
    44. Glaser P,
    45. Goffeau A,
    46. Golightly EJ,
    47. Grandi G,
    48. Guiseppi G,
    49. Guy BJ,
    50. Haga K
    , et al. 1997. The complete genome sequence of the gram-positive bacterium Bacillus subtilis. Nature 390:249–256. doi:10.1038/36786.
    OpenUrlCrossRefPubMedWeb of Science
  4. 4.↵
    1. Reyes-Lamothe R,
    2. Nicolas E,
    3. Sherratt DJ
    . 2012. Chromosome replication and segregation in bacteria. Annu Rev Genet 46:121–143. doi:10.1146/annurev-genet-110711-155421.
    OpenUrlCrossRefPubMedWeb of Science
  5. 5.↵
    1. Lesterlin C,
    2. Mercier R,
    3. Boccard F,
    4. Barre FX,
    5. Cornet F
    . 2005. Roles for replichores and macrodomains in segregation of the Escherichia coli chromosome. EMBO Rep 6:557–562. doi:10.1038/sj.embor.7400428.
    OpenUrlAbstract/FREE Full Text
  6. 6.↵
    1. Egan ES,
    2. Fogel MA,
    3. Waldor MK
    . 2005. Divided genomes: negotiating the cell cycle in prokaryotes with multiple chromosomes. Mol Microbiol 56:1129–1138. doi:10.1111/j.1365-2958.2005.04622.x.
    OpenUrlCrossRefPubMedWeb of Science
  7. 7.↵
    1. Jha JK,
    2. Baek JH,
    3. Venkova-Canova T,
    4. Chattoraj DK
    . 2012. Chromosome dynamics in multichromosome bacteria. Biochim Biophys Acta 1819:826–829. doi:10.1016/j.bbagrm.2012.01.012.
    OpenUrlCrossRefPubMed
  8. 8.↵
    1. Heidelberg JF,
    2. Eisen JA,
    3. Nelson WC,
    4. Clayton RA,
    5. Gwinn ML,
    6. Dodson RJ,
    7. Haft DH,
    8. Hickey EK,
    9. Peterson JD,
    10. Umayam L,
    11. Gill SR,
    12. Nelson KE,
    13. Read TD,
    14. Tettelin H,
    15. Richardson D,
    16. Ermolaeva MD,
    17. Vamathevan J,
    18. Bass S,
    19. Qin H,
    20. Dragoi I,
    21. Sellers P,
    22. McDonald L,
    23. Utterback T,
    24. Fleishmann RD,
    25. Nierman WC,
    26. White O,
    27. Salzberg SL,
    28. Smith HO,
    29. Colwell RR,
    30. Mekalanos JJ,
    31. Venter JC,
    32. Fraser CM
    . 2000. DNA sequence of both chromosomes of the cholera pathogen Vibrio cholerae. Nature 406:477–483. doi:10.1038/35020000.
    OpenUrlCrossRefPubMedWeb of Science
  9. 9.↵
    1. Harrison PW,
    2. Lower RP,
    3. Kim NK,
    4. Young JP
    . 2010. Introducing the bacterial 'chromid': not a chromosome, not a plasmid. Trends Microbiol 18:141–148. doi:10.1016/j.tim.2009.12.010.
    OpenUrlCrossRefPubMedWeb of Science
  10. 10.↵
    1. Val ME,
    2. Marbouty M,
    3. de Lemos Martins F,
    4. Kennedy SP,
    5. Kemble H,
    6. Bland MJ,
    7. Possoz C,
    8. Koszul R,
    9. Skovgaard O,
    10. Mazel D
    . 2016. A checkpoint control orchestrates the replication of the two chromosomes of Vibrio cholerae. Sci Adv 2:e1501914. doi:10.1126/sciadv.1501914.
    OpenUrlFREE Full Text
  11. 11.↵
    1. Rocha EP
    . 2004. The replication-related organization of bacterial genomes. Microbiology (Reading) 150:1609–1627. doi:10.1099/mic.0.26974-0.
    OpenUrlCrossRefPubMedWeb of Science
  12. 12.↵
    1. Dillon MM,
    2. Sung W,
    3. Lynch M,
    4. Cooper VS
    . 2018. Periodic variation of mutation rates in bacterial genomes associated with replication timing. mBio 9:e01371-18. doi:10.1128/mBio.01371-18.
    OpenUrlAbstract/FREE Full Text
  13. 13.↵
    1. Foster PL,
    2. Hanson AJ,
    3. Lee H,
    4. Popodi EM,
    5. Tang H
    . 2013. On the mutational topology of the bacterial genome. G3 (Bethesda) 3:399–407. doi:10.1534/g3.112.005355.
    OpenUrlCrossRef
  14. 14.↵
    1. Mira A,
    2. Ochman H
    . 2002. Gene location and bacterial sequence divergence. Mol Biol Evol 19:1350–1358. doi:10.1093/oxfordjournals.molbev.a004196.
    OpenUrlCrossRefPubMedWeb of Science
  15. 15.↵
    1. Cooper VS,
    2. Vohr SH,
    3. Wrocklage SC,
    4. Hatcher PJ
    . 2010. Why genes evolve faster on secondary chromosomes in bacteria. PLoS Comput Biol 6:e1000732. doi:10.1371/journal.pcbi.1000732.
    OpenUrlCrossRefPubMed
  16. 16.↵
    1. Lobry JR
    . 1996. Asymmetric substitution patterns in the two DNA strands of bacteria. Mol Biol Evol 13:660–665. doi:10.1093/oxfordjournals.molbev.a025626.
    OpenUrlCrossRefPubMedWeb of Science
  17. 17.↵
    1. Kono N,
    2. Tomita M,
    3. Arakawa K
    . 2018. Accelerated laboratory evolution reveals the influence of replication on the GC skew in Escherichia coli. Genome Biol Evol 10:3110–3117. doi:10.1093/gbe/evy237.
    OpenUrlCrossRef
  18. 18.↵
    1. Grigoriev A
    . 1998. Analyzing genomes with cumulative skew diagrams. Nucleic Acids Res 26:2286–2290. doi:10.1093/nar/26.10.2286.
    OpenUrlCrossRefPubMedWeb of Science
  19. 19.↵
    1. Medigue C,
    2. Krin E,
    3. Pascal G,
    4. Barbe V,
    5. Bernsel A,
    6. Bertin PN,
    7. Cheung F,
    8. Cruveiller S,
    9. D'Amico S,
    10. Duilio A,
    11. Fang G,
    12. Feller G,
    13. Ho C,
    14. Mangenot S,
    15. Marino G,
    16. Nilsson J,
    17. Parrilli E,
    18. Rocha EP,
    19. Rouy Z,
    20. Sekowska A,
    21. Tutino ML,
    22. Vallenet D,
    23. von Heijne G,
    24. Danchin A
    . 2005. Coping with cold: the genome of the versatile marine Antarctica bacterium Pseudoalteromonas haloplanktis TAC125. Genome Res 15:1325–1335. doi:10.1101/gr.4126905.
    OpenUrlAbstract/FREE Full Text
  20. 20.↵
    1. Bowman JP
    . 2007. Bioactive compound synthetic capacity and ecological significance of marine bacterial genus pseudoalteromonas. Mar Drugs 5:220–241. doi:10.3390/md504220.
    OpenUrlCrossRefPubMedWeb of Science
  21. 21.↵
    1. Wietz M,
    2. Gram L,
    3. Jørgensen B,
    4. Schramm A
    . 2010. Latitudinal patterns in the abundance of major marine bacterioplankton groups. Aquat Microb Ecol 61:179–189. doi:10.3354/ame01443.
    OpenUrlCrossRef
  22. 22.↵
    1. Techkarnjanaruk S,
    2. Goodman AE
    . 1999. Multiple genes involved in chitin degradation from the marine bacterium Pseudoalteromonas sp. strain S91. Microbiology 145:925–934. doi:10.1099/13500872-145-4-925.
    OpenUrlCrossRefPubMed
  23. 23.↵
    1. Qin QL,
    2. Li Y,
    3. Zhang YJ,
    4. Zhou ZM,
    5. Zhang WX,
    6. Chen XL,
    7. Zhang XY,
    8. Zhou BC,
    9. Wang L,
    10. Zhang YZ
    . 2011. Comparative genomics reveals a deep-sea sediment-adapted life style of Pseudoalteromonas sp. SM9913. ISME J 5:274–284. doi:10.1038/ismej.2010.103.
    OpenUrlCrossRefPubMed
  24. 24.↵
    1. Chen XL,
    2. Dong S,
    3. Xu F,
    4. Dong F,
    5. Li PY,
    6. Zhang XY,
    7. Zhou BC,
    8. Zhang YZ,
    9. Xie BB
    . 2016. Characterization of a new cold-adapted and salt-activated polysaccharide lyase family 7 alginate lyase from Pseudoalteromonas sp. SM0524. Front Microbiol 7:1120. doi:10.3389/fmicb.2016.01120.
    OpenUrlCrossRefPubMed
  25. 25.↵
    1. Zhao GY,
    2. Chen XL,
    3. Zhao HL,
    4. Xie BB,
    5. Zhou BC,
    6. Zhang YZ
    . 2008. Hydrolysis of insoluble collagen by deseasin MCP-01 from deep-sea Pseudoalteromonas sp. SM9913: collagenolytic characters, collagen-binding ability of C-terminal polycystic kidney disease domain, and implication for its novel role in deep-sea sedimentary particulate organic nitrogen degradation. J Biol Chem 283:36100–36107. doi:10.1074/jbc.M804438200.
    OpenUrlAbstract/FREE Full Text
  26. 26.↵
    1. Zhao HL,
    2. Chen XL,
    3. Xie BB,
    4. Zhou MY,
    5. Gao X,
    6. Zhang XY,
    7. Zhou BC,
    8. Weiss AS,
    9. Zhang YZ
    . 2012. Elastolytic mechanism of a novel M23 metalloprotease pseudoalterin from deep-sea Pseudoalteromonas sp. CF6-2: cleaving not only glycyl bonds in the hydrophobic regions but also peptide bonds in the hydrophilic regions involved in cross-linking. J Biol Chem 287:39710–39720. doi:10.1074/jbc.M112.405076.
    OpenUrlAbstract/FREE Full Text
  27. 27.↵
    1. Tang BL,
    2. Yang J,
    3. Chen XL,
    4. Wang P,
    5. Zhao HL,
    6. Su HN,
    7. Li CY,
    8. Yu Y,
    9. Zhong S,
    10. Wang L,
    11. Lidbury I,
    12. Ding H,
    13. Wang M,
    14. McMinn A,
    15. Zhang XY,
    16. Chen Y,
    17. Zhang YZ
    . 2020. A predator-prey interaction between a marine Pseudoalteromonas sp. and Gram-positive bacteria. Nat Commun 11:285. doi:10.1038/s41467-019-14133-x.
    OpenUrlCrossRef
  28. 28.↵
    1. Gauthier G,
    2. Gauthier M,
    3. Christen R
    . 1995. Phylogenetic analysis of the genera Alteromonas, Shewanella, and Moritella using genes coding for small-subunit rRNA sequences and division of the genus Alteromonas into two genera, Alteromonas (emended) and Pseudoalteromonas gen. nov., and proposal of twelve new species combinations. Int J Syst Bacteriol 45:755–761. doi:10.1099/00207713-45-4-755.
    OpenUrlCrossRefPubMed
  29. 29.↵
    1. Giovannoni SJ,
    2. Tripp HJ,
    3. Givan S,
    4. Podar M,
    5. Vergin KL,
    6. Baptista D,
    7. Bibbs L,
    8. Eads J,
    9. Richardson TH,
    10. Noordewier M,
    11. Rappe MS,
    12. Short JM,
    13. Carrington JC,
    14. Mathur EJ
    . 2005. Genome streamlining in a cosmopolitan oceanic bacterium. Science 309:1242–1245. doi:10.1126/science.1114057.
    OpenUrlAbstract/FREE Full Text
  30. 30.↵
    1. Thomas T,
    2. Evans FF,
    3. Schleheck D,
    4. Mai-Prochnow A,
    5. Burke C,
    6. Penesyan A,
    7. Dalisay DS,
    8. Stelzer-Braid S,
    9. Saunders N,
    10. Johnson J,
    11. Ferriera S,
    12. Kjelleberg S,
    13. Egan S
    . 2008. Analysis of the Pseudoalteromonas tunicata genome reveals properties of a surface-associated life style in the marine environment. PLoS One 3:e3252. doi:10.1371/journal.pone.0003252.
    OpenUrlCrossRefPubMed
  31. 31.↵
    1. Galli E,
    2. Ferat JL,
    3. Desfontaines JM,
    4. Val ME,
    5. Skovgaard O,
    6. Barre FX,
    7. Possoz C
    . 2019. Replication termination without a replication fork trap. Sci Rep 9:8315. doi:10.1038/s41598-019-43795-2.
    OpenUrlCrossRef
  32. 32.↵
    1. Duigou S,
    2. Knudsen KG,
    3. Skovgaard O,
    4. Egan ES,
    5. Lobner-Olesen A,
    6. Waldor MK
    . 2006. Independent control of replication initiation of the two Vibrio cholerae chromosomes by DnaA and RctB. J Bacteriol 188:6419–6424. doi:10.1128/JB.00565-06.
    OpenUrlAbstract/FREE Full Text
  33. 33.↵
    1. Leslie NR,
    2. Sherratt DJ
    . 1995. Site-specific recombination in the replication terminus region of Escherichia coli: functional replacement of dif. EMBO J 14:1561–1570. doi:10.1002/j.1460-2075.1995.tb07142.x.
    OpenUrlCrossRefPubMedWeb of Science
  34. 34.↵
    1. Bhowmik BK,
    2. Clevenger AL,
    3. Zhao H,
    4. Rybenkov VV
    . 2018. Segregation but not replication of the Pseudomonas aeruginosa chromosome terminates at Dif. mBio 9:e01088-18. doi:10.1128/mBio.01088-18.
    OpenUrlCrossRef
  35. 35.↵
    1. Simmons LA,
    2. Breier AM,
    3. Cozzarelli NR,
    4. Kaguni JM
    . 2004. Hyperinitiation of DNA replication in Escherichia coli leads to replication fork collapse and inviability. Mol Microbiol 51:349–358. doi:10.1046/j.1365-2958.2003.03842.x.
    OpenUrlCrossRefPubMedWeb of Science
  36. 36.↵
    1. Skovgaard O,
    2. Bak M,
    3. Lobner-Olesen A,
    4. Tommerup N
    . 2011. Genome-wide detection of chromosomal rearrangements, indels, and mutations in circular chromosomes by short read sequencing. Genome Res 21:1388–1393. doi:10.1101/gr.117416.110.
    OpenUrlAbstract/FREE Full Text
  37. 37.↵
    1. Lundgren M,
    2. Andersson A,
    3. Chen L,
    4. Nilsson P,
    5. Bernander R
    . 2004. Three replication origins in Sulfolobus species: synchronous initiation of chromosome replication and asynchronous termination. Proc Natl Acad Sci U S A 101:7046–7051. doi:10.1073/pnas.0400656101.
    OpenUrlAbstract/FREE Full Text
  38. 38.↵
    1. Du WL,
    2. Dubarry N,
    3. Passot FM,
    4. Kamgoue A,
    5. Murray H,
    6. Lane D,
    7. Pasta F
    . 2016. Orderly replication and segregation of the four replicons of Burkholderia cenocepacia J2315. PLoS Genet 12:e1006172. doi:10.1371/journal.pgen.1006172.
    OpenUrlCrossRef
  39. 39.↵
    1. Tang BC,
    2. Rong JC,
    3. Dang YR,
    4. Xie BB,
    5. Chen XL,
    6. Zhang XY
    . 2018. Complete genomic sequence of Pseudoalteromonas sp. SAO4-4, a protease-producing bacterium isolated from seawater of the Atlantic Ocean. Genome Announc 6:e00284-18. doi:10.1128/genomeA.00284-18.
    OpenUrlCrossRef
  40. 40.↵
    1. Beurmann S,
    2. Ushijima B,
    3. Videau P,
    4. Svoboda CM,
    5. Smith AM,
    6. Rivers OS,
    7. Aeby GS,
    8. Callahan SM
    . 2017. Pseudoalteromonas piratica strain OCN003 is a coral pathogen that causes a switch from chronic to acute Montipora white syndrome in Montipora capitata. PLoS One 12:e0188319. doi:10.1371/journal.pone.0188319.
    OpenUrlCrossRef
  41. 41.↵
    1. Rasmussen T,
    2. Jensen RB,
    3. Skovgaard O
    . 2007. The two chromosomes of Vibrio cholerae are initiated at different time points in the cell cycle. EMBO J 26:3124–3131. doi:10.1038/sj.emboj.7601747.
    OpenUrlCrossRefPubMedWeb of Science
  42. 42.↵
    1. Kemter FS,
    2. Messerschmidt SJ,
    3. Schallopp N,
    4. Sobetzko P,
    5. Lang E,
    6. Bunk B,
    7. Sproer C,
    8. Teschler JK,
    9. Yildiz FH,
    10. Overmann J,
    11. Waldminghaus T
    . 2018. Synchronous termination of replication of the two chromosomes is an evolutionary selected feature in Vibrionaceae. PLoS Genet 14:e1007251. doi:10.1371/journal.pgen.1007251.
    OpenUrlCrossRef
  43. 43.↵
    1. del Solar G,
    2. Giraldo R,
    3. Ruiz-Echevarria MJ,
    4. Espinosa M,
    5. Diaz-Orejas R
    . 1998. Replication and control of circular bacterial plasmids. Microbiol Mol Biol Rev 62:434–464. doi:10.1128/MMBR.62.2.434-464.1998.
    OpenUrlAbstract/FREE Full Text
  44. 44.↵
    1. Duggin IG,
    2. Bell SD
    . 2009. Termination structures in the Escherichia coli chromosome replication fork trap. J Mol Biol 387:532–539. doi:10.1016/j.jmb.2009.02.027.
    OpenUrlCrossRefPubMed
  45. 45.↵
    1. Das B,
    2. Martinez E,
    3. Midonet C,
    4. Barre FX
    . 2013. Integrative mobile elements exploiting Xer recombination. Trends Microbiol 21:23–30. doi:10.1016/j.tim.2012.10.003.
    OpenUrlCrossRefPubMedWeb of Science
  46. 46.↵
    1. Val ME,
    2. Kennedy SP,
    3. El Karoui M,
    4. Bonne L,
    5. Chevalier F,
    6. Barre FX
    . 2008. FtsK-dependent dimer resolution on multiple chromosomes in the pathogen Vibrio cholerae. PLoS Genet 4:e1000201. doi:10.1371/journal.pgen.1000201.
    OpenUrlCrossRefPubMed
  47. 47.↵
    1. Venkova-Canova T,
    2. Chattoraj DK
    . 2011. Transition from a plasmid to a chromosomal mode of replication entails additional regulators. Proc Natl Acad Sci U S A 108:6199–6204. doi:10.1073/pnas.1013244108.
    OpenUrlAbstract/FREE Full Text
  48. 48.↵
    1. Baek JH,
    2. Chattoraj DK
    . 2014. Chromosome I controls chromosome II replication in Vibrio cholerae. PLoS Genet 10:e1004184. doi:10.1371/journal.pgen.1004184.
    OpenUrlCrossRefPubMed
  49. 49.↵
    1. Tillier ER,
    2. Collins RA
    . 2000. Genome rearrangement by replication-directed translocation. Nat Genet 26:195–197. doi:10.1038/79918.
    OpenUrlCrossRefPubMedWeb of Science
  50. 50.↵
    1. Rocha EP
    . 2008. The organization of the bacterial genome. Annu Rev Genet 42:211–233. doi:10.1146/annurev.genet.42.110807.091653.
    OpenUrlCrossRefPubMedWeb of Science
  51. 51.↵
    1. Couturier E,
    2. Rocha EP
    . 2006. Replication-associated gene dosage effects shape the genomes of fast-growing bacteria but only for transcription and translation genes. Mol Microbiol 59:1506–1518. doi:10.1111/j.1365-2958.2006.05046.x.
    OpenUrlCrossRefPubMedWeb of Science
  52. 52.↵
    1. Mackiewicz P,
    2. Zakrzewska-Czerwinska J,
    3. Zawilak A,
    4. Dudek MR,
    5. Cebrat S
    . 2004. Where does bacterial replication start? Rules for predicting the oriC region. Nucleic Acids Res 32:3781–3791. doi:10.1093/nar/gkh699.
    OpenUrlCrossRefPubMedWeb of Science
  53. 53.↵
    1. Dillon MM,
    2. Sung W,
    3. Sebra R,
    4. Lynch M,
    5. Cooper VS
    . 2017. Genome-wide biases in the rate and molecular spectrum of spontaneous mutations in Vibrio cholerae and Vibrio fischeri. Mol Biol Evol 34:93–109. doi:10.1093/molbev/msw224.
    OpenUrlCrossRefPubMed
  54. 54.↵
    1. Long H,
    2. Sung W,
    3. Miller SF,
    4. Ackerman MS,
    5. Doak TG,
    6. Lynch M
    . 2014. Mutation rate, spectrum, topology, and context-dependency in the DNA mismatch repair-deficient Pseudomonas fluorescens ATCC948. Genome Biol Evol 7:262–271. doi:10.1093/gbe/evu284.
    OpenUrlCrossRefPubMed
  55. 55.↵
    1. Sung W,
    2. Ackerman MS,
    3. Gout JF,
    4. Miller SF,
    5. Williams E,
    6. Foster PL,
    7. Lynch M
    . 2015. Asymmetric context-dependent mutation patterns revealed through mutation-accumulation experiments. Mol Biol Evol 32:1672–1683. doi:10.1093/molbev/msv055.
    OpenUrlCrossRefPubMed
  56. 56.↵
    1. Dettman JR,
    2. Sztepanacz JL,
    3. Kassen R
    . 2016. The properties of spontaneous mutations in the opportunistic pathogen Pseudomonas aeruginosa. BMC Genomics 17:27. doi:10.1186/s12864-015-2244-3.
    OpenUrlCrossRef
  57. 57.↵
    1. Li RQ,
    2. Zhu HM,
    3. Ruan J,
    4. Qian WB,
    5. Fang XD,
    6. Shi ZB,
    7. Li YR,
    8. Li ST,
    9. Shan G,
    10. Kristiansen K,
    11. Li SG,
    12. Yang HM,
    13. Wang J,
    14. Wang J
    . 2010. De novo assembly of human genomes with massively parallel short read sequencing. Genome Res 20:265–272. doi:10.1101/gr.097261.109.
    OpenUrlAbstract/FREE Full Text
  58. 58.↵
    1. Kurtz S,
    2. Phillippy A,
    3. Delcher AL,
    4. Smoot M,
    5. Shumway M,
    6. Antonescu C,
    7. Salzberg SL
    . 2004. Versatile and open software for comparing large genomes. Genome Biol 5:R12. doi:10.1186/gb-2004-5-2-r12.
    OpenUrlCrossRefPubMed
  59. 59.↵
    1. Delcher AL,
    2. Bratke KA,
    3. Powers EC,
    4. Salzberg SL
    . 2007. Identifying bacterial genes and endosymbiont DNA with Glimmer. Bioinformatics 23:673–679. doi:10.1093/bioinformatics/btm009.
    OpenUrlCrossRefPubMedWeb of Science
  60. 60.↵
    1. Tatusov RL,
    2. Koonin EV,
    3. Lipman DJ
    . 1997. A genomic perspective on protein families. Science 278:631–637. doi:10.1126/science.278.5338.631.
    OpenUrlAbstract/FREE Full Text
  61. 61.↵
    1. Carver T,
    2. Thomson N,
    3. Bleasby A,
    4. Berriman M,
    5. Parkhill J
    . 2009. DNAPlotter: circular and linear interactive genome visualization. Bioinformatics 25:119–120. doi:10.1093/bioinformatics/btn578.
    OpenUrlCrossRefPubMedWeb of Science
  62. 62.↵
    1. Kono N,
    2. Arakawa K,
    3. Tomita M
    . 2011. Comprehensive prediction of chromosome dimer resolution sites in bacterial genomes. BMC Genomics 12:19. doi:10.1186/1471-2164-12-19.
    OpenUrlCrossRefPubMed
  63. 63.↵
    1. Yu ZC,
    2. Zhao DL,
    3. Ran LY,
    4. Mi ZH,
    5. Wu ZY,
    6. Pang X,
    7. Zhang XY,
    8. Su HN,
    9. Shi M,
    10. Song XY,
    11. Xie BB,
    12. Qin QL,
    13. Zhou BC,
    14. Chen XL,
    15. Zhang YZ
    . 2014. Development of a genetic system for the deep-sea psychrophilic bacterium Pseudoalteromonas sp. SM9913. Microb Cell Fact 13:13. doi:10.1186/1475-2859-13-13.
    OpenUrlCrossRef
  64. 64.↵
    1. Wang P,
    2. Yu Z,
    3. Li B,
    4. Cai X,
    5. Zeng Z,
    6. Chen X,
    7. Wang X
    . 2015. Development of an efficient conjugation-based genetic manipulation system for Pseudoalteromonas. Microb Cell Fact 14:11. doi:10.1186/s12934-015-0194-8.
    OpenUrlCrossRefPubMed
  65. 65.↵
    1. Li H,
    2. Durbin R
    . 2010. Fast and accurate long-read alignment with Burrows-Wheeler transform. Bioinformatics 26:589–595. doi:10.1093/bioinformatics/btp698.
    OpenUrlCrossRefPubMedWeb of Science
  66. 66.↵
    1. Li H
    . 2011. A statistical framework for SNP calling, mutation discovery, association mapping and population genetical parameter estimation from sequencing data. Bioinformatics 27:2987–2993. doi:10.1093/bioinformatics/btr509.
    OpenUrlCrossRefPubMedWeb of Science
  67. 67.↵
    1. Emms DM,
    2. Kelly S
    . 2019. OrthoFinder: phylogenetic orthology inference for comparative genomics. Genome Biol 20:238. doi:10.1186/s13059-019-1832-y.
    OpenUrlCrossRef
  68. 68.↵
    1. Katoh K,
    2. Standley DM
    . 2013. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol Biol Evol 30:772–780. doi:10.1093/molbev/mst010.
    OpenUrlCrossRefPubMedWeb of Science
  69. 69.↵
    1. Capella-Gutierrez S,
    2. Silla-Martinez JM,
    3. Gabaldon T
    . 2009. trimAl: a tool for automated alignment trimming in large-scale phylogenetic analyses. Bioinformatics 25:1972–1973. doi:10.1093/bioinformatics/btp348.
    OpenUrlCrossRefPubMedWeb of Science
  70. 70.↵
    1. Foster PG
    . 2004. Modeling compositional heterogeneity. Syst Biol 53:485–495. doi:10.1080/10635150490445779.
    OpenUrlCrossRefPubMedWeb of Science
  71. 71.↵
    1. Nguyen LT,
    2. Schmidt HA,
    3. von Haeseler A,
    4. Minh BQ
    . 2015. IQ-TREE: a fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. Mol Biol Evol 32:268–274. doi:10.1093/molbev/msu300.
    OpenUrlCrossRefPubMed
  72. 72.↵
    1. Wang HC,
    2. Minh BQ,
    3. Susko E,
    4. Roger AJ
    . 2018. Modeling site heterogeneity with posterior mean site frequency profiles accelerates accurate phylogenomic estimation. Syst Biol 67:216–235. doi:10.1093/sysbio/syx068.
    OpenUrlCrossRef
  73. 73.↵
    1. Quang Le S,
    2. Gascuel O,
    3. Lartillot N
    . 2008. Empirical profile mixture models for phylogenetic reconstruction. Bioinformatics 24:2317–2323. doi:10.1093/bioinformatics/btn445.
    OpenUrlCrossRefPubMedWeb of Science
  74. 74.↵
    1. Lartillot N,
    2. Rodrigue N,
    3. Stubbs D,
    4. Richer J
    . 2013. PhyloBayes MPI: phylogenetic reconstruction with infinite mixtures of profiles in a parallel environment. Syst Biol 62:611–615. doi:10.1093/sysbio/syt022.
    OpenUrlCrossRefPubMed
  75. 75.↵
    1. Wang HC,
    2. Susko E,
    3. Roger AJ
    . 2019. The relative importance of modeling site pattern heterogeneity versus partition-wise heterotachy in phylogenomic inference. Syst Biol 68:1003–1019. doi:10.1093/sysbio/syz021.
    OpenUrlCrossRef
  76. 76.↵
    1. dosReis M,
    2. Yang Z
    . 2019. Bayesian molecular clock dating using genome-scale datasets. Methods Mol Biol 1910:309–330. doi:10.1007/978-1-4939-9074-0_10.
    OpenUrlCrossRef
  77. 77.↵
    1. dos Reis M,
    2. Yang Z
    . 2011. Approximate likelihood calculation on a phylogeny for Bayesian estimation of divergence times. Mol Biol Evol 28:2161–2172. doi:10.1093/molbev/msr045.
    OpenUrlCrossRefPubMedWeb of Science
  78. 78.↵
    1. Battistuzzi FU,
    2. Hedges SB
    . 2009. A major clade of prokaryotes with ancient adaptations to life on land. Mol Biol Evol 26:335–343. doi:10.1093/molbev/msn247.
    OpenUrlCrossRefPubMedWeb of Science
  79. 79.↵
    1. Luo H,
    2. Csűros M,
    3. Hughes AL,
    4. Moran MA
    . 2013. Evolution of divergent life history strategies in marine alphaproteobacteria. mBio 4:e00373-13. doi:10.1128/mBio.00373-13.
    OpenUrlAbstract/FREE Full Text
  80. 80.↵
    1. Kumar S,
    2. Stecher G,
    3. Suleski M,
    4. Hedges SB
    . 2017. TimeTree: a resource for timelines, timetrees, and divergence times. Mol Biol Evol 34:1812–1819. doi:10.1093/molbev/msx116.
    OpenUrlCrossRefPubMed
  81. 81.↵
    1. Crowe SA,
    2. Dossing LN,
    3. Beukes NJ,
    4. Bau M,
    5. Kruger SJ,
    6. Frei R,
    7. Canfield DE
    . 2013. Atmospheric oxygenation three billion years ago. Nature 501:535–538. doi:10.1038/nature12426.
    OpenUrlCrossRefPubMedWeb of Science
  82. 82.↵
    1. Brocks JJ,
    2. Buick R,
    3. Summons RE,
    4. Logan GA
    . 2003. A reconstruction of Archean biological diversity based on molecular fossils from the 2.78 to 2.45 billion-year-old Mount Bruce supergroup, Hamersley Basin, Western Australia. Geochim Cosmochim Acta 67:4321–4335. doi:10.1016/S0016-7037(03)00209-6.
    OpenUrlCrossRefWeb of Science
  83. 83.↵
    1. Sanchez-Baracaldo P,
    2. Raven JA,
    3. Pisani D,
    4. Knoll AH
    . 2017. Early photosynthetic eukaryotes inhabited low-salinity habitats. Proc Natl Acad Sci U S A 114:E7737–E7745. doi:10.1073/pnas.1620089114.
    OpenUrlAbstract/FREE Full Text
  84. 84.↵
    1. Sanchez-Baracaldo P,
    2. Ridgwell A,
    3. Raven JA
    . 2014. A neoproterozoic transition in the marine nitrogen cycle. Curr Biol 24:652–657. doi:10.1016/j.cub.2014.01.041.
    OpenUrlCrossRefPubMed
  85. 85.↵
    1. Varani AM,
    2. Siguier P,
    3. Gourbeyre E,
    4. Charneau V,
    5. Chandler M
    . 2011. ISsaga is an ensemble of web-based methods for high throughput identification and semi-automatic annotation of insertion sequences in prokaryotic genomes. Genome Biol 12:R30. doi:10.1186/gb-2011-12-3-r30.
    OpenUrlCrossRefPubMed
  86. 86.↵
    1. Cury J,
    2. Jove T,
    3. Touchon M,
    4. Neron B,
    5. Rocha EP
    . 2016. Identification and analysis of integrons and cassette arrays in bacterial genomes. Nucleic Acids Res 44:4539–4550. doi:10.1093/nar/gkw319.
    OpenUrlCrossRefPubMed
  87. 87.↵
    1. Akhter S,
    2. Aziz RK,
    3. Edwards RA
    . 2012. PhiSpy: a novel algorithm for finding prophages in bacterial genomes that combines similarity- and composition-based strategies. Nucleic Acids Res 40:e126. doi:10.1093/nar/gks406.
    OpenUrlCrossRefPubMed
  88. 88.↵
    1. Arndt D,
    2. Grant JR,
    3. Marcu A,
    4. Sajed T,
    5. Pon A,
    6. Liang Y,
    7. Wishart DS
    . 2016. PHASTER: a better, faster version of the PHAST phage search tool. Nucleic Acids Res 44:W16–W21. doi:10.1093/nar/gkw387.
    OpenUrlCrossRefPubMed
  89. 89.↵
    1. Goto N,
    2. Prins P,
    3. Nakao M,
    4. Bonnal R,
    5. Aerts J,
    6. Katayama T
    . 2010. BioRuby: bioinformatics software for the Ruby programming language. Bioinformatics 26:2617–2619. doi:10.1093/bioinformatics/btq475.
    OpenUrlCrossRefPubMedWeb of Science
PreviousNext
Back to top
Download PDF
Citation Tools
Evolutionary Trajectory of the Replication Mode of Bacterial Replicons
Bin-Bin Xie, Jin-Cheng Rong, Bai-Lu Tang, Sishuo Wang, Guiming Liu, Qi-Long Qin, Xi-Ying Zhang, Weipeng Zhang, Qunxin She, Yin Chen, Fuchuan Li, Shengying Li, Xiu-Lan Chen, Haiwei Luo, Yu-Zhong Zhang
mBio Jan 2021, 12 (1) e02745-20; DOI: 10.1128/mBio.02745-20

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Print

Alerts
Sign In to Email Alerts with your Email Address
Email

Thank you for sharing this mBio article.

NOTE: We request your email address only to inform the recipient that it was you who recommended this article, and that it is not junk mail. We do not retain these email addresses.

Enter multiple addresses on separate lines or separate them with commas.
Evolutionary Trajectory of the Replication Mode of Bacterial Replicons
(Your Name) has forwarded a page to you from mBio
(Your Name) thought you would be interested in this article in mBio.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Share
Evolutionary Trajectory of the Replication Mode of Bacterial Replicons
Bin-Bin Xie, Jin-Cheng Rong, Bai-Lu Tang, Sishuo Wang, Guiming Liu, Qi-Long Qin, Xi-Ying Zhang, Weipeng Zhang, Qunxin She, Yin Chen, Fuchuan Li, Shengying Li, Xiu-Lan Chen, Haiwei Luo, Yu-Zhong Zhang
mBio Jan 2021, 12 (1) e02745-20; DOI: 10.1128/mBio.02745-20
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Top
  • Article
    • ABSTRACT
    • INTRODUCTION
    • RESULTS
    • DISCUSSION
    • MATERIALS AND METHODS
    • ACKNOWLEDGMENTS
    • FOOTNOTES
    • REFERENCES
  • Figures & Data
  • Info & Metrics
  • PDF

KEYWORDS

chromid
chromosome replication
unidirectional replication
chromosome evolution
Pseudoalteromonas

Related Articles

Cited By...

About

  • About mBio
  • Editor in Chief
  • Board of Editors
  • AAM Fellows
  • Policies
  • For Reviewers
  • For the Media
  • For Librarians
  • For Advertisers
  • Alerts
  • RSS
  • FAQ
  • Permissions
  • Journal Announcements

Authors

  • ASM Author Center
  • Submit a Manuscript
  • Author Warranty
  • Article Types
  • Ethics
  • Contact Us

Follow #mBio

@ASMicrobiology

       

ASM Journals

ASM journals are the most prominent publications in the field, delivering up-to-date and authoritative coverage of both basic and clinical microbiology.

About ASM | Contact Us | Press Room

 

ASM is a member of

Scientific Society Publisher Alliance

 

American Society for Microbiology
1752 N St. NW
Washington, DC 20036
Phone: (202) 737-3600

Copyright © 2021 American Society for Microbiology | Privacy Policy | Website feedback

Online ISSN: 2150-7511