Recent Advances in the Identification of Replication Origins Based on the Z-curve Method

Unveiling human origins of replication using deep learning: accurate prediction and comprehensive analysis

Accurate identification of replication origins (ORIs) is crucial for a comprehensive investigation into the progression of human cell growth and cancer therapy. Here, we proposed a computational approach Ori-FinderH, which can efficiently and precisely predict the human ORIs of various lengths by combining the Z-curve method with deep learning approach. Compared with existing methods, Ori-FinderH exhibits superior performance, achieving an area under the receiver operating characteristic curve (AUC) of 0.9616 for K562 cell line in 10-fold cross-validation. In addition, we also established a cross-cell-line predictive model, which yielded a further improved AUC of 0.9706. The model was subsequently employed as a fitness function to support genetic algorithm for generating artificial ORIs. Sequence analysis through iORI-Euk revealed that a vast majority of the created sequences, specifically 98% or more, incorporate at least one ORI for three cell lines (Hela, MCF7 and K562). This innovative approach could provide more efficient, accurate and comprehensive information for experimental investigation, thereby further advancing the development of this field.

Ori-Finder 2022: A Comprehensive Web Server for Prediction and Analysis of Bacterial Replication Origins

The replication of DNA is a complex biological process that is essential for life. Bacterial DNA replication is initiated at genomic loci referred to as replication origins (oriCs). Integrating the Z-curve method, DnaA box distribution, and comparative genomic analysis, we developed a web server to predict bacterial oriCs in 2008 called Ori-Finder, which contributes to clarify the characteristics of bacterial oriCs. The oriCs of hundreds of sequenced bacterial genomes have been annotated in the genome reports using Ori-Finder and the predicted results have been deposited in DoriC, a manually curated database of oriCs. This has facilitated large-scale data mining of functional elements in oriCs and strand-biased analysis. Here, we describe Ori-Finder 2022 with updated prediction framework, interactive visualization module, new analysis module, and user-friendly interface. More species-specific indicator genes and functional elements of oriCs are integrated into the updated framework, which has also been redesigned to predict oriCs in draft genomes. The interactive visualization module displays more genomic information related to oriCs and their functional elements. The analysis module includes regulatory protein annotation, repeat sequence discovery, homologous oriC search, and strand-biased analyses. The redesigned interface provides additional customization options for oriC prediction. Ori-Finder 2022 is freely available at http://tubic.tju.edu.cn/Ori-Finder/ and https://tubic.org/Ori-Finder/.

Computational characterization and analysis of molecular sequence data of Elizabethkingia meningoseptica

OBJECTIVE

Elizabethkingia meningoseptica is a multidrug resistance strain which primarily causes meningitis in neonates and immunocompromised patients. Being a nosocomial infection causing agent, less information is available in literature, specifically, about its genomic makeup and associated features. An attempt is made to study them through bioinformatics tools with respect to compositions, embedded periodicities, open reading frames, origin of replication, phylogeny, orthologous gene clusters analysis and pathways.

RESULTS

Complete DNA and protein sequence pertaining to E. meningoseptica were thoroughly analyzed as part of the study. E. meningoseptica G4076 genome showed 7593 ORFs it is GC rich. Fourier based analysis showed the presence of typical three base periodicity at the genome level. Putative origin of replication has been identified. Phylogenetically, E. meningoseptica is relatively closer to E. anophelis compared to other Elizabethkingia species. A total of 2606 COGs were shared by all five Elizabethkingia species. Out of 3391 annotated proteins, we could identify 18 unique ones involved in metabolic pathway of E. meningoseptica and this can be an initiation point for drug designing and development. Our study is novel in the aspect in characterizing and analyzing the whole genome data of E. meningoseptica.

Cyanobacterial multi-copy chromosomes and their replication

While the model bacteria Escherichia coli and Bacillus subtilis harbor single chromosomes, which is known as monoploidy, some freshwater cyanobacteria contain multiple chromosome copies per cell throughout their cell cycle, which is known as polyploidy. In the model cyanobacteria Synechococcus elongatus PCC 7942 and Synechocystis sp. PCC 6803, chromosome copy number (ploidy) is regulated in response to growth phase and environmental factors. In S. elongatus 7942, chromosome replication is asynchronous both among cells and chromosomes. Comparative analysis of S. elongatus 7942 and S. sp. 6803 revealed a variety of DNA replication mechanisms. In this review, the current knowledge of ploidy and DNA replication mechanisms in cyanobacteria is summarized together with information on the features common with plant chloroplasts. It is worth noting that the occurrence of polyploidy and its regulation are correlated with certain cyanobacterial lifestyles and are shared between some cyanobacteria and chloroplasts.

ABBREVIATIONS

NGS: next-generation sequencing; Repli-seq: replication sequencing; BrdU: 5-bromo-2'-deoxyuridine; TK: thymidine kinase; GCSI: GC skew index; PET: photosynthetic electron transport; RET: respiration electron transport; Cyt b₆f complex: cytochrome b₆f complex; PQ: plastoquinone; PC: plastocyanin.

Quantitative analysis and assessment of base composition asymmetry and gene orientation bias in bacterial genomes

Base composition asymmetry and gene orientation bias are two common genomic structures in bacterial genomes. Here, correlation coefficients between nucleotide disparities and coding sequence (CDS) skew have been calculated, which provides insights into their relationship from an individual genome perspective. Consequently, we find GC and RY disparities correlate significantly with CDS skew, since around 60% of the bacterial genomes under study have correlation coefficients > 0.9. Then, we present a model for quantitative assessment of nucleotide disparity and CDS skew in which a numerical index R² is used for evaluation. We find that skew curves with higher R² perform better on the prediction of replication origins in bacteria.

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.

Quantitative analysis of correlation between AT and GC biases among bacterial genomes

Due to different replication mechanisms between the leading and lagging strands, nucleotide composition asymmetries widely exist in bacterial genomes. A general consideration reveals that the leading strand is enriched in Guanine (G) and Thymine (T), and the lagging strand shows richness in Adenine (A) and Cytosine (C). However, some bacteria like Bacillus subtilis have been discovered composing more A than T in the leading strand. To investigate the difference, we analyze the nucleotide asymmetry from the aspect of AT and GC bias correlations. In this study, we propose a windowless method, the Z-curve Correlation Coefficient (ZCC) index, based on the Z-curve method, and analyzed more than 2000 bacterial genomes. We find that the majority of bacteria reveal negative correlations between AT and GC biases, while most genomes in Firmicutes and Tenericutes have positive ZCC indexes. The presence of PolC, purine asymmetry and stronger genes preference in the leading strand are not confined to Firmicutes, but also likely to happen in other phyla dominated by positive ZCC indexes. This method also provides a new insight into other relevant features like aerobism, and can be applied to analyze the correlation between RY (Purine and Pyrimidine) and MK (Amino and Keto) bias and so on.

Context based computational analysis and characterization of ARS consensus sequences (ACS) of Saccharomyces cerevisiae genome

Genome-wide experimental studies in Saccharomyces cerevisiae reveal that autonomous replicating sequence (ARS) requires an essential consensus sequence (ACS) for replication activity. Computational studies identified thousands of ACS like patterns in the genome. However, only a few hundreds of these sites act as replicating sites and the rest are considered as dormant or evolving sites. In a bid to understand the sequence makeup of replication sites, a content and context-based analysis was performed on a set of replicating ACS sequences that binds to origin-recognition complex (ORC) denoted as ORC-ACS and non-replicating ACS sequences (nrACS), that are not bound by ORC. In this study, DNA properties such as base composition, correlation, sequence dependent thermodynamic and DNA structural profiles, and their positions have been considered for characterizing ORC-ACS and nrACS. Analysis reveals that ORC-ACS depict marked differences in nucleotide composition and context features in its vicinity compared to nrACS. Interestingly, an A-rich motif was also discovered in ORC-ACS sequences within its nucleosome-free region. Profound changes in the conformational features, such as DNA helical twist, inclination angle and stacking energy between ORC-ACS and nrACS were observed. Distribution of ACS motifs in the non-coding segments points to the locations of ORC-ACS which are found far away from the adjacent gene start position compared to nrACS thereby enabling an accessible environment for ORC-proteins. Our attempt is novel in considering the contextual view of ACS and its flanking region along with nucleosome positioning in the S. cerevisiae genome and may be useful for any computational prediction scheme.

Bioinformatics-Based Molecular Classification of Arthrobacter Plasmids

The omnipresence of Arthrobacter species in polluted and toxic soils indicates their great potential in environmental biotechnologies, but practical applications of these bacteria are scarce mainly due to the availability of useful genetic engineering tools. Although many fully sequenced Arthrobacter genomes have been deposited in GenBank, little is known about the biology of their plasmids, especially the core functions: replication and partition. In this study the available Arthrobacter plasmid sequences were analyzed in order to identify their putative replication origin. At least the oris from the cryptic plasmids pXZ10142, pCG1, and pBL1 appear to work in this genus. Based on ParA homolog sequences, the Arthrobacter specific plasmids were classified into 4 clades. Iteron-like sequences were identified on most of the plasmids, indicating the position of the putative Arthrobacter specific oris. Although attempts were made to identify the core gene set required for plasmid replication in this genus, it was not possible. The plasmid proteomes showed a rather low similarity.

Acidianus Tailed Spindle Virus: a New Archaeal Large Tailed Spindle Virus Discovered by Culture-Independent Methods

The field of viral metagenomics has expanded our understanding of viral diversity from all three domains of life (Archaea, Bacteria, and Eukarya). Traditionally, viral metagenomic studies provide information about viral gene content but rarely provide knowledge about virion morphology and/or cellular host identity. Here we describe a new virus, Acidianus tailed spindle virus (ATSV), initially identified by bioinformatic analysis of viral metagenomic data sets from a high-temperature (80°C) acidic (pH 2) hot spring located in Yellowstone National Park, followed by more detailed characterization using only environmental samples without dependency on culturing. Characterization included the identification of the large tailed spindle virion morphology, determination of the complete 70.8-kb circular double-stranded DNA (dsDNA) viral genome content, and identification of its cellular host. Annotation of the ATSV genome revealed a potential three-domain gene product containing an N-terminal leucine-rich repeat domain, followed by a likely posttranslation regulatory region consisting of high serine and threonine content, and a C-terminal ESCRT-III domain, suggesting interplay with the host ESCRT system. The host of ATSV, which is most closely related to Acidianus hospitalis, was determined by a combination of analysis of cellular clustered regularly interspaced short palindromic repeat (CRISPR)/Cas loci and dual viral and cellular fluorescence in situ hybridization (viral FISH) analysis of environmental samples and confirmed by culture-based infection studies. This work provides an expanded pathway for the discovery, isolation, and characterization of new viruses using culture-independent approaches and provides a platform for predicting and confirming virus hosts.

IMPORTANCE

Virus discovery and characterization have been traditionally accomplished by using culture-based methods. While a valuable approach, it is limited by the availability of culturable hosts. In this research, we report a virus-centered approach to virus discovery and characterization, linking viral metagenomic sequences to a virus particle, its sequenced genome, and its host directly in environmental samples, without using culture-dependent methods. This approach provides a pathway for the discovery, isolation, and characterization of new viruses. While this study used an acidic hot spring environment to characterize a new archaeal virus, Acidianus tailed spindle virus (ATSV), the approach can be generally applied to any environment to expand knowledge of virus diversity in all three domains of life.

Identification of the Replication Origins from Cyanothece ATCC 51142 and Their Interactions with the DnaA Protein: From In Silico to In Vitro Studies

Based on the complete genome of Cyanothece ATCC 51142, the oriCs of both the circular and linear chromosomes in Cyanothece ATCC 51142 have been predicted by utilizing a web-based system Ori-Finder. Here, we provide experimental support for the results of Ori-Finder to identify the replication origins of Cyanothece ATCC 51142 and their interactions with the initiator protein, DnaA. The two replication origins are composed of three characteristically arranged DnaA boxes and an AT-rich stretch, and the oriC in the circular chromosome is followed by the dnaN gene. The dnaA gene is located downstream of the origin of the circular chromosome and it expresses a typical DnaA protein that is divided into four domains (I, II, III, IV), as with other members of the DnaA protein family. We purify DnaA (IV) and characterize the interaction of the purified protein with the replication origins, so as to offer experimental support for the prediction. The results of the electrophoretic mobility shift assay and DNase I footprint assay demonstrate that the C-terminal domain of the DnaA protein from Cyanothece ATCC 51142 specifically binds the oriCs of both the circular and linear chromosomes, and the DNase I footprint assay demonstrates that DnaA (IV) exhibits hypersensitive affinity with DnaA boxes in both oriCs.

Diversification of DnaA dependency for DNA replication in cyanobacterial evolution

Regulating DNA replication is essential for all living cells. The DNA replication initiation factor DnaA is highly conserved in prokaryotes and is required for accurate initiation of chromosomal replication at oriC. DnaA-independent free-living bacteria have not been identified. The dnaA gene is absent in plastids and some symbiotic bacteria, although it is not known when or how DnaA-independent mechanisms were acquired. Here, we show that the degree of dependency of DNA replication on DnaA varies among cyanobacterial species. Deletion of the dnaA gene in Synechococcus elongatus PCC 7942 shifted DNA replication from oriC to a different site as a result of the integration of an episomal plasmid. Moreover, viability during the stationary phase was higher in dnaA disruptants than in wild-type cells. Deletion of dnaA did not affect DNA replication or cell growth in Synechocystis sp. PCC 6803 or Anabaena sp. PCC 7120, indicating that functional dependency on DnaA was already lost in some nonsymbiotic cyanobacterial lineages during diversification. Therefore, we proposed that cyanobacteria acquired DnaA-independent replication mechanisms before symbiosis and such an ancestral cyanobacterium was the sole primary endosymbiont to form a plastid precursor.

Choosing a suitable method for the identification of replication origins in microbial genomes

As the replication of genomic DNA is arguably the most important task performed by a cell and given that it is controlled at the initiation stage, the events that occur at the replication origin play a central role in the cell cycle. Making sense of DNA replication origins is important for improving our capacity to study cellular processes and functions in the regulation of gene expression, genome integrity in much finer detail. Thus, clearly comprehending the positions and sequences of replication origins which are fundamental to chromosome organization and duplication is the first priority of all. In view of such important roles of replication origins, tremendous work has been aimed at identifying and testing the specificity of replication origins. A number of computational tools based on various skew types have been developed to predict replication origins. Using various in silico approaches such as Ori-Finder, and databases such as DoriC, researchers have predicted the locations of replication origins sites for thousands of bacterial chromosomes and archaeal genomes. Based on the predicted results, we should choose an effective method for identifying and confirming the interactions at origins of replication. Here we describe the main existing experimental methods that aimed to determine the replication origin regions and list some of the many the practical applications of these methods.

Ori-Finder 2, an integrated tool to predict replication origins in the archaeal genomes

DNA replication is one of the most basic processes in all three domains of cellular life. With the advent of the post-genomic era, the increasing number of complete archaeal genomes has created an opportunity for exploration of the molecular mechanisms for initiating cellular DNA replication by in vivo experiments as well as in silico analysis. However, the location of replication origins (oriCs) in many sequenced archaeal genomes remains unknown. We present a web-based tool Ori-Finder 2 to predict oriCs in the archaeal genomes automatically, based on the integrated method comprising the analysis of base composition asymmetry using the Z-curve method, the distribution of origin recognition boxes identified by FIMO tool, and the occurrence of genes frequently close to oriCs. The web server is also able to analyze the unannotated genome sequences by integrating with gene prediction pipelines and BLAST software for gene identification and function annotation. The result of the predicted oriCs is displayed as an HTML table, which offers an intuitive way to browse the result in graphical and tabular form. The software presented here is accurate for the genomes with single oriC, but it does not necessarily find all the origins of replication for the genomes with multiple oriCs. Ori-Finder 2 aims to become a useful platform for the identification and analysis of oriCs in the archaeal genomes, which would provide insight into the replication mechanisms in archaea. The web server is freely available at http://tubic.tju.edu.cn/Ori-Finder2/.