Comprehensive Analysis of Replication Origins in Saccharomyces cerevisiae Genomes

Identification of 1600 replication origins in S. cerevisiae

There are approximately 500 known origins of replication in the yeast genome, and the process by which DNA replication initiates at these locations is well understood. In particular, these sites are made competent to initiate replication by loading of the Mcm replicative helicase prior to the start of S phase; thus, 'a site that binds Mcm in G1' might be considered to provide an operational definition of a replication origin. By fusing a subunit of Mcm to micrococcal nuclease, we previously showed that known origins are typically bound by a single Mcm double hexamer, loaded adjacent to the ARS consensus sequence (ACS). Here, we extend this analysis from known origins to the entire genome, identifying candidate Mcm binding sites whose signal intensity varies over at least three orders of magnitude. Published data quantifying single-stranded DNA (ssDNA) during S phase revealed replication initiation among the most abundant 1600 of these sites, with replication activity decreasing with Mcm abundance and disappearing at the limit of detection of ssDNA. Three other hallmarks of replication origins were apparent among the most abundant 5500 sites. Specifically, these sites: (1) appeared in intergenic nucleosome-free regions flanked on one or both sides by well-positioned nucleosomes; (2) were flanked by ACSs; and (3) exhibited a pattern of GC skew characteristic of replication initiation. We conclude that, if sites at which Mcm double hexamers are loaded can function as replication origins, then DNA replication origins are at least threefold more abundant than previously assumed, and we suggest that replication may occasionally initiate in essentially every intergenic region. These results shed light on recent reports that as many as 15% of replication events initiate outside of known origins, and this broader distribution of replication origins suggest that S phase in yeast may be less distinct from that in humans than widely assumed.

Identification of 1600 replication origins in S. cerevisiae

There are approximately 500 known origins of replication in the yeast genome, and the process by which DNA replication initiates at these locations is well understood. In particular, these sites are made competent to initiate replication by loading of the Mcm replicative helicase prior to the start of S phase; thus, "a site to which MCM is bound in G1" might be considered to provide an operational definition of a replication origin. By fusing a subunit of Mcm to micrococcal nuclease, a technique referred to as "Chromatin Endogenous Cleavage", we previously showed that known origins are typically bound by a single Mcm double hexamer, loaded adjacent to the ARS consensus sequence (ACS). Here we extend this analysis from known origins to the entire genome, identifying candidate Mcm binding sites whose signal intensity varies over at least 3 orders of magnitude. Published data quantifying the production of ssDNA during S phase showed clear evidence of replication initiation among the most abundant 1600 of these sites, with replication activity decreasing in concert with Mcm abundance and disappearing at the limit of detection of ssDNA. Three other hallmarks of replication origins were apparent among the most abundant 5,500 sites. Specifically, these sites (1) appeared in intergenic nucleosome-free regions that were flanked on one or both sides by well-positioned nucleosomes; (2) were flanked by ACSs; and (3) exhibited a pattern of GC skew characteristic of replication initiation. Furthermore, the high resolution of this technique allowed us to demonstrate a strong bias for detecting Mcm double-hexamers downstream rather than upstream of the ACS, which is consistent with the directionality of Mcm loading by Orc that has been observed in vitro. We conclude that, if sites at which Mcm double-hexamers are loaded can function as replication origins, then DNA replication origins are at least 3-fold more abundant than previously assumed, and we suggest that replication may occasionally initiate in essentially every intergenic region. These results shed light on recent reports that as many as 15% of replication events initiate outside of known origins, and this broader distribution of replication origins suggest that S phase in yeast may be less distinct from that in humans than is widely assumed.

DAPI staining and DNA content estimation of nuclei in uncultivable microbial eukaryotes (Arcellinida and Ciliates)

Though representing a major component of eukaryotic biodiversity, many microbial eukaryotes remain poorly studied, including the focus of the present work, testate amoebae of the order Arcellinida (Amoebozoa) and non-model lineages of ciliates (Alveolata). In particular, knowledge of genome structures and changes in genome content over the often-complex life cycles of these lineages remains enigmatic. However, the limited available knowledge suggests that microbial eukaryotes have the potential to challenge our textbook views on eukaryotic genomes and genome evolution. In this study, we developed protocols for DAPI (4',6-diamidino-2-phenylindole) staining of Arcellinida nuclei and adapted protocols for ciliates. In addition, image analysis software was used to estimate the DNA content in the nuclei of Arcellinida and ciliates, and the measurements of target organisms were compared to those  of well-known model organisms.The results demonstrate that the methods we have developed for nuclear staining in these lineages are effective and can be applied to other microbial eukaryotic groups by adjusting certain stages in the protocols.

OriC-ENS: A sequence-based ensemble classifier for predicting origin of replication in S. cerevisiae

DNA Replication plays the most crucial part in biological inheritance, ensuring an even flow of genetic information from parent to offspring. The beginning site of DNA Replication which is called the Origin of Replication (ORI), plays a significant role in understanding the molecular mechanisms and genomic analysis of DNA. Hence, it is paramount to accurately identify the origin of replication to gain a more accurate understanding of the biochemical and genomic properties of DNA. In this paper, We have proposed a new approach named OriC-ENS that uses sequence-based feature extraction techniques, K-mer, K-gapped Mono-Di, and Di Mono, and an ensemble classification technique that uses majority voting for the identification of Origin of Replication. We have used three SVM classifiers, one for the K-mer features and two more for K-Gapped Mono-Di and K-Gapped Di-mono features. Finally, we used majority voting to combine the prediction by each predictor. Experimental results on the S. Cerevisiae dataset have shown that our method achieves an accuracy of 91.62 % which outperforms other state-of-the-art methods by a significant margin. We have also tested our method using other evaluation metrics such as Matthews Correlation Coefficient (MCC), Area Under Curve(AUC), Sensitivity, and Specificity, where it has achieved a score of 0.83, 0.98, 0.90, and 0.92 respectively. We have further evaluated our model on an independent test set collected from OriDB, consisting of the sequences of Schizosaccharomyces pombe where we have seen that our model can predict the origin of replication efficiently and with great precision. We have made our python-based source code available at https://github.com/MehediAzim/OriC-ENS.

Refactoring the Conjugation Machinery of Promiscuous Plasmid RP4 into a Device for Conversion of Gram-Negative Isolates to Hfr Strains

Chromosomal exchange and subsequent recombination of the cognate DNA between bacteria was one of the most useful genetic tools (e.g., Hfr strains) for genetic analyses of E. coli before the genomic era. In this paper, yeast assembly has been used to recruit the conjugation machinery of environmentally promiscuous RP4 plasmid into a minimized, synthetic construct that enables transfer of chromosomal segments between donor/recipient strains of P. putida KT2440 and potentially many other Gram-negative bacteria. The synthetic device features [i] a R6K suicidal plasmid backbone, [ii] a mini-Tn5 transposon vector, and [iii] the minimal set of genes necessary for active conjugation (RP4 Tra1 and Tra2 clusters) loaded as cargo in the mini-Tn5 mobile element. Upon insertion of the transposon in different genomic locations, the ability of P. putida-TRANS (transference of RP4-activated nucleotide segments) donor strains to mobilize genomic stretches of DNA into neighboring bacteria was tested. To this end, a P. putida double mutant ΔpyrF (uracil auxotroph) Δedd (unable to grow on glucose) was used as recipient in mating experiments, and the restoration of the pyrF+/edd+ phenotypes allowed for estimation of chromosomal transfer efficiency. Cells with the inserted transposon behaved in a manner similar to Hfr-like strains and were able to transfer up to 23% of their genome at frequencies close to 10-6 exconjugants per recipient cell. The hereby described TRANS device not only expands the molecular toolbox for P. putida, but it also enables a suite of genomic manipulations which were thus far only possible with domesticated laboratory strains and species.

Identification and Characterization of the Mitochondrial Replication Origin for Stable and Episomal Expression in Yarrowia lipolytica

Episomal plasmids are crucial expression tools for recombinant protein production and genome editing. In Saccharomyces cerevisiae, 2-μm artificial plasmids with a high copy number have been developed and used in metabolic engineering and synthetic biology. However, in unconventional yeasts such as Yarrowia lipolytica, episomal expression relies on a chromosome replication system; this system has the disadvantages of genetic instability and low copy numbers. In this study, we identified and characterized replication origins from the mitochondrial DNA (mtDNA) of Y. lipolytica. A 516-bp mtDNA sequence, mtORI, was confirmed to mediate the autonomous replication of circular plasmids with high protein expression levels and hereditary stability. However, the nonhomologous end-joining pathway could interfere with mtORI plasmid replication and engender genetic heterogeneity. In the Po 1f ΔKu70 strain, the homogeneity of the mtORI plasmid was significantly improved, and the highest copy number reached 5.0 per cell. Overall, mitochondrial-origin sequences can be used to establish highly stable and autonomously replicating plasmids, which can be a powerful supplement to the current synthetic biology tool library and promote the development of Y. lipolytica as a microbial cell factory.

A deep learning framework combined with word embedding to identify DNA replication origins

The DNA replication influences the inheritance of genetic information in the DNA life cycle. As the distribution of replication origins (ORIs) is the major determinant to precisely regulate the replication process, the correct identification of ORIs is significant in giving an insightful understanding of DNA replication mechanisms and the regulatory mechanisms of genetic expressions. For eukaryotes in particular, multiple ORIs exist in each of their gene sequences to complete the replication in a reasonable period of time. To simplify the identification process of eukaryote's ORIs, most of existing methods are developed by traditional machine learning algorithms, and target to the gene sequences with a fixed length. Consequently, the identification results are not satisfying, i.e. there is still great room for improvement. To break through the limitations in previous studies, this paper develops sequence segmentation methods, and employs the word embedding technique, 'Word2vec', to convert gene sequences into word vectors, thereby grasping the inner correlations of gene sequences with different lengths. Then, a deep learning framework to perform the ORI identification task is constructed by a convolutional neural network with an embedding layer. On the basis of the analysis of similarity reduction dimensionality diagram, Word2vec can effectively transform the inner relationship among words into numerical feature. For four species in this study, the best models are obtained with the overall accuracy of 0.975, 0.765, 0.885, 0.967, the Matthew's correlation coefficient of 0.940, 0.530, 0.771, 0.934, and the AUC of 0.975, 0.800, 0.888, 0.981, which indicate that the proposed predictor has a stable ability and provide a high confidence coefficient to classify both of ORIs and non-ORIs. Compared with state-of-the-art methods, the proposed predictor can achieve ORI identification with significant improvement. It is therefore reasonable to anticipate that the proposed method will make a useful high throughput tool for genome analysis.

Strong preference for the integration of transforming DNA via homologous recombination in Trichoderma atroviride

Trichoderma species play important roles in nature as plant growth promotors and antagonists of phytopathogenic fungi, and are used as models to study photomorphogenesis. Molecular tools have been implemented to manipulate and improve these fungi. However, instability of transformants or very low frequency of homologous recombination has been reported. Here, we report the fate of transforming DNA, demonstrating that it can follow two different fates. When a vector contains sequences also present in the Trichodermaatroviride genome, it mainly integrates by homologous recombination generating stable recombinant strains. In contrast, vectors with no sequence homology to the T. atroviride genome generate unstable transformants, losing the transforming DNA in the first generation of conidia produced without selection where, surprisingly, the vector behaves as autoreplicative. Integration by homologous recombination was demonstrated when transformants were generated with a truncated version of the blr2 gene, resulting in insertional mutants with phenotypes identical to those of knockout mutants. Our results indicate that T. atroviride is highly efficient in integrating DNA by homologous recombination and that plasmid vectors with no sequence homology to the genome are maintained for several generations in T. atroviride if kept under selective pressure even though they lacked fungal autonomous replication sequences.

Ori-Finder 3: a web server for genome-wide prediction of replication origins in Saccharomyces cerevisiae

DNA replication is a fundamental process in all organisms; this event initiates at sites termed origins of replication. The characteristics of eukaryotic replication origins are best understood in Saccharomyces cerevisiae. For this species, origin prediction algorithms or web servers have been developed based on the sequence features of autonomously replicating sequences (ARSs). However, their performances are far from satisfactory. By utilizing the Z-curve methodology, we present a novel pipeline, Ori-Finder 3, for the computational prediction of replication origins in S. cerevisiae at the genome-wide level based solely on DNA sequences. The ARS exhibiting both an AT-rich stretch and ARS consensus sequence element can be predicted at the single-nucleotide level. For the identified ARSs in the S. cerevisiae reference genome, 83 and 60% of the top 100 and top 300 predictions matched the known ARS records, respectively. Based on Ori-Finder 3, we subsequently built a database of the predicted ARSs identified in more than a hundred S. cerevisiae genomes. Consequently, we developed a user-friendly web server including the ARS prediction pipeline and the predicted ARSs database, which can be freely accessed at http://tubic.tju.edu.cn/Ori-Finder3.

Repair characteristics and time-dependent effects in response to heavy-ion beam irradiation in Saccharomyces cerevisiae: a comparison with X-ray irradiation

Heavy-ion beam (HIB) irradiation has been widely used in microbial mutation breeding. However, a global cellular response to such radiation remains mostly uncharacterised. In this study, we used transcriptomics to analyse the damage repair response in Saccharomyces cerevisiae following a semi-lethal HIB irradiation (80 Gy), which induced a significant number of DNA double-strand breaks. Our analysis of differentially expressed genes (DEGs) from 50 to 150 min post-irradiation revealed that upregulated genes were significantly enriched for gene ontology and Kyoto encyclopaedia of genes and genomes terms related to damage repair response. Based on the number of DEGs, their annotation, and their relative expression, we established that the peak of the damage repair response occurred 75 min post-irradiation. Moreover, we exploited the data from our recent study on X-ray irradiation-induced repair to compare the transcriptional patterns induced by semi-lethal HIB and X-ray irradiations. Although these two radiations have different properties, we found a significant overlap (> 50%) for the DEGs associated with five typical DNA repair pathways and, in both cases, identified homologous recombination repair (HRR) as the predominant repair pathway. Nevertheless, when we compared the relative enrichment of the five DNA repair pathways at the key time point of the repair process, we found that the relative enrichment of HRR was higher after HIB irradiation than after X-ray irradiation. Additionally, the peak stage of HRR following HIB irradiation was ahead of that following X-ray irradiation. Since mutations occur during the DNA repair process, uncovering detailed repair characteristics should further the understanding of the associated mutagenesis features.