Capturing the dynamics of genome replication on individual ultra-long nanopore sequence reads

Quantifying DNA replication speeds in single cells by scEdU-seq

In a human cell, thousands of replication forks simultaneously coordinate duplication of the entire genome. The rate at which this process occurs might depend on the epigenetic state of the genome and vary between, or even within, cell types. To accurately measure DNA replication speeds, we developed single-cell 5-ethynyl-2'-deoxyuridine sequencing to detect nascent replicated DNA. We observed that the DNA replication speed is not constant but increases during S phase of the cell cycle. Using genetic and pharmacological perturbations we were able to alter this acceleration of replication and conclude that DNA damage inflicted by the process of transcription limits the speed of replication during early S phase. In late S phase, during which less-transcribed regions replicate, replication accelerates and approaches its maximum speed.

Replication initiation sites and zones in the mammalian genome: Where are they located and how are they defined?

Eukaryotic DNA replication is a tightly controlled process that occurs in two main steps, i.e., licensing and firing, which take place in the G1 and S phases of the cell cycle, respectively. In Saccharomyces cerevisiae, the budding yeast, replication origins contain consensus sequences that are recognized and bound by the licensing factor Orc1-6, which then recruits the replicative Mcm2-7 helicase. By contrast, mammalian initiation sites lack such consensus sequences, and the mammalian ORC does not exhibit sequence specificity. Studies performed over the past decades have identified replication initiation sites in the mammalian genome using sequencing-based assays, raising the question of whether replication initiation occurs at confined sites or in broad zones across the genome. Although recent reports have shown that the licensed MCMs in mammalian cells are broadly distributed, suggesting that ORC-dependent licensing may not determine the initiation sites/zones, they are predominantly located upstream of actively transcribed genes. This review compares the mechanism of replication initiation in yeast and mammalian cells, summarizes the sequencing-based technologies used for the identification of initiation sites/zones, and proposes a possible mechanism of initiation-site/zone selection in mammalian cells. Future directions and challenges in this field are also discussed.

Monitoring and quantifying replication fork dynamics with high-throughput methods

Before each cell division, eukaryotic cells must replicate their chromosomes to ensure the accurate transmission of genetic information. Chromosome replication involves more than just DNA duplication; it also includes chromatin assembly, inheritance of epigenetic marks, and faithful resumption of all genomic functions after replication. Recent progress in quantitative technologies has revolutionized our understanding of the complexity and dynamics of DNA replication forks at both molecular and genomic scales. Here, we highlight the pivotal role of these novel methods in uncovering the principles and mechanisms of chromosome replication. These technologies have illuminated the regulation of genome replication programs, quantified the impact of DNA replication on genomic mutations and evolutionary processes, and elucidated the mechanisms of replication-coupled chromatin assembly and epigenome maintenance.

Single-Macromolecule Studies of Eukaryotic Genomic Maintenance

Genomes are self-organized and self-maintained as long, complex macromolecules of chromatin. The inherent heterogeneity, stochasticity, phase separation, and chromatin dynamics of genome operation make it challenging to study genomes using ensemble methods. Various single-molecule force-, fluorescent-, and sequencing-based techniques rooted in different disciplines have been developed to fill critical gaps in the capabilities of bulk measurements, each providing unique, otherwise inaccessible, insights into the structure and maintenance of the genome. Capable of capturing molecular-level details about the organization, conformational changes, and packaging of genetic material, as well as processive and stochastic movements of maintenance factors, a single-molecule toolbox provides an excellent opportunity for collaborative research to understand how genetic material functions in health and malfunctions in disease. In this review, we discuss novel insights brought to genomic sciences by single-molecule techniques and their potential to continue to revolutionize the field-one molecule at a time.

Detection of ribonucleotides embedded in DNA by Nanopore sequencing

Ribonucleotides represent the most common non-canonical nucleotides found in eukaryotic genomes. The sources of chromosome-embedded ribonucleotides and the mechanisms by which unrepaired rNMPs trigger genome instability and human pathologies are not fully understood. The available sequencing technologies only allow to indirectly deduce the genomic location of rNMPs. Oxford Nanopore Technologies (ONT) may overcome such limitation, revealing the sites of rNMPs incorporation in genomic DNA directly from raw sequencing signals. We synthesized two types of DNA molecules containing rNMPs at known or random positions and we developed data analysis pipelines for DNA-embedded ribonucleotides detection by ONT. We report that ONT can identify all four ribonucleotides incorporated in DNA by capturing rNMPs-specific alterations in nucleotide alignment features, current intensity, and dwell time. We propose that ONT may be successfully employed to directly map rNMPs in genomic DNA and we suggest a strategy to build an ad hoc basecaller to analyse native genomes.

Integrating high-throughput analysis to create an atlas of replication origins in Trypanosoma cruzi in the context of genome structure and variability

Trypanosoma cruzi is the etiologic agent of the most prevalent human parasitic disease in Latin America, Chagas disease. Its genome is rich in multigenic families that code for virulent antigens and are present in the rapidly evolving genomic compartment named Disruptive. DNA replication is a meticulous biological process in which flaws can generate mutations and changes in chromosomal and gene copy numbers. Here, integrating high-throughput and single-molecule analyses, we were able to identify Predominant, Flexible, and Dormant Orc1Cdc6-dependent origins as well as Orc1Cdc6-independent origins. Orc1Cdc6-dependent origins were found in multigenic family loci, while independent origins were found in the Core compartment that contains conserved and hypothetical protein-coding genes, in addition to multigenic families. In addition, we found that Orc1Cdc6 density is related to the firing of origins and that Orc1Cdc6-binding sites within fired origins are depleted of a specific class of nucleosomes that we previously categorized as dynamic. Together, these data suggest that Orc1Cdc6-dependent origins may contribute to the rapid evolution of the Disruptive compartment and, therefore, to the success of T. cruzi infection and that the local epigenome landscape is also involved in this process.IMPORTANCETrypanosoma cruzi, responsible for Chagas disease, affects millions globally, particularly in Latin America. Lack of vaccine or treatment underscores the need for research. Parasite's genome, with virulent antigen-coding multigenic families, resides in the rapidly evolving Disruptive compartment. Study sheds light on the parasite's dynamic DNA replication, discussing the evolution of the Disruptive compartment. Therefore, the findings represent a significant stride in comprehending T. cruzi's biology and the molecular bases that contribute to the success of infection caused by this parasite.

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.

Combing of picogram level DNA equivalent to genomic DNA present in single human cell by self propelled droplet motion over a stable gradient surface

DNA combing is a powerful technique for studying replication profile, fork-directionality and fork velocity. At present, there is requirement of a methodology to comb DNA present in a single human cell for studying replication dynamics at early embryonic stage. In our study, a surface having dual characteristics i.e., affinity towards negatively charged single DNA molecules and a hydrophobic gradient for self propelled droplet motion of combing solution was developed. The surface was made by coating of TCOS (trichloro-octylsilane) by vapor diffusion on APTES (Aminopropyl-triethoxysilane) coated glass slides. A gradient surface having high deposition efficiency (DE) was developed on which 5 picogram DNA equivalent to genomic DNA present in one single human cell can be combed. The gradient surface was thermostable in nature having the ability to sustain boiling temperature for two hours and sustain anisotropy in 70 % ethanol for 80 h. Applicability for multiple runs was enhanced such that the surface can be used for 13-14 times. Factors associated with gradient surface are unidirectional movement of combing solution droplet over the gradient surface for combing straight DNA molecules and a longer gradient surface of more than 1 cm such that long size DNA molecules can be combed. Ellipsometry and contact angle hysteresis confirmed the presence of hydrophobic gradient. XPS (X-ray photoelectron spectroscopy) and FTIR (Fourier Transform Infrared Spectroscopy) confirmed the presence of characteristic affinity towards negatively charged DNA molecules on the gradient surface. Combing solution was optimized for increasing deposition efficiency and for increasing the applicability of gradient surface for multiple runs. High temperature of combing solution was found to increase Deposition Efficiency. Combing solution was also optimized for combing single DNA molecules over the gradient surface. Single DNA molecules were combed by reducing pH and lowering concentration of triton-X in the combing solution. Dye: bp ratio was optimized for high fluorescent intensity and low surface background.

Adapting Nanopore Sequencing Basecalling Models for Modification Detection via Incremental Learning and Anomaly Detection

We leverage machine learning approaches to adapt nanopore sequencing basecallers for nucleotide modification detection. We first apply the incremental learning technique to improve the basecalling of modification-rich sequences, which are usually of high biological interests. With sequence backbones resolved, we further run anomaly detection on individual nucleotides to determine their modification status. By this means, our pipeline promises the single-molecule, single-nucleotide and sequence context-free detection of modifications. We benchmark the pipeline using control oligos, further apply it in the basecalling of densely-modified yeast tRNAs and E.coli genomic DNAs, the cross-species detection of N6-methyladenosine (m6A) in mammalian mRNAs, and the simultaneous detection of N1-methyladenosine (m1A) and m6A in human mRNAs. Our IL-AD workflow is available at: https://github.com/wangziyuan66/IL-AD.

Single molecule MATAC-seq reveals key determinants of DNA replication origin efficiency

Stochastic origin activation gives rise to significant cell-to-cell variability in the pattern of genome replication. The molecular basis for heterogeneity in efficiency and timing of individual origins is a long-standing question. Here, we developed Methylation Accessibility of TArgeted Chromatin domain Sequencing (MATAC-Seq) to determine single-molecule chromatin accessibility of four specific genomic loci. MATAC-Seq relies on preferential modification of accessible DNA by methyltransferases combined with Nanopore-Sequencing for direct readout of methylated DNA-bases. Applying MATAC-Seq to selected early-efficient and late-inefficient yeast replication origins revealed large heterogeneity of chromatin states. Disruption of INO80 or ISW2 chromatin remodeling complexes leads to changes at individual nucleosomal positions that correlate with changes in their replication efficiency. We found a chromatin state with an accessible nucleosome-free region in combination with well-positioned +1 and +2 nucleosomes as a strong predictor for efficient origin activation. Thus, MATAC-Seq identifies the large spectrum of alternative chromatin states that co-exist on a given locus previously masked in population-based experiments and provides a mechanistic basis for origin activation heterogeneity during eukaryotic DNA replication. Consequently, our single-molecule chromatin accessibility assay will be ideal to define single-molecule heterogeneity across many fundamental biological processes such as transcription, replication, or DNA repair in vitro and ex vivo.

Where and when to start: Regulating DNA replication origin activity in eukaryotic genomes

In eukaryotic genomes, hundreds to thousands of potential start sites of DNA replication named origins are dispersed across each of the linear chromosomes. During S-phase, only a subset of origins is selected in a stochastic manner to assemble bidirectional replication forks and initiate DNA synthesis. Despite substantial progress in our understanding of this complex process, a comprehensive 'identity code' that defines origins based on specific nucleotide sequences, DNA structural features, the local chromatin environment, or 3D genome architecture is still missing. In this article, we review the genetic and epigenetic features of replication origins in yeast and metazoan chromosomes and highlight recent insights into how this flexibility in origin usage contributes to nuclear organization, cell growth, differentiation, and genome stability.

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.

Genomics in the long-read sequencing era

Long-read sequencing (LRS) technologies have provided extremely powerful tools to explore genomes. While in the early years these methods suffered technical limitations, they have recently made significant progress in terms of read length, throughput, and accuracy and bioinformatics tools have strongly improved. Here, we aim to review the current status of LRS technologies, the development of novel methods, and the impact on genomics research. We will explore the most impactful recent findings made possible by these technologies focusing on high-resolution sequencing of genomes and transcriptomes and the direct detection of DNA and RNA modifications. We will also discuss how LRS methods promise a more comprehensive understanding of human genetic variation, transcriptomics, and epigenetics for the coming years.

Single-copy locus proteomics of early- and late-firing DNA replication origins identifies a role of Ask1/DASH complex in replication timing control

The chromatin environment at origins of replication is thought to influence DNA replication initiation in eukaryotic genomes. However, it remains unclear how and which chromatin features control the firing of early-efficient (EE) or late-inefficient (LI) origins. Here, we use site-specific recombination and single-locus chromatin isolation to purify EE and LI replication origins in Saccharomyces cerevisiae. Using mass spectrometry, we define the protein composition of native chromatin regions surrounding the EE and LI replication start sites. In addition to known origin interactors, we find the microtubule-binding Ask1/DASH complex as an origin-regulating factor. Strikingly, tethering of Ask1 to individual origin sites advances replication timing (RT) of the targeted chromosomal domain. Targeted degradation of Ask1 globally changes RT of a subset of origins, which can be reproduced by inhibiting microtubule dynamics. Thus, our findings mechanistically connect RT and chromosomal organization via Ask1/DASH with the microtubule cytoskeleton.

Origins of DNA replication in eukaryotes

Errors occurring during DNA replication can result in inaccurate replication, incomplete replication, or re-replication, resulting in genome instability that can lead to diseases such as cancer or disorders such as autism. A great deal of progress has been made toward understanding the entire process of DNA replication in eukaryotes, including the mechanism of initiation and its control. This review focuses on the current understanding of how the origin recognition complex (ORC) contributes to determining the location of replication initiation in the multiple chromosomes within eukaryotic cells, as well as methods for mapping the location and temporal patterning of DNA replication. Origin specification and configuration vary substantially between eukaryotic species and in some cases co-evolved with gene-silencing mechanisms. We discuss the possibility that centromeres and origins of DNA replication were originally derived from a common element and later separated during evolution.

Genome-wide measurement of DNA replication fork directionality and quantification of DNA replication initiation and termination with Okazaki fragment sequencing

Studying the dynamics of genome replication in mammalian cells has been historically challenging. To reveal the location of replication initiation and termination in the human genome, we developed Okazaki fragment sequencing (OK-seq), a quantitative approach based on the isolation and strand-specific sequencing of Okazaki fragments, the lagging strand replication intermediates. OK-seq quantitates the proportion of leftward- and rightward-oriented forks at every genomic locus and reveals the location and efficiency of replication initiation and termination events. Here we provide the detailed experimental procedures for performing OK-seq in unperturbed cultured human cells and budding yeast and the bioinformatics pipelines for data processing and computation of replication fork directionality. Furthermore, we present the analytical approach based on a hidden Markov model, which allows automated detection of ascending, descending and flat replication fork directionality segments revealing the zones of replication initiation, termination and unidirectional fork movement across the entire genome. These tools are essential for the accurate interpretation of human and yeast replication programs. The experiments and the data processing can be accomplished within six days. Besides revealing the genome replication program in fine detail, OK-seq has been instrumental in numerous studies unravelling mechanisms of genome stability, epigenome maintenance and genome evolution.

Unscheduled DNA replication in G1 causes genome instability and damage signatures indicative of replication collisions

DNA replicates once per cell cycle. Interfering with the regulation of DNA replication initiation generates genome instability through over-replication and has been linked to early stages of cancer development. Here, we engineer genetic systems in budding yeast to induce unscheduled replication in a G1-like cell cycle state. Unscheduled G1 replication initiates at canonical S-phase origins. We quantifiy the composition of replisomes in G1- and S-phase and identified firing factors, polymerase α, and histone supply as factors that limit replication outside S-phase. G1 replication per se does not trigger cellular checkpoints. Subsequent replication during S-phase, however, results in over-replication and leads to chromosome breaks and chromosome-wide, strand-biased occurrence of RPA-bound single-stranded DNA, indicating head-to-tail replication collisions as a key mechanism generating genome instability upon G1 replication. Low-level, sporadic induction of G1 replication induces an identical response, indicating findings from synthetic systems are applicable to naturally occurring scenarios of unscheduled replication initiation.

Targeted profiling of human extrachromosomal DNA by CRISPR-CATCH

Extrachromosomal DNA (ecDNA) is a common mode of oncogene amplification but is challenging to analyze. Here, we adapt CRISPR-CATCH, in vitro CRISPR-Cas9 treatment and pulsed field gel electrophoresis of agarose-entrapped genomic DNA, previously developed for bacterial chromosome segments, to isolate megabase-sized human ecDNAs. We demonstrate strong enrichment of ecDNA molecules containing EGFR, FGFR2 and MYC from human cancer cells and NRAS ecDNA from human metastatic melanoma with acquired therapeutic resistance. Targeted enrichment of ecDNA versus chromosomal DNA enabled phasing of genetic variants, identified the presence of an EGFRvIII mutation exclusively on ecDNAs and supported an excision model of ecDNA genesis in a glioblastoma model. CRISPR-CATCH followed by nanopore sequencing enabled single-molecule ecDNA methylation profiling and revealed hypomethylation of the EGFR promoter on ecDNAs. We distinguished heterogeneous ecDNA species within the same sample by size and sequence with base-pair resolution and discovered functionally specialized ecDNAs that amplify select enhancers or oncogene-coding sequences.

Modification mapping by nanopore sequencing

Next generation sequencing (NGS) has provided biologists with an unprecedented view into biological processes and their regulation over the past 2 decades, fueling a wave of development of high throughput methods based on short read DNA and RNA sequencing. For nucleic acid modifications, NGS has been coupled with immunoprecipitation, chemical treatment, enzymatic treatment, and/or the use of reverse transcriptase enzymes with fortuitous activities to enrich for and to identify covalent modifications of RNA and DNA. However, the majority of nucleic acid modifications lack commercial monoclonal antibodies, and mapping techniques that rely on chemical or enzymatic treatments to manipulate modification signatures add additional technical complexities to library preparation. Moreover, such approaches tend to be specific to a single class of RNA or DNA modification, and generate only indirect readouts of modification status. Third generation sequencing technologies such as the commercially available "long read" platforms from Pacific Biosciences and Oxford Nanopore Technologies are an attractive alternative for high throughput detection of nucleic acid modifications. While the former can indirectly sense modified nucleotides through changes in the kinetics of reverse transcription reactions, nanopore sequencing can in principle directly detect any nucleic acid modification that produces a signal distortion as the nucleic acid passes through a nanopore sensor embedded within a charged membrane. To date, more than a dozen endogenous DNA and RNA modifications have been interrogated by nanopore sequencing, as well as a number of synthetic nucleic acid modifications used in metabolic labeling, structure probing, and other emerging applications. This review is intended to introduce the reader to nanopore sequencing and key principles underlying its use in direct detection of nucleic acid modifications in unamplified DNA or RNA samples, and outline current approaches for detecting and quantifying nucleic acid modifications by nanopore sequencing. As this technology matures, we anticipate advances in both sequencing chemistry and analysis methods will lead to rapid improvements in the identification and quantification of these epigenetic marks.

Replicon-seq: seeing is believing

Claussin et al. introduce Replicon-seq, a new genome-wide DNA sequencing technology that monitors the progression of individual replisomes at high resolution in vivo.

Replication dependent and independent mechanisms of GAA repeat instability

Trinucleotide repeat instability is a driver of human disease. Large expansions of (GAA)n repeats in the first intron of the FXN gene are the cause Friedreich's ataxia (FRDA), a progressive degenerative disorder which cannot yet be prevented or treated. (GAA)n repeat instability arises during both replication-dependent processes, such as cell division and intergenerational transmission, as well as in terminally differentiated somatic tissues. Here, we provide a brief historical overview on the discovery of (GAA)n repeat expansions and their association to FRDA, followed by recent advances in the identification of triplex H-DNA formation and replication fork stalling. The main body of this review focuses on the last decade of progress in understanding the mechanism of (GAA)n repeat instability during DNA replication and/or DNA repair. We propose that the discovery of additional mechanisms of (GAA)n repeat instability can be achieved via both comparative approaches to other repeat expansion diseases and genome-wide association studies. Finally, we discuss the advances towards FRDA prevention or amelioration that specifically target (GAA)n repeat expansions.

Genome-wide mapping of individual replication fork velocities using nanopore sequencing

Little is known about replication fork velocity variations along eukaryotic genomes, since reference techniques to determine fork speed either provide no sequence information or suffer from low throughput. Here we present NanoForkSpeed, a nanopore sequencing-based method to map and extract the velocity of individual forks detected as tracks of the thymidine analogue bromodeoxyuridine incorporated during a brief pulse-labelling of asynchronously growing cells. NanoForkSpeed retrieves previous Saccharomyces cerevisiae mean fork speed estimates (≈2 kb/min) in the BT1 strain exhibiting highly efficient bromodeoxyuridine incorporation and wild-type growth, and precisely quantifies speed changes in cells with altered replisome progression or exposed to hydroxyurea. The positioning of >125,000 fork velocities provides a genome-wide map of fork progression based on individual fork rates, showing a uniform fork speed across yeast chromosomes except for a marked slowdown at known pausing sites.

Theulot et al. introduce NanoForkSpeed, a nanopore sequencing-based method to map individual replication fork velocities on entire genomes. NFS shows that fork speed is uniform across yeast chromosomes except for a marked slowdown at pausing sites.

DeepRepeat: direct quantification of short tandem repeats on signal data from nanopore sequencing

Despite recent improvements in basecalling accuracy, nanopore sequencing still has higher error rates on short-tandem repeats (STRs). Instead of using basecalled reads, we developed DeepRepeat which converts ionic current signals into red-green-blue channels, thus transforming the repeat detection problem into an image recognition problem. DeepRepeat identifies and accurately quantifies telomeric repeats in the CHM13 cell line and achieves higher accuracy in quantifying repeats in long STRs than competing methods. We also evaluate DeepRepeat on genome-wide or candidate region datasets from seven different sources. In summary, DeepRepeat enables accurate quantification of long STRs and complements existing methods relying on basecalled reads.

Purification of mammalian telomeric DNA for single-molecule analysis

Here we provide a detailed protocol for the enrichment of telomeric repeats from mouse and human cells. The procedure consists of two successive rounds of digestion with frequently cutting restriction enzymes followed by size fractionation. Around 2 mg of genomic DNA is required, and the procedure lasts 5-6 d and yields preparations enriched >800-fold in telomeres. The purified material is suitable for single-molecule analysis of telomere structure, visualizing telomere replication and recombination intermediates by electron microscopy or performing molecular combing at telomeric repeats. No special skills are required for the enrichment procedure, while some assistance is needed in harvesting a large number of plates in a timely fashion at the beginning of the procedure. A smaller-scale version of the protocol that involves one round of digestion and purification requires 200 µg of DNA and enriches telomeres ~50-fold in 4 d is also provided. The latter can be combined with specific labeling for single-molecule analysis of replicating DNA or for long-read sequencing analysis of telomeric repeats. The procedure described here can be adapted to the enrichment of other repetitive elements, based on the use of restriction enzymes that do not cut into the repeat of interest.

Mapping replication forks, one replicon at a time

Claussin et al. (2022) present an elegant approach to replication fork mapping that combines single-molecule resolution with genome-wide coverage to provide unprecedented insight into the robust nature of DNA replication.

Single-molecule mapping of replisome progression

Fundamental aspects of DNA replication, such as the anatomy of replication stall sites, how replisomes are influenced by gene transcription, and whether the progression of sister replisomes is coordinated, are poorly understood. Available techniques do not allow the precise mapping of the positions of individual replisomes on chromatin. We have developed a method called Replicon-seq that entails the excision of full-length replicons by controlled nuclease cleavage at replication forks. Replicons are sequenced using Nanopore, which provides a single-molecule readout of long DNA. Using Replicon-seq, we found that sister replisomes function autonomously and yet progress through chromatin with remarkable consistency. Replication forks that encounter obstacles pause for a short duration but rapidly resume synthesis. The helicase Rrm3 plays a critical role both in mitigating the effect of protein barriers and with facilitating efficient termination. Replicon-seq provides a high-resolution means of defining how individual replisomes move across the genome.

FORK-seq: Single-Molecule Profiling of DNA Replication

Most genome replication mapping methods profile cell populations, masking cell-to-cell heterogeneity. Here, we describe FORK-seq, a nanopore sequencing method to map replication of single DNA molecules at 200 nucleotide resolution using a nanopore current interpretation tool allowing the quantification of BrdU incorporation. Along pulse-chased replication intermediates from Saccharomyces cerevisiae, we can orient replication tracks and reproduce population-based replication directionality profiles. Additionally, we can map individual initiation and termination events. Thus, FORK-seq reveals the full extent of cell-to-cell heterogeneity in DNA replication.

Intellectual disability genomics: current state, pitfalls and future challenges

Intellectual disability (ID) can be caused by non-genetic and genetic factors, the latter being responsible for more than 1700 ID-related disorders. The broad ID phenotypic and genetic heterogeneity, as well as the difficulty in the establishment of the inheritance pattern, often result in a delay in the diagnosis. It has become apparent that massive parallel sequencing can overcome these difficulties. In this review we address: (i) ID genetic aetiology, (ii) clinical/medical settings testing, (iii) massive parallel sequencing, (iv) variant filtering and prioritization, (v) variant classification guidelines and functional studies, and (vi) ID diagnostic yield. Furthermore, the need for a constant update of the methodologies and functional tests, is essential. Thus, international collaborations, to gather expertise, data and resources through multidisciplinary contributions, are fundamental to keep track of the fast progress in ID gene discovery.

A blood drop through the pore: nanopore sequencing in hematology

The development of new sequencing platforms, technologies, and bioinformatics tools in the past decade fostered key discoveries in human genomics. Among the most recent sequencing technologies, nanopore sequencing (NS) has caught the interest of researchers for its intriguing potential and flexibility. This up-to-date review highlights the recent application of NS in the hematology field, focusing on progress and challenges of the technological approaches employed for the identification of pathologic alterations. The molecular and analytic pipelines developed for the analysis of the whole-genome, target regions, and transcriptomics provide a proof of evidence of the unparalleled amount of information that could be retrieved by an innovative approach based on long-read sequencing.

Application of third-generation sequencing in cancer research

In the past several years, nanopore sequencing technology from Oxford Nanopore Technologies (ONT) and single-molecule real-time (SMRT) sequencing technology from Pacific BioSciences (PacBio) have become available to researchers and are currently being tested for cancer research. These methods offer many advantages over most widely used high-throughput short-read sequencing approaches and allow the comprehensive analysis of transcriptomes by identifying full-length splice isoforms and several other posttranscriptional events. In addition, these platforms enable structural variation characterization at a previously unparalleled resolution and direct detection of epigenetic marks in native DNA and RNA. Here, we present a comprehensive summary of important applications of these technologies in cancer research, including the identification of complex structure variants, alternatively spliced isoforms, fusion transcript events, and exogenous RNA. Furthermore, we discuss the impact of the newly developed nanopore direct RNA sequencing (RNA-Seq) approach in advancing epitranscriptome research in cancer. Although the unique challenges still present for these new single-molecule long-read methods, they will unravel many aspects of cancer genome complexity in unprecedented ways and present an encouraging outlook for continued application in an increasing number of different cancer research settings.

Variation and Evolution of Human Centromeres: A Field Guide and Perspective

We are entering a new era in genomics where entire centromeric regions are accurately represented in human reference assemblies. Access to these high-resolution maps will enable new surveys of sequence and epigenetic variation in the population and offer new insight into satellite array genomics and centromere function. Here, we focus on the sequence organization and evolution of alpha satellites, which are credited as the genetic and genomic definition of human centromeres due to their interaction with inner kinetochore proteins and their importance in the development of human artificial chromosome assays. We provide an overview of alpha satellite repeat structure and array organization in the context of these high-quality reference data sets; discuss the emergence of variation-based surveys; and provide perspective on the role of this new source of genetic and epigenetic variation in the context of chromosome biology, genome instability, and human disease.

Towards inferring nanopore sequencing ionic currents from nucleotide chemical structures

The characteristic ionic currents of nucleotide kmers are commonly used in analyzing nanopore sequencing readouts. We present a graph convolutional network-based deep learning framework for predicting kmer characteristic ionic currents from corresponding chemical structures. We show such a framework can generalize the chemical information of the 5-methyl group from thymine to cytosine by correctly predicting 5-methylcytosine-containing DNA 6mers, thus shedding light on the de novo detection of nucleotide modifications.

Nanopore sequencing technology, bioinformatics and applications

Rapid advances in nanopore technologies for sequencing single long DNA and RNA molecules have led to substantial improvements in accuracy, read length and throughput. These breakthroughs have required extensive development of experimental and bioinformatics methods to fully exploit nanopore long reads for investigations of genomes, transcriptomes, epigenomes and epitranscriptomes. Nanopore sequencing is being applied in genome assembly, full-length transcript detection and base modification detection and in more specialized areas, such as rapid clinical diagnoses and outbreak surveillance. Many opportunities remain for improving data quality and analytical approaches through the development of new nanopores, base-calling methods and experimental protocols tailored to particular applications.

Barapost: Binning of Nucleotide Sequences According to Taxonomic Annotation

Contemporary sequencing technologies, Oxford Nanopore in particular, provide a way to sequence multiple samples during single run using molecular barcodes. Specific circumstances, however, can make barcoding undesirable or unaffordable. Here, we introduce Barapost: a command-line toolkit that demultiplexes long nucleotide sequences relying on taxonomic annotation instead of barcoding.

Nanopore sequencing and de novo assembly of a misidentified Camelpox vaccine reveals putative epigenetic modifications and alternate protein signal peptides

DNA viruses can exploit host cellular epigenetic processes to their advantage; however, the epigenome status of most DNA viruses remains undetermined. Third generation sequencing technologies allow for the identification of modified nucleotides from sequencing experiments without specialized sample preparation, permitting the detection of non-canonical epigenetic modifications that may distinguish viral nucleic acid from that of their host, thus identifying attractive targets for advanced therapeutics and diagnostics. We present a novel nanopore de novo assembly pipeline used to assemble a misidentified Camelpox vaccine. Two confirmed deletions of this vaccine strain in comparison to the closely related Vaccinia virus strain modified vaccinia Ankara make it one of the smallest non-vector derived orthopoxvirus genomes to be reported. Annotation of the assembly revealed a previously unreported signal peptide at the start of protein A38 and several predicted signal peptides that were found to differ from those previously described. Putative epigenetic modifications around various motifs have been identified and the assembly confirmed previous work showing the vaccine genome to most closely resemble that of Vaccinia virus strain Modified Vaccinia Ankara. The pipeline may be used for other DNA viruses, increasing the understanding of DNA virus evolution, virulence, host preference, and epigenomics.

Towards population-scale long-read sequencing

Long-read sequencing technologies have now reached a level of accuracy and yield that allows their application to variant detection at a scale of tens to thousands of samples. Concomitant with the development of new computational tools, the first population-scale studies involving long-read sequencing have emerged over the past 2 years and, given the continuous advancement of the field, many more are likely to follow. In this Review, we survey recent developments in population-scale long-read sequencing, highlight potential challenges of a scaled-up approach and provide guidance regarding experimental design. We provide an overview of current long-read sequencing platforms, variant calling methodologies and approaches for de novo assemblies and reference-based mapping approaches. Furthermore, we summarize strategies for variant validation, genotyping and predicting functional impact and emphasize challenges remaining in achieving long-read sequencing at a population scale.

Promoters are key organizers of the duplication of vertebrate genomes

In vertebrates, single cell analyses of replication timing patterns brought to light a very well controlled program suggesting a tight regulation on initiation sites. Mapping of replication origins with different methods has revealed discrete preferential sites, enriched in promoters and potential G-quadruplex motifs, which can aggregate into initiation zones spanning several tens of kilobases (kb). Another characteristic of replication origins is a nucleosome-free region (NFR). A modified yeast strain containing a humanized origin recognition complex (ORC) fires new origins at NFRs revealing their regulatory role. In cooperation with NFRs, the histone variant H2A.Z facilitates ORC loading through di-methylation of lysine 20 of histone H4. Recent studies using genome editing methods show that efficient initiation sites associated with transcriptional activity can synergize over several tens of kb by establishing physical contacts and lead to the formation of early domains of DNA replication demonstrating a co-regulation between replication initiation and transcription.

DNA molecular combing-based replication fork directionality profiling

The replication strategy of metazoan genomes is still unclear, mainly because definitive maps of replication origins are missing. High-throughput methods are based on population average and thus may exclusively identify efficient initiation sites, whereas inefficient origins go undetected. Single-molecule analyses of specific loci can detect both common and rare initiation events along the targeted regions. However, these usually concentrate on positioning individual events, which only gives an overview of the replication dynamics. Here, we computed the replication fork directionality (RFD) profiles of two large genes in different transcriptional states in chicken DT40 cells, namely untranscribed and transcribed DMD and CCSER1 expressed at WT levels or overexpressed, by aggregating hundreds of oriented replication tracks detected on individual DNA fibres stretched by molecular combing. These profiles reconstituted RFD domains composed of zones of initiation flanking a zone of termination originally observed in mammalian genomes and were highly consistent with independent population-averaging profiles generated by Okazaki fragment sequencing. Importantly, we demonstrate that inefficient origins do not appear as detectable RFD shifts, explaining why dispersed initiation has remained invisible to population-based assays. Our method can both generate quantitative profiles and identify discrete events, thereby constituting a comprehensive approach to study metazoan genome replication.

Single-Molecule Techniques to Study Chromatin

Besides the basic organization in nucleosome core particles (NCPs), eukaryotic chromatin is further packed through interactions with numerous protein complexes including transcription factors, chromatin remodeling and modifying enzymes. This nucleoprotein complex provides the template for many important biological processes, such as DNA replication, transcription, and DNA repair. Thus, to understand the molecular basis of these DNA transactions, it is critical to define individual changes of the chromatin structure at precise genomic regions where these machineries assemble and drive biological reactions. Single-molecule approaches provide the only possible solution to overcome the heterogenous nature of chromatin and monitor the behavior of individual chromatin transactions in real-time. In this review, we will give an overview of currently available single-molecule methods to obtain mechanistic insights into nucleosome positioning, histone modifications and DNA replication and transcription analysis-previously unattainable with population-based assays.

Genome-wide mapping of human DNA replication by optical replication mapping supports a stochastic model of eukaryotic replication

The heterogeneous nature of eukaryotic replication kinetics and the low efficiency of individual initiation sites make mapping the location and timing of replication initiation in human cells difficult. To address this challenge, we have developed optical replication mapping (ORM), a high-throughput single-molecule approach, and used it to map early-initiation events in human cells. The single-molecule nature of our data and a total of >2,500-fold coverage of the human genome on 27 million fibers averaging ∼300 kb in length allow us to identify initiation sites and their firing probability with high confidence. We find that the distribution of human replication initiation is consistent with inefficient, stochastic activation of heterogeneously distributed potential initiation complexes enriched in accessible chromatin. These observations are consistent with stochastic models of initiation-timing regulation and suggest that stochastic regulation of replication kinetics is a fundamental feature of eukaryotic replication, conserved from yeast to humans.

DNAscent v2: detecting replication forks in nanopore sequencing data with deep learning

Background

Measuring DNA replication dynamics with high throughput and single-molecule resolution is critical for understanding both the basic biology behind how cells replicate their DNA and how DNA replication can be used as a therapeutic target for diseases like cancer. In recent years, the detection of base analogues in Oxford Nanopore Technologies (ONT) sequencing reads has become a promising new method to supersede existing single-molecule methods such as DNA fibre analysis: ONT sequencing yields long reads with high throughput, and sequenced molecules can be mapped to the genome using standard sequence alignment software.

Results

This paper introduces DNAscent v2, software that uses a residual neural network to achieve fast, accurate detection of the thymidine analogue BrdU with single-nucleotide resolution. DNAscent v2 also comes equipped with an autoencoder that interprets the pattern of BrdU incorporation on each ONT-sequenced molecule into replication fork direction to call the location of replication origins termination sites. DNAscent v2 surpasses previous versions of DNAscent in BrdU calling accuracy, origin calling accuracy, speed, and versatility across different experimental protocols. Unlike NanoMod, DNAscent v2 positively identifies BrdU without the need for sequencing unmodified DNA. Unlike RepNano, DNAscent v2 calls BrdU with single-nucleotide resolution and detects more origins than RepNano from the same sequencing data. DNAscent v2 is open-source and available at https://github.com/MBoemo/DNAscent.

Conclusions

This paper shows that DNAscent v2 is the new state-of-the-art in the high-throughput, single-molecule detection of replication fork dynamics. These improvements in DNAscent v2 mark an important step towards measuring DNA replication dynamics in large genomes with single-molecule resolution. Looking forward, the increase in accuracy in single-nucleotide resolution BrdU calls will also allow DNAscent v2 to branch out into other areas of genome stability research, particularly the detection of DNA repair.

Supplementary Information

The online version contains supplementary material available at (10.1186/s12864-021-07736-6).

Approaching Protein Barriers: Emerging Mechanisms of Replication Pausing in Eukaryotes

During nuclear DNA replication multiprotein replisome machines have to jointly traverse and duplicate the total length of each chromosome during each cell cycle. At certain genomic locations replisomes encounter tight DNA-protein complexes and slow down. This fork pausing is an active process involving recognition of a protein barrier by the approaching replisome via an evolutionarily conserved Fork Pausing/Protection Complex (FPC). Action of the FPC protects forks from collapse at both programmed and accidental protein barriers, thus promoting genome integrity. In addition, FPC stimulates the DNA replication checkpoint and regulates topological transitions near the replication fork. Eukaryotic cells have been proposed to employ physiological programmed fork pausing for various purposes, such as maintaining copy number at repetitive loci, precluding replication-transcription encounters, regulating kinetochore assembly, or controlling gene conversion events during mating-type switching. Here we review the growing number of approaches used to study replication pausing in vivo and in vitro as well as the characterization of additional factors recently reported to modulate fork pausing in different systems. Specifically, we focus on the positive role of topoisomerases in fork pausing. We describe a model where replisome progression is inherently cautious, which ensures general preservation of fork stability and genome integrity but can also carry out specialized functions at certain loci. Furthermore, we highlight classical and novel outstanding questions in the field and propose venues for addressing them. Given how little is known about replisome pausing at protein barriers in human cells more studies are required to address how conserved these mechanisms are.

An evolutionary model identifies the main evolutionary biases for the evolution of genome-replication profiles

Recent results comparing the temporal program of genome replication of yeast species belonging to the Lachancea clade support the scenario that the evolution of the replication timing program could be mainly driven by correlated acquisition and loss events of active replication origins. Using these results as a benchmark, we develop an evolutionary model defined as birth-death process for replication origins and use it to identify the evolutionary biases that shape the replication timing profiles. Comparing different evolutionary models with data, we find that replication origin birth and death events are mainly driven by two evolutionary pressures, the first imposes that events leading to higher double-stall probability of replication forks are penalized, while the second makes less efficient origins more prone to evolutionary loss. This analysis provides an empirically grounded predictive framework for quantitative evolutionary studies of the replication timing program.

Gaussian mixture model-based unsupervised nucleotide modification number detection using nanopore-sequencing readouts

MOTIVATION

Nucleotide modification status can be decoded from the Oxford Nanopore Technologies nanopore-sequencing ionic current signals. Although various algorithms have been developed for nanopore-sequencing-based modification analysis, more detailed characterizations, such as modification numbers, corresponding signal levels and proportions are still lacking.

RESULTS

We present a framework for the unsupervised determination of the number of nucleotide modifications from nanopore-sequencing readouts. We demonstrate the approach can effectively recapitulate the number of modifications, the corresponding ionic current signal levels, as well as mixing proportions under both DNA and RNA contexts. We further show, by integrating information from multiple detected modification regions, that the modification status of DNA and RNA molecules can be inferred. This method forms a key step of de novo characterization of nucleotide modifications, shedding light on the interpretation of various biological questions.

AVAILABILITY AND IMPLEMENTATION

Modified nanopolish: https://github.com/adbailey4/nanopolish/tree/cigar_output. All other codes used to reproduce the results: https://github.com/hd2326/ModificationNumber.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Genome-wide analysis of DNA replication and DNA double-strand breaks using TrAEL-seq

Faithful replication of the entire genome requires replication forks to copy large contiguous tracts of DNA, and sites of persistent replication fork stalling present a major threat to genome stability. Understanding the distribution of sites at which replication forks stall, and the ensuing fork processing events, requires genome-wide methods that profile replication fork position and the formation of recombinogenic DNA ends. Here, we describe Transferase-Activated End Ligation sequencing (TrAEL-seq), a method that captures single-stranded DNA 3' ends genome-wide and with base pair resolution. TrAEL-seq labels both DNA breaks and replication forks, providing genome-wide maps of replication fork progression and fork stalling sites in yeast and mammalian cells. Replication maps are similar to those obtained by Okazaki fragment sequencing; however, TrAEL-seq is performed on asynchronous populations of wild-type cells without incorporation of labels, cell sorting, or biochemical purification of replication intermediates, rendering TrAEL-seq far simpler and more widely applicable than existing replication fork direction profiling methods. The specificity of TrAEL-seq for DNA 3' ends also allows accurate detection of double-strand break sites after the initiation of DNA end resection, which we demonstrate by genome-wide mapping of meiotic double-strand break hotspots in a dmc1Δ mutant that is competent for end resection but not strand invasion. Overall, TrAEL-seq provides a flexible and robust methodology with high sensitivity and resolution for studying DNA replication and repair, which will be of significant use in determining mechanisms of genome instability.

A streamlined clinical metagenomic sequencing protocol for rapid pathogen identification

Metagenomic next-generation sequencing (mNGS) holds promise as a diagnostic tool for unbiased pathogen identification and precision medicine. However, its medical utility depends largely on assay simplicity and reproducibility. In the current study, we aimed to develop a streamlined Illumina and Oxford Nanopore-based DNA/RNA library preparation protocol and rapid data analysis pipeline. The Illumina sequencing-based mNGS method was first developed and evaluated using a set of samples with known aetiology. Its sensitivity for RNA viruses (influenza A, H1N1) was < 6.4 × 10² EID50/mL, and a good correlation between viral loads and mapped reads was observed. Then, the rapid turnaround time of Nanopore sequencing was tested by sequencing influenza A virus and adenoviruses. Furthermore, 11 respiratory swabs or sputum samples pre-tested for a panel of pathogens were analysed, and the pathogens identified by Illumina sequencing showed 81.8% concordance with qPCR results. Additional sequencing of cerebrospinal fluid (CSF) samples from HIV-1-positive patients with meningitis/encephalitis detected HIV-1 RNA and Toxoplasma gondii sequences. In conclusion, we have developed a simplified protocol that realizes efficient metagenomic sequencing of a variety of clinical samples and pathogen identification in a clinically meaningful time frame.

Mapping replication timing domains genome wide in single mammalian cells with single-cell DNA replication sequencing

Replication timing (RT) domains are stable units of chromosome structure that are regulated in the context of development and disease. Conventional genome-wide RT mapping methods require many S-phase cells for either the effective enrichment of replicating DNA through bromodeoxyuridine (BrdU) immunoprecipitation or the determination of copy-number differences during S-phase, which precludes their application to non-abundant cell types and single cells. Here, we provide a simple, cost-effective, and robust protocol for single-cell DNA replication sequencing (scRepli-seq). The scRepli-seq methodology relies on whole-genome amplification (WGA) of genomic DNA (gDNA) from single S-phase cells and next-generation sequencing (NGS)-based determination of copy-number differences that arise between replicated and unreplicated DNA. Haplotype-resolved scRepli-seq, which distinguishes pairs of homologous chromosomes within a single cell, is feasible by using single-nucleotide polymorphism (SNP)/indel information. We also provide computational pipelines for quality control, normalization, and binarization of the scRepli-seq data. The experimental portion of this protocol (before sequencing) takes 3 d.

Detection of base analogs incorporated during DNA replication by nanopore sequencing

DNA synthesis is a fundamental requirement for cell proliferation and DNA repair, but no single method can identify the location, direction and speed of replication forks with high resolution. Mammalian cells have the ability to incorporate thymidine analogs along with the natural A, T, G and C bases during DNA synthesis, which allows for labeling of replicating or repaired DNA. Here, we demonstrate the use of the Oxford Nanopore Technologies MinION to detect 11 different thymidine analogs including CldU, BrdU, IdU as well as EdU alone or coupled to Biotin and other bulky adducts in synthetic DNA templates. We also show that the large adduct Biotin can be distinguished from the smaller analog IdU, which opens the possibility of using analog combinations to identify the location and direction of DNA synthesis. Furthermore, we detect IdU label on single DNA molecules in the genome of mouse pluripotent stem cells and using CRISPR/Cas9-mediated enrichment, determine replication rates using newly synthesized DNA strands in human mitochondrial DNA. We conclude that this novel method, termed Replipore sequencing, has the potential for on target examination of DNA replication in a wide range of biological contexts.

DeepSimulator1.5: a more powerful, quicker and lighter simulator for Nanopore sequencing

MOTIVATION

Nanopore sequencing is one of the leading third-generation sequencing technologies. A number of computational tools have been developed to facilitate the processing and analysis of the Nanopore data. Previously, we have developed DeepSimulator1.0 (DS1.0), which is the first simulator for Nanopore sequencing to produce both the raw electrical signals and the reads. However, although DS1.0 can produce high-quality reads, for some sequences, the divergence between the simulated raw signals and the real signals can be large. Furthermore, the Nanopore sequencing technology has evolved greatly since DS1.0 was released. It is thus necessary to update DS1.0 to accommodate those changes.

RESULTS

We propose DeepSimulator1.5 (DS1.5), all three modules of which have been updated substantially from DS1.0. As for the sequence generator, we updated the sample read length distribution to reflect the newest real reads' features. In terms of the signal generator, which is the core of DeepSimulator, we added one more pore model, the context-independent pore model, which is much faster than the previous context-dependent one. Furthermore, to make the generated signals more similar to the real ones, we added a low-pass filter to post-process the pore model signals. Regarding the basecaller, we added the support for the newest official basecaller, Guppy, which can support both GPU and CPU. In addition, multiple optimizations, related to multiprocessing control, memory and storage management, have been implemented to make DS1.5 a much more amenable and lighter simulator than DS1.0.

AVAILABILITY AND IMPLEMENTATION

The main program and the data are available at https://github.com/lykaust15/DeepSimulator.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.