Identification of replication origins in prokaryotic genomes

High Nucleotide Skew Palindromic DNA Sequences Function as Potential Replication Origins due to their Unzipping Propensity

Locations of DNA replication initiation in prokaryotes, called "origins of replication", are well-characterized. However, a mechanistic understanding of the sequence dependence of the local unzipping of double-stranded DNA, the first step towards replication initiation, is lacking. Here, utilizing a Markov chain model that was created to address the directional nature of DNA unzipping and replication, we model the sequence dependence of local melting of double-stranded linear DNA segments. We show that generalized palindromic sequences with high nucleotide skews have a low kinetic barrier for local melting near melting temperatures. This allows for such sequences to function as potential replication origins. We support our claim with evidence for high-skew palindromic sequences within the replication origins of mitochondrial DNA, bacteria, archaea and plasmids.

OUP accepted manuscript

Conflicts between transcription and replication machinery are a potent source of replication stress and genome instability; however, no technique currently exists to identify endogenous genomic locations prone to transcription–replication interactions. Here, we report a novel method to identify genomic loci prone to transcription–replication interactions termed transcription–replication immunoprecipitation on nascent DNA sequencing, TRIPn-Seq. TRIPn-Seq employs the sequential immunoprecipitation of RNA polymerase 2 phosphorylated at serine 5 (RNAP2s5) followed by enrichment of nascent DNA previously labeled with bromodeoxyuridine. Using TRIPn-Seq, we mapped 1009 unique transcription–replication interactions (TRIs) in mouse primary B cells characterized by a bimodal pattern of RNAP2s5, bidirectional transcription, an enrichment of RNA:DNA hybrids, and a high probability of forming G-quadruplexes. TRIs are highly enriched at transcription start sites and map to early replicating regions. TRIs exhibit enhanced Replication Protein A association and TRI-associated genes exhibit higher replication fork termination than control transcription start sites, two marks of replication stress. TRIs colocalize with double-strand DNA breaks, are enriched for deletions, and accumulate mutations in tumors. We propose that replication stress at TRIs induces mutations potentially contributing to age-related disease, as well as tumor formation and development.

Virus-associated organosulfur metabolism in human and environmental systems

Viruses influence the fate of nutrients and human health by killing microorganisms and altering metabolic processes. Organosulfur metabolism and biologically derived hydrogen sulfide play dynamic roles in manifestation of diseases, infrastructure degradation, and essential biological processes. Although microbial organosulfur metabolism is well studied, the role of viruses in organosulfur metabolism is unknown. Here, we report the discovery of 39 gene families involved in organosulfur metabolism encoded by 3,749 viruses from diverse ecosystems, including human microbiomes. The viruses infect organisms from all three domains of life. Six gene families encode for enzymes that degrade organosulfur compounds into sulfide, whereas others manipulate organosulfur compounds and may influence sulfide production. We show that viral metabolic genes encode key enzymatic domains, are translated into protein, and are maintained after recombination, and sulfide provides a fitness advantage to viruses. Our results reveal viruses as drivers of organosulfur metabolism with important implications for human and environmental health.

Recombineering and MAGE

Recombination-mediated genetic engineering, also known as recombineering, is the genomic incorporation of homologous single-stranded or double-stranded DNA into bacterial genomes. Recombineering and its derivative methods have radically improved genome engineering capabilities, perhaps none more so than multiplex automated genome engineering (MAGE). MAGE is representative of a set of highly multiplexed single-stranded DNA-mediated technologies. First described in Escherichia coli, both MAGE and recombineering are being rapidly translated into diverse prokaryotes and even into eukaryotic cells. Together, this modern set of tools offers the promise of radically improving the scope and throughput of experimental biology by providing powerful new methods to ease the genetic manipulation of model and non-model organisms. In this Primer, we describe recombineering and MAGE, their optimal use, their diverse applications and methods for pairing them with other genetic editing tools. We then look forward to the future of genetic engineering.

Probabilistic model based on circular statistics for quantifying coverage depth dynamics originating from DNA replication

BACKGROUND

With the development of DNA sequencing technology, static omics profiling in microbial communities, such as taxonomic and functional gene composition determination, has become possible. Additionally, the recently proposed in situ growth rate estimation method allows the applicable range of current comparative metagenomics to be extended to dynamic profiling. However, with this method, the applicable target range is presently limited. Furthermore, the characteristics of coverage depth during replication have not been sufficiently investigated.

RESULTS

We developed a probabilistic model that mimics coverage depth dynamics. This statistical model explains the bias that occurs in the coverage depth due to DNA replication and errors that arise from coverage depth observation. Although our method requires a complete genome sequence, it involves a stable to low coverage depth (>0.01×). We also evaluated the estimation using real whole-genome sequence datasets and reproduced the growth dynamics observed in previous studies. By utilizing a circular distribution in the model, our method facilitates the quantification of unmeasured coverage depth features, including peakedness, skewness, and degree of density, around the replication origin. When we applied the model to time-series culture samples, the skewness parameter, which indicates the asymmetry, was stable over time; however, the peakedness and degree of density parameters, which indicate the concentration level at the replication origin, changed dynamically. Furthermore, we demonstrated the activity measurement of multiple replication origins in a single chromosome.

CONCLUSIONS

We devised a novel framework for quantifying coverage depth dynamics. Our study is expected to serve as a basis for replication activity estimation from a broader perspective using the statistical model.

Two-step chromosome segregation in the stalked budding bacterium Hyphomonas neptunium

Chromosome segregation typically occurs after replication has finished in eukaryotes but during replication in bacteria. Here, we show that the alphaproteobacterium Hyphomonas neptunium, which proliferates by bud formation at the tip of a stalk-like cellular extension, segregates its chromosomes in a unique two-step process. First, the two sister origin regions are targeted to opposite poles of the mother cell, driven by the ParABS partitioning system. Subsequently, once the bulk of chromosomal DNA has been replicated and the bud exceeds a certain threshold size, the cell initiates a second segregation step during which it transfers the stalk-proximal origin region through the stalk into the nascent bud compartment. Thus, while chromosome replication and segregation usually proceed concurrently in bacteria, the two processes are largely uncoupled in H. neptunium, reminiscent of eukaryotic mitosis. These results indicate that stalked budding bacteria have evolved specific mechanisms to adjust chromosome segregation to their unusual life cycle.

Reconsidering plasmid maintenance factors for computational plasmid design

Plasmids are genetic parasites of microorganisms. The genomes of naturally occurring plasmids are expected to be polished via natural selection to achieve long-term persistence in the microbial cell population. However, plasmid genomes are extremely diverse, and the rules governing plasmid genomes are not fully understood. Therefore, computationally designing plasmid genomes optimized for model and nonmodel organisms remains challenging. Here, we summarize current knowledge of the plasmid genome organization and the factors that can affect plasmid persistence, with the aim of constructing synthetic plasmids for use in gram-negative bacteria. Then, we introduce publicly available resources, plasmid data, and bioinformatics tools that are useful for computational plasmid design.

Euclidean Distance Analysis Enables Nucleotide Skew Analysis in Viral Genomes

Nucleotide skew analysis is a versatile method to study the nucleotide composition of RNA/DNA molecules, in particular to reveal characteristic sequence signatures. For instance, skew analysis of the nucleotide bias of several viral RNA genomes indicated that it is enriched in the unpaired, single-stranded genome regions, thus creating an even more striking virus-specific signature. The comparison of skew graphs for many virus isolates or families is difficult, time-consuming, and nonquantitative. Here, we present a procedure for a more simple identification of similarities and dissimilarities between nucleotide skew data of coronavirus, flavivirus, picornavirus, and HIV-1 RNA genomes. Window and step sizes were normalized to correct for differences in length of the viral genome. Cumulative skew data are converted into pairwise Euclidean distance matrices, which can be presented as neighbor-joining trees. We present skew value trees for the four virus families and show that closely related viruses are placed in small clusters. Importantly, the skew value trees are similar to the trees constructed by a “classical” model of evolutionary nucleotide substitution. Thus, we conclude that the simple calculation of Euclidean distances between nucleotide skew data allows an easy and quantitative comparison of characteristic sequence signatures of virus genomes. These results indicate that the Euclidean distance analysis of nucleotide skew data forms a nice addition to the virology toolbox.

Segregation but Not Replication of the Pseudomonas aeruginosa Chromosome Terminates at Dif

Coordination between chromosome replication and segregation is essential for equal partitioning of genetic material between daughter cells. In bacteria, this is achieved through the proximity of the origin of replication, oriC, and the chromosome partitioning site, parS We report here that in Pseudomonas aeruginosa, segregation but not replication is also controlled at the terminus region of the chromosome. Using the fluorescent repressor operator system (FROS), we investigated chromosome segregation in P. aeruginosa strain PAO1-UW, wherein the chromosome dimer resolution site, dif, is asymmetrically positioned relative to oriC In these cells, segregation proceeded sequentially along the two chromosomal arms and terminated at dif In contrast, chromosome replication terminated elsewhere, opposite from oriC We further found two large domains on the longer arm of the chromosome, wherein DNA segregated simultaneously. Notably, GC-skew, which reflects a bias in nucleotide usage between the leading and lagging strands of the chromosome, switches polarity at the dif locus but not necessarily at the terminus of replication. These data demonstrate that termination of chromosome replication and segregation can be physically separated without adverse effects on bacterial fitness. They also reveal the critical role of the dif region in defining the global layout of the chromosome and the progression of chromosome segregation and suggest that chromosome packing adapts to its subcellular layout.IMPORTANCE Segregation of genetic information is a central event in cellular life. In bacteria, chromosome segregation occurs concurrently with replication, sequentially along the two arms from oriC to dif How the two processes are coordinated is unknown. We explored here chromosome segregation in an opportunistic human pathogen, Pseudomonas aeruginosa, using its strain with markedly unequal chromosomal arms. We found that replication and segregation diverge in this strain and terminate at very different locations, whereas the longer chromosomal arm folds into large domains to align itself with the shorter arm. The significance of this research is in establishing that segregation and replication of bacterial chromosomes are largely uncoupled from each other and that the large-scale structure of the chromosome adapts to its subcellular layout.

Domestication of Lambda Phage Genes into a Putative Third Type of Replicative Helicase Matchmaker

At the onset of the initiation of chromosome replication, bacterial replicative helicases are recruited and loaded on the DnaA-oriC nucleoprotein platform, assisted by proteins like DnaC/DnaI or DciA. Two orders of bacteria appear, however, to lack either of these factors, raising the question of the essentiality of these factors in bacteria. Through a phylogenomic approach, we identified a pair of genes that could have substituted for dciA. The two domesticated genes are specific of the dnaC/dnaI- and dciA-lacking organisms and apparently domesticated from lambdoid phage genes. They derive from λO and λP and were renamed dopC and dopE, respectively. DopE is expected to bring the replicative helicase to the bacterial origin of replication, while DopC might assist DopE in this function. The confirmation of the implication of DopCE in the handling of the replicative helicase at the onset of replication in these organisms would generalize to all bacteria and therefore to all living organisms the need for specific factors dedicated to this function.

Recent development of Ori-Finder system and DoriC database for microbial replication origins

DNA replication begins at replication origins in all three domains of life. Identification and characterization of replication origins are important not only in providing insights into the structure and function of the replication origins but also in understanding the regulatory mechanisms of the initiation step in DNA replication. The Z-curve method has been used in the identification of replication origins in archaeal genomes successfully since 2002. Furthermore, the Web servers of Ori-Finder and Ori-Finder 2 have been developed to predict replication origins in both bacterial and archaeal genomes based on the Z-curve method, and the replication origins with manual curation have been collected into an online database, DoriC. Ori-Finder system and DoriC database are currently used in the research field of DNA replication origins in prokaryotes, including: (i) identification of oriC regions in bacterial and archaeal genomes; (ii) discovery and analysis of the conserved sequences within oriC regions; and (iii) strand-biased analysis of bacterial genomes.

Up to now, more and more predicted results by Ori-Finder system were supported by subsequent experiments, and Ori-Finder system has been used to identify the replication origins in > 100 newly sequenced prokaryotes in their genome reports. In addition, the data in DoriC database have been widely used in the large-scale analyses of replication origins and strand bias in prokaryotic genomes. Here, we review the development of Ori-Finder system and DoriC database as well as their applications. Some future directions and aspects for extending the application of Ori-Finder and DoriC are also presented.

Cultivation of Peptidiphaga gingivicola from subgingival plaque: The first representative of a novel genus of Actinomycetaceae

A novel bacterium was isolated from the subgingival plaque of a patient with periodontal disease. Bacterial strain BA112^T is a facultative Gram-positive coccus. It metabolizes alanine, arginine, glycine, histidine, leucine, proline, serine and tyrosine, but does not appear to use carbohydrates. Urease, esculin, indole, catalase and nitrate reduction tests were all negative. Major cellular fatty acids were C_18:0 , C_12:0 , C_16:0 , C_18:1 w9c and C_20:0 . The genome was sequenced and is 2.4 Mbp in length and has 64% GC content. Based on phylogenetics of the 16S rRNA sequence and concatenated alignments of 37 conserved proteins, BA112^T belongs to the family Actinomycetaceae but is located on a branch of the tree without currently named members. Based on our phenotypic and phylogenetic studies, we propose that BA112^T is the first known representative of a new genus, for which the name Peptidiphaga gingivicola gen. nov., sp. nov. is proposed. The type strain is BA112^T .

Measurement of bacterial replication rates in microbial communities

Culture-independent microbiome studies have increased our understanding of the complexity and metabolic potential of microbial communities. However, to understand the contribution of individual microbiome members to community functions, it is important to determine which bacteria are actively replicating. We developed an algorithm, iRep, that uses draft-quality genome sequences and single time-point metagenome sequencing to infer microbial population replication rates. The algorithm calculates an index of replication (iRep) based on the sequencing coverage trend that results from bi-directional genome replication from a single origin of replication. We apply this method to show that microbial replication rates increase after antibiotic administration in human infants. We also show that uncultivated, groundwater-associated, Candidate Phyla Radiation bacteria only rarely replicate quickly in subsurface communities undergoing substantial changes in geochemistry. Our method can be applied to any genome-resolved microbiome study to track organism responses to varying conditions, identify actively growing populations and measure replication rates for use in modeling studies.

Chromosomal Replication Complexity: A Novel DNA Metrics and Genome Instability Factor

As the ratio of the copy number of the most replicated to the unreplicated regions in the same chromosome, the definition of chromosomal replication complexity (CRC) appears to leave little room for variation, being either two during S-phase or one otherwise. However, bacteria dividing faster than they replicate their chromosome spike CRC to four and even eight. A recent experimental inquiry about the limits of CRC in Escherichia coli revealed two major reasons to avoid elevating it further: (i) increased chromosomal fragmentation and (ii) complications with subsequent double-strand break repair. Remarkably, examples of stable elevated CRC in eukaryotic chromosomes are well known under various terms like "differential replication," "underreplication," "DNA puffs," "onion-skin replication," or "re-replication" and highlight the phenomenon of static replication fork (sRF). To accurately describe the resulting "amplification by overinitiation," I propose a new term: "replification" (subchromosomal overreplication). In both prokaryotes and eukaryotes, replification, via sRF processing, causes double-strand DNA breaks and, with their repair elevating chromosomal rearrangements, represents a novel genome instability factor. I suggest how static replication bubbles could be stabilized and speculate that some tandem duplications represent such persistent static bubbles. Moreover, I propose how static replication bubbles could be transformed into tandem duplications, double minutes, or inverted triplications. Possible experimental tests of these models are discussed.

Neocentromeres Provide Chromosome Segregation Accuracy and Centromere Clustering to Multiple Loci along a Candida albicans Chromosome

Assembly of kinetochore complexes, involving greater than one hundred proteins, is essential for chromosome segregation and genome stability. Neocentromeres, or new centromeres, occur when kinetochores assemble de novo, at DNA loci not previously associated with kinetochore proteins, and they restore chromosome segregation to chromosomes lacking a functional centromere. Neocentromeres have been observed in a number of diseases and may play an evolutionary role in adaptation or speciation. However, the consequences of neocentromere formation on chromosome missegregation rates, gene expression, and three-dimensional (3D) nuclear structure are not well understood. Here, we used Candida albicans, an organism with small, epigenetically-inherited centromeres, as a model system to study the functions of twenty different neocentromere loci along a single chromosome, chromosome 5. Comparison of neocentromere properties relative to native centromere functions revealed that all twenty neocentromeres mediated chromosome segregation, albeit to different degrees. Some neocentromeres also caused reduced levels of transcription from genes found within the neocentromere region. Furthermore, like native centromeres, neocentromeres clustered in 3D with active/functional centromeres, indicating that formation of a new centromere mediates the reorganization of 3D nuclear architecture. This demonstrates that centromere clustering depends on epigenetically defined function and not on the primary DNA sequence, and that neocentromere function is independent of its distance from the native centromere position. Together, the results show that a neocentromere can form at many loci along a chromosome and can support the assembly of a functional kinetochore that exhibits native centromere functions including chromosome segregation accuracy and centromere clustering within the nucleus.

The accurate segregation of chromosomes during cell division is essential for maintaining genome integrity. The centromere is the DNA region on each chromosome where assembly of a large protein complex, the kinetochore, is required to maintain proper chromosome segregation. In addition, active centromeres exhibit a specific three-dimensional organization within the nucleus: the centromeres associate with one another in a clustered manner. Neocentromeres, or new centromeres, appear at new places along the chromosome when a native centromere becomes non-functional. We used a yeast model, Candida albicans, and isolated twenty instances in which neocentromeres had formed at different positions. All of these neocentromeres were able to direct chromosome segregation, but some had increased error rates. Like native centromeres, these neocentromeres cluster in the nucleus with the other active centromeres. This implies that formation of a neocentromere leads to reorganization of the three-dimensional structure of the nucleus so that different regions of the chromosome are in closer contact to regions of other chromosomes. Recent work suggests that approximately 3% of cancers may contain chromosomes with neocentromeres. Our observations that many neocentromeres have increased error rates provides insight into genome instability in cancer cells. Changes in chromosome copy number may benefit the cancer cells by increasing numbers of oncogenes and/or drug resistance genes, but may also sensitize the cells to chemotherapy approaches that target chromosome segregation mechanisms.

CRISPR-mediated control of the bacterial initiation of replication

Programmable control of the cell cycle has been shown to be a powerful tool in cell-biology studies. Here, we develop a novel system for controlling the bacterial cell cycle, based on binding of CRISPR/dCas9 to the origin-of-replication locus. Initiation of replication of bacterial chromosomes is accurately regulated by the DnaA protein, which promotes the unwinding of DNA at oriC We demonstrate that the binding of CRISPR/dCas9 to any position within origin or replication blocks the initiation of replication. Serial-dilution plating, single-cell fluorescence microscopy, and flow-cytometry experiments show that ongoing rounds of chromosome replication are finished upon CRISPR/dCas9 binding, but no new rounds are initiated. Upon arrest, cells stay metabolically active and accumulate cell mass. We find that elevating the temperature from 37 to 42°C releases the CRISR/dCas9 replication inhibition, and we use this feature to recover cells from the arrest. Our simple and robust method of controlling the bacterial cell cycle is a useful asset for synthetic biology and DNA-replication studies in particular. The inactivation of CRISPR/dCas9 binding at elevated temperatures may furthermore be of wide interest for CRISPR/Cas9 applications in genomic engineering.

Static and Dynamic Factors Limit Chromosomal Replication Complexity in Escherichia coli, Avoiding Dangers of Runaway Overreplication

We define chromosomal replication complexity (CRC) as the ratio of the copy number of the most replicated regions to that of unreplicated regions on the same chromosome. Although a typical CRC of eukaryotic or bacterial chromosomes is 2, rapidly growing Escherichia coli cells induce an extra round of replication in their chromosomes (CRC = 4). There are also E. coli mutants with stable CRC∼6. We have investigated the limits and consequences of elevated CRC in E. coli and found three limits: the "natural" CRC limit of ∼8 (cells divide more slowly); the "functional" CRC limit of ∼22 (cells divide extremely slowly); and the "tolerance" CRC limit of ∼64 (cells stop dividing). While the natural limit is likely maintained by the eclipse system spacing replication initiations, the functional limit might reflect the capacity of the chromosome segregation system, rather than dedicated mechanisms, and the tolerance limit may result from titration of limiting replication factors. Whereas recombinational repair is beneficial for cells at the natural and functional CRC limits, we show that it becomes detrimental at the tolerance CRC limit, suggesting recombinational misrepair during the runaway overreplication and giving a rationale for avoidance of the latter.

Regional centromeres in the yeast Candida lusitaniae lack pericentromeric heterochromatin

Point centromeres are specified by a short consensus sequence that seeds kinetochore formation, whereas regional centromeres lack a conserved sequence and instead are epigenetically inherited. Regional centromeres are generally flanked by heterochromatin that ensures high levels of cohesin and promotes faithful chromosome segregation. However, it is not known whether regional centromeres require pericentromeric heterochromatin. In the yeast Candida lusitaniae, we identified a distinct type of regional centromere that lacks pericentromeric heterochromatin. Centromere locations were determined by ChIP-sequencing of two key centromere proteins, Cse4 and Mif2, and are consistent with bioinformatic predictions. The centromeric DNA sequence was unique for each chromosome and spanned 4-4.5 kbp, consistent with regional epigenetically inherited centromeres. However, unlike other regional centromeres, there was no evidence of pericentromeric heterochromatin in C. lusitaniae. In particular, flanking genes were expressed at a similar level to the rest of the genome, and a URA3 reporter inserted adjacent to a centromere was not repressed. In addition, regions flanking the centromeric core were not associated with hypoacetylated histones or a sirtuin deacetylase that generates heterochromatin in other yeast. Interestingly, the centromeric chromatin had a distinct pattern of histone modifications, being enriched for methylated H3K79 and H3R2 but lacking methylation of H3K4, which is found at other regional centromeres. Thus, not all regional centromeres require flanking heterochromatin.

Genome Sequence of Mycoplasma meleagridis Type Strain 17529

Mycoplasma meleagridis is a prominent turkey bacterial pathogen associated with airsacculitis and reproductive disorders. Notwithstanding the economic losses caused by M. meleagridis, its genome has still not been sequenced. For a better understanding of its genetic background and pathogenicity mechanisms, we sequenced the genome of M. meleagridis type strain ATCC 25294.

Strand-biased gene distribution, purine assymetry and environmental factors influence protein evolution in Bacillus

A strong purine asymmetry, along with strand-biased gene distribution and the presence of PolC, prevails in Bacillus and some other members of Firmicutes, Fusobacteria and Tenericutes. The analysis of protein features in 21 Bacillus species of diverse metabolic, virulence and ecological traits revealed that purine asymmetry in conjunction with lineage/niche specific constraints significantly influences protein evolution in Bacillus. All Bacillus species, except for Se-respiring Bacillus selenitireducens, display distinct strand-specific biases in amino acid usage, which may affect the isoelectric point or surface charge distribution of proteins with prevalence of acidic and basic residues in the leading and lagging strand proteins, respectively.

Association of purine asymmetry, strand-biased gene distribution and PolC within Firmicutes and beyond: a new appraisal

Background

The Firmicutes often possess three conspicuous genome features: marked Purine Asymmetry (PAS) across two strands of replication, Strand-biased Gene Distribution (SGD) and presence of two isoforms of DNA polymerase III alpha subunit, PolC and DnaE. Despite considerable research efforts, it is not clear whether the co-existence of PAS, PolC and/or SGD is an essential and exclusive characteristic of the Firmicutes. The nature of correlations, if any, between these three features within and beyond the lineages of Firmicutes has also remained elusive. The present study has been designed to address these issues.

Results

A large-scale analysis of diverse bacterial genomes indicates that PAS, PolC and SGD are neither essential nor exclusive features of the Firmicutes. PolC prevails in four bacterial phyla: Firmicutes, Fusobacteria, Tenericutes and Thermotogae, while PAS occurs only in subsets of Firmicutes, Fusobacteria and Tenericutes. There are five major compositional trends in Firmicutes: (I) an explicit PAS or G + A-dominance along the entire leading strand (II) only G-dominance in the leading strand, (III) alternate stretches of purine-rich and pyrimidine-rich sequences, (IV) G + T dominance along the leading strand, and (V) no identifiable patterns in base usage. Presence of strong SGD has been observed not only in genomes having PAS, but also in genomes with G-dominance along their leading strands – an observation that defies the notion of co-occurrence of PAS and SGD in Firmicutes. The PolC-containing non-Firmicutes organisms often have alternate stretches of R-dominant and Y-dominant sequences along their genomes and most of them show relatively weak, but significant SGD. Firmicutes having G + A-dominance or G-dominance along LeS usually show distinct base usage patterns in three codon sites of genes. Probable molecular mechanisms that might have incurred such usage patterns have been proposed.

Conclusion

Co-occurrence of PAS, strong SGD and PolC should not be regarded as a genome signature of the Firmicutes. Presence of PAS in a species may warrant PolC and strong SGD, but PolC and/or SGD not necessarily implies PAS.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-430) contains supplementary material, which is available to authorized users.

The precarious prokaryotic chromosome

Evolutionary selection for optimal genome preservation, replication, and expression should yield similar chromosome organizations in any type of cells. And yet, the chromosome organization is surprisingly different between eukaryotes and prokaryotes. The nuclear versus cytoplasmic accommodation of genetic material accounts for the distinct eukaryotic and prokaryotic modes of genome evolution, but it falls short of explaining the differences in the chromosome organization. I propose that the two distinct ways to organize chromosomes are driven by the differences between the global-consecutive chromosome cycle of eukaryotes and the local-concurrent chromosome cycle of prokaryotes. Specifically, progressive chromosome segregation in prokaryotes demands a single duplicon per chromosome, while other "precarious" features of the prokaryotic chromosomes can be viewed as compensations for this severe restriction.

Accelerated growth in the absence of DNA replication origins

DNA replication initiates at defined sites called origins, which serve as binding sites for initiator proteins that recruit the replicative machinery. Origins differ in number and structure across the three domains of life and their properties determine the dynamics of chromosome replication. Bacteria and some archaea replicate from single origins, whereas most archaea and all eukaryotes replicate using multiple origins. Initiation mechanisms that rely on homologous recombination operate in some viruses. Here we show that such mechanisms also operate in archaea. We use deep sequencing to study replication in Haloferax volcanii and identify four chromosomal origins of differing activity. Deletion of individual origins results in perturbed replication dynamics and reduced growth. However, a strain lacking all origins has no apparent defects and grows significantly faster than wild type. Origin-less cells initiate replication at dispersed sites rather than at discrete origins and have an absolute requirement for the recombinase RadA, unlike strains lacking individual origins. Our results demonstrate that homologous recombination alone can efficiently initiate the replication of an entire cellular genome. This raises the question of what purpose replication origins serve and why they have evolved.

Genome analysis coupled with physiological studies reveals a diverse nitrogen metabolism in Methylocystis sp. strain SC2

BACKGROUND

Methylocystis sp. strain SC2 can adapt to a wide range of methane concentrations. This is due to the presence of two isozymes of particulate methane monooxygenase exhibiting different methane oxidation kinetics. To gain insight into the underlying genetic information, its genome was sequenced and found to comprise a 3.77 Mb chromosome and two large plasmids.

PRINCIPAL FINDINGS

We report important features of the strain SC2 genome. Its sequence is compared with those of seven other methanotroph genomes, comprising members of the Alphaproteobacteria, Gammaproteobacteria, and Verrucomicrobia. While the pan-genome of all eight methanotroph genomes totals 19,358 CDS, only 154 CDS are shared. The number of core genes increased with phylogenetic relatedness: 328 CDS for proteobacterial methanotrophs and 1,853 CDS for the three alphaproteobacterial Methylocystaceae members, Methylocystis sp. strain SC2 and strain Rockwell, and Methylosinus trichosporium OB3b. The comparative study was coupled with physiological experiments to verify that strain SC2 has diverse nitrogen metabolism capabilities. In correspondence to a full complement of 34 genes involved in N2 fixation, strain SC2 was found to grow with atmospheric N2 as the sole nitrogen source, preferably at low oxygen concentrations. Denitrification-mediated accumulation of 0.7 nmol (30)N2/hr/mg dry weight of cells under anoxic conditions was detected by tracer analysis. N2 production is related to the activities of plasmid-borne nitric oxide and nitrous oxide reductases.

CONCLUSIONS/PERSPECTIVES

Presence of a complete denitrification pathway in strain SC2, including the plasmid-encoded nosRZDFYX operon, is unique among known methanotrophs. However, the exact ecophysiological role of this pathway still needs to be elucidated. Detoxification of toxic nitrogen compounds and energy conservation under oxygen-limiting conditions are among the possible roles. Relevant features that may stimulate further research are, for example, absence of CRISPR/Cas systems in strain SC2, high number of iron acquisition systems in strain OB3b, and large number of transposases in strain Rockwell.

The chromosome cycle of prokaryotes

In both eukaryotes and prokaryotes, chromosomal DNA undergoes replication, condensation-decondensation and segregation, sequentially, in some fixed order. Other conditions, like sister-chromatid cohesion (SCC), may span several chromosomal events. One set of these chromosomal transactions within a single cell cycle constitutes the 'chromosome cycle'. For many years it was generally assumed that the prokaryotic chromosome cycle follows major phases of the eukaryotic one: -replication-condensation-segregation-(cell division)-decondensation-, with SCC of unspecified length. Eventually it became evident that, in contrast to the strictly consecutive chromosome cycle of eukaryotes, all stages of the prokaryotic chromosome cycle run concurrently. Thus, prokaryotes practice 'progressive' chromosome segregation separated from replication by a brief SCC, and all three transactions move along the chromosome at the same fast rate. In other words, in addition to replication forks, there are 'segregation forks' in prokaryotic chromosomes. Moreover, the bulk of prokaryotic DNA outside the replication-segregation transition stays compacted. I consider possible origins of this concurrent replication-segregation and outline the 'nucleoid administration' system that organizes the dynamic part of the prokaryotic chromosome cycle.

Neocentromeres and epigenetically inherited features of centromeres

Neocentromeres are ectopic sites where new functional kinetochores assemble and permit chromosome segregation. Neocentromeres usually form following genomic alterations that remove or disrupt centromere function. The ability to form neocentromeres is conserved in eukaryotes ranging from fungi to mammals. Neocentromeres that rescue chromosome fragments in cells with gross chromosomal rearrangements are found in several types of human cancers, and in patients with developmental disabilities. In this review, we discuss the importance of neocentromeres to human health and evaluate recently developed model systems to study neocentromere formation, maintenance, and function in chromosome segregation. Additionally, studies of neocentromeres provide insight into native centromeres; analysis of neocentromeres found in human clinical samples and induced in model organisms distinguishes features of centromeres that are dependent on centromere DNA from features that are epigenetically inherited together with the formation of a functional kinetochore.

DNA replication and strand asymmetry in prokaryotic and mitochondrial genomes

Different patterns of strand asymmetry have been documented in a variety of prokaryotic genomes as well as mitochondrial genomes. Because different replication mechanisms often lead to different patterns of strand asymmetry, much can be learned of replication mechanisms by examining strand asymmetry. Here I summarize the diverse patterns of strand asymmetry among different taxonomic groups to suggest that (1) the single-origin replication may not be universal among bacterial species as the endosymbionts Wigglesworthia glossinidia, Wolbachia species, cyanobacterium Synechocystis 6803 and Mycoplasma pulmonis genomes all exhibit strand asymmetry patterns consistent with the multiple origins of replication, (2) different replication origins in some archaeal genomes leave quite different patterns of strand asymmetry, suggesting that different replication origins in the same genome may be differentially used, (3) mitochondrial genomes from representative vertebrate species share one strand asymmetry pattern consistent with the strand-displacement replication documented in mammalian mtDNA, suggesting that the mtDNA replication mechanism in mammals may be shared among all vertebrate species, and (4) mitochondrial genomes from primitive forms of metazoans such as the sponge and hydra (representing Porifera and Cnidaria, respectively), as well as those from plants, have strand asymmetry patterns similar to single-origin or multi-origin replications observed in prokaryotes and are drastically different from mitochondrial genomes from other metazoans. This may explain why sponge and hydra mitochondrial genomes, as well as plant mitochondrial genomes, evolves much slower than those from other metazoans.

Helicobacter pylori oriC--the first bipartite origin of chromosome replication in Gram-negative bacteria

Binding of the DnaA protein to oriC leads to DNA melting within the DNA unwinding element (DUE) and initiates replication of the bacterial chromosome. Helicobacter pylori oriC was previously identified as a region localized upstream of dnaA and containing a cluster of DnaA boxes bound by DnaA protein with a high affinity. However, no unwinding within the oriC sequence has been detected. Comprehensive in silico analysis presented in this work allowed us to identify an additional region (oriC2), separated from the original one (oriC1) by the dnaA gene. DnaA specifically binds both regions, but DnaA-dependent DNA unwinding occurs only within oriC2. Surprisingly, oriC2 is bound exclusively as supercoiled DNA, which directly shows the importance of the DNA topology in DnaA-oriC interactions, similarly as previously presented only for initiator-origin interactions in Archaea and some Eukaryota. We conclude that H. pylori oriC exhibits bipartite structure, being the first such origin discovered in a Gram-negative bacterium. The H. pylori mode of initiator-oriC interactions, with the loop formation between the subcomplexes of the discontinuous origin, resembles those discovered in Bacillus subtilis chromosome and in many plasmids, which might suggest a similar way of controlling initiation of replication.

Flexibility of centromere and kinetochore structures

Centromeres, and the kinetochores that assemble on them, are essential for accurate chromosome segregation. Diverse centromere organization patterns and kinetochore structures have evolved in eukaryotes ranging from yeast to humans. In addition, centromere DNA and kinetochore position can vary even within individual cells. This flexibility is manifested in several ways: centromere DNA sequences evolve rapidly, kinetochore positions shift in response to altered chromosome structure, and kinetochore complex numbers change in response to fluctuations in kinetochore protein levels. Despite their differences, all of these diverse structures promote efficient chromosome segregation. This robustness is inherent to chromosome segregation mechanisms and balances genome stability with adaptability. In this review, we explore the mechanisms and consequences of centromere and kinetochore flexibility as well as the benefits and limitations of different experimental model systems for their study.

CAGO: a software tool for dynamic visual comparison and correlation measurement of genome organization

CAGO (Comparative Analysis of Genome Organization) is developed to address two critical shortcomings of conventional genome atlas plotters: lack of dynamic exploratory functions and absence of signal analysis for genomic properties. With dynamic exploratory functions, users can directly manipulate chromosome tracks of a genome atlas and intuitively identify distinct genomic signals by visual comparison. Signal analysis of genomic properties can further detect inconspicuous patterns from noisy genomic properties and calculate correlations between genomic properties across various genomes. To implement dynamic exploratory functions, CAGO presents each genome atlas in Scalable Vector Graphics (SVG) format and allows users to interact with it using a SVG viewer through JavaScript. Signal analysis functions are implemented using R statistical software and a discrete wavelet transformation package waveslim. CAGO is not only a plotter for generating complex genome atlases, but also a platform for exploring genome atlases with dynamic exploratory functions for visual comparison and with signal analysis for comparing genomic properties across multiple organisms. The web-based application of CAGO, its source code, user guides, video demos, and live examples are publicly available and can be accessed at http://cbs.ym.edu.tw/cago.

Complete genome and proteome of Acholeplasma laidlawii

We present the complete genome sequence and proteogenomic map for Acholeplasma laidlawii PG-8A (class Mollicutes, order Acholeplasmatales, family Acholeplasmataceae). The genome of A. laidlawii is represented by a single 1,496,992-bp circular chromosome with an average G+C content of 31 mol%. This is the longest genome among the Mollicutes with a known nucleotide sequence. It contains genes of polymerase type I, SOS response, and signal transduction systems, as well as RNA regulatory elements, riboswitches, and T boxes. This demonstrates a significant capability for the regulation of gene expression and mutagenic response to stress. Acholeplasma laidlawii and phytoplasmas are the only Mollicutes known to use the universal genetic code, in which UGA is a stop codon. Within the Mollicutes group, only the sterol-nonrequiring Acholeplasma has the capacity to synthesize saturated fatty acids de novo. Proteomic data were used in the primary annotation of the genome, validating expression of many predicted proteins. We also detected posttranslational modifications of A. laidlawii proteins: phosphorylation and acylation. Seventy-four candidate phosphorylated proteins were found: 16 candidates are proteins unique to A. laidlawii, and 11 of them are surface-anchored or integral membrane proteins, which implies the presence of active signaling pathways. Among 20 acylated proteins, 14 contained palmitic chains, and six contained stearic chains. No residue of linoleic or oleic acid was observed. Acylated proteins were components of mainly sugar and inorganic ion transport systems and were surface-anchored proteins with unknown functions.

Mapping of origin of replication in Themococcales

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

Epigenetically-inherited centromere and neocentromere DNA replicates earliest in S-phase

Eukaryotic centromeres are maintained at specific chromosomal sites over many generations. In the budding yeast Saccharomyces cerevisiae, centromeres are genetic elements defined by a DNA sequence that is both necessary and sufficient for function; whereas, in most other eukaryotes, centromeres are maintained by poorly characterized epigenetic mechanisms in which DNA has a less definitive role. Here we use the pathogenic yeast Candida albicans as a model organism to study the DNA replication properties of centromeric DNA. By determining the genome-wide replication timing program of the C. albicans genome, we discovered that each centromere is associated with a replication origin that is the first to fire on its respective chromosome. Importantly, epigenetic formation of new ectopic centromeres (neocentromeres) was accompanied by shifts in replication timing, such that a neocentromere became the first to replicate and became associated with origin recognition complex (ORC) components. Furthermore, changing the level of the centromere-specific histone H3 isoform led to a concomitant change in levels of ORC association with centromere regions, further supporting the idea that centromere proteins determine origin activity. Finally, analysis of centromere-associated DNA revealed a replication-dependent sequence pattern characteristic of constitutively active replication origins. This strand-biased pattern is conserved, together with centromere position, among related strains and species, in a manner independent of primary DNA sequence. Thus, inheritance of centromere position is correlated with a constitutively active origin of replication that fires at a distinct early time. We suggest a model in which the distinct timing of DNA replication serves as an epigenetic mechanism for the inheritance of centromere position.

Centromeres form at the same chromosomal position from generation to generation, yet in most species this inheritance occurs in a DNA sequence–independent manner that is not well understood. Here, we determine the timing of DNA replication across the genome of the human fungal pathogen Candida albicans and find that centromeric DNA is the first locus to replicate on each chromosome. Furthemore, this unique replication timing may be important for centromere inheritance, based on several observations. First, DNA sequence patterns at centromeres indicate that, despite high levels of primary sequence divergence, the region has served as a replication origin for millions of years; second, formation of a neocentromere (a new centromere formed at an ectopic locus following deletion of the native centromere DNA) results in the establishment of a new, early-firing origin of replication; and third, a centromere-specific protein, Cse4p, recruits origin replication complex proteins in a concentration-dependent manner. Thus, centromere position is inherited by an epigenetic mechanism that appears to be defined by a distinctively early firing DNA replication origin.

Nitrosopumilus maritimus genome reveals unique mechanisms for nitrification and autotrophy in globally distributed marine crenarchaea

Ammonia-oxidizing archaea are ubiquitous in marine and terrestrial environments and now thought to be significant contributors to carbon and nitrogen cycling. The isolation of Candidatus "Nitrosopumilus maritimus" strain SCM1 provided the opportunity for linking its chemolithotrophic physiology with a genomic inventory of the globally distributed archaea. Here we report the 1,645,259-bp closed genome of strain SCM1, revealing highly copper-dependent systems for ammonia oxidation and electron transport that are distinctly different from known ammonia-oxidizing bacteria. Consistent with in situ isotopic studies of marine archaea, the genome sequence indicates N. maritimus grows autotrophically using a variant of the 3-hydroxypropionate/4-hydroxybutryrate pathway for carbon assimilation, while maintaining limited capacity for assimilation of organic carbon. This unique instance of archaeal biosynthesis of the osmoprotectant ectoine and an unprecedented enrichment of multicopper oxidases, thioredoxin-like proteins, and transcriptional regulators points to an organism responsive to environmental cues and adapted to handling reactive copper and nitrogen species that likely derive from its distinctive biochemistry. The conservation of N. maritimus gene content and organization within marine metagenomes indicates that the unique physiology of these specialized oligophiles may play a significant role in the biogeochemical cycles of carbon and nitrogen.

Crystal structure of ArgP from Mycobacterium tuberculosis confirms two distinct conformations of full-length LysR transcriptional regulators and reveals its function in DNA binding and transcriptional regulation

Mycobacterium tuberculosis presents a challenging medical problem partly due to its persistent nonreplicative state. The inhibitor of chromosomal replication (iciA) protein encoded by M. tuberculosis has been suggested to inhibit chromosome replication initiation in vitro. However, iciA has also been identified as arginine permease (ArgP), a regulatory transcription factor for arginine outward transport. In order to understand the function of ArgP, we have determined its crystal structure by X-ray crystallography to a resolution of 2.7 A. ArgP is a member of the LysR-type transcriptional regulators (LTTRs) and forms a homodimer with each subunit containing two domains: a DNA binding domain (DBD) and a regulatory domain (RD). Two conformationally distinct subunits were identified: closed subunit and open subunit. This phenomenon was first observed in LTTR CbnR, but not in LTTR CrgA, and might be common in LTTRs. We identified two forms of dimers: DBD-type dimers and RD-type dimers. The former is confirmed in solution, and the latter is considered to form oligomers during function. We provide the first structural insights into the interaction of the extreme C-terminal residues with the DBD, which is confirmed by mutagenesis and analytical ultracentrifugation to be important for stability of the functional dimer. The structure serves as a model to suggest how three critical aspects, namely, DNA binding, homo-oligomerization, and interaction with RNAP, are mediated during regulation processing. A model is proposed for the LysR family of dimeric regulators.

A study in entire chromosomes of violations of the intra-strand parity of complementary nucleotides (Chargaff's second parity rule)

Chargaff's rule of intra-strand parity (ISP) between complementary mono/oligonucleotides in chromosomes is well established in the scientific literature. Although a large numbers of papers have been published citing works and discussions on ISP in the genomic era, scientists are yet to find all the factors responsible for such a universal phenomenon in the chromosomes. In the present work, we have tried to address the issue from a new perspective, which is a parallel feature to ISP. The compositional abundance values of mono/oligonucleotides were determined in all non-overlapping sub-chromosomal regions of specific size. Also the frequency distributions of the mono/oligonucleotides among the regions were compared using the Kolmogorov–Smirnov test. Interestingly, the frequency distributions between the complementary mono/oligonucleotides revealed statistical similarity, which we named as intra-strand frequency distribution parity (ISFDP). ISFDP was observed as a general feature in chromosomes of bacteria, archaea and eukaryotes. Violation of ISFDP was also observed in several chromosomes. Chromosomes of different strains belonging a species in bacteria/archaea (Haemophilus influenza, Xylella fastidiosa etc.) and chromosomes of a eukaryote are found to be different among each other with respect to ISFDP violation. ISFDP correlates weakly with ISP in chromosomes suggesting that the latter one is not entirely responsible for the former. Asymmetry of replication topography and composition of forward-encoded sequences between the strands in chromosomes are found to be insufficient to explain the ISFDP feature in all chromosomes. This suggests that multiple factors in chromosomes are responsible for establishing ISFDP.