Biology:Ribose-seq

From HandWiki
Short description: Genetic mapping technique

Ribose-seq is a mapping technique used in genetics research to determine the full profile of embedded ribonucleotides, specifically ribonucleoside monophosphates (rNMPs), in genomic DNA.[1] Embedded ribonucleotides are thought to be the most common alteration to DNA in cells,[2][3] and their presence in genomic DNA can affect genome stability.[4][5] As recent studies have suggested that ribonucleotides in mouse DNA may affect disease pathology,[1] ribonucleotide incorporation in genomic DNA has become an important target of medical genetics research.[1] Ribose-seq allows scientists to determine the precise location and type of ribonucleotides that have been incorporated into eukaryotic or prokaryotic DNA.[1]

The technique exploits the presence of the extra hydroxyl groups (OH) found in the 2’ end of ribonucleotides, which can distort and destabilize DNA. The technique was developed through a collaboration with a group of researchers at the Georgia Institute of Technology including Francesca Storici and Kyung Duk Koh (now at the University of California San Francisco), and Jay Hesselberth at the University of Colorado Anschutz Medical School.[6]

A schematic diagram illustrating the workflow of the ribose-seq technique. Adapted from Koh et al. (2015).

History

Nucleic acids are the essential macromolecules that carry genetic information in all life forms.[7] These biopolymers consist of nucleotide monomers, which are organic molecules that consist of a phosphate group, a nitrogen-containing base, and a five-carbon sugar (ribose in ribonucleotides and deoxyribose in deoxyribonucleotides). While ribonucleotides are typically located in ribonucleic acid (RNA), they can also be found in deoxyribonucleic acid (DNA).

The first discovery of ribonucleoside monophosphates (rNMPs) embedded into DNA was in mitochondrial DNA from mouse and human cells in 1973 by Lawrence Grossman, Robert Watson, and Jerome Vinograd.[8] In 2006, the presence of rNMPs were confirmed in nuclear DNA of Schizosaccharomyces pombe ("fission yeast"), and exploration in other species has continued since. [9][10]

Incorporation of ribonucleotides into DNA

Multiple studies have identified mechanisms by which ribonucleotides can be added to, removed from, or generated from DNA.[4] Primarily, ribonucleotides are incorporated into DNA during the DNA synthesis process.[11] To initiate DNA replication, short RNA primers synthesized on complementary DNA to allow for the subsequent binding and replication action of the primary replication enzymes, DNA polymerases. Typically, these primers are then removed by nuclease enzymes, RNases H[12] or Flap Structure-specific Endonuclease 1 (FEN1).[13] As such, there is a fairly substantial volume of transient rNMPs that are present in DNA in the form of RNA primers and then later removed. After DNA replication is initiated, rNMPs can be fully incorporated into DNA by DNA polymerases and this represents the major mechanism of ribonucleotide incorporation in DNA with over 1 ribonucleotide per 1,000 deoxyribonucleotides incorporated.[14][15][16] The embedded rNMPs are often targeted for removal by the Ribonucleotide Excision Repair (RER) mechanism.[17] However, if RER is not working properly, this can become a more persistent source of embedded rNMPs in DNA.[14] Another mechanism by which rNMPs may arise in DNA is through oxidative stress, which alters the deoxyribose unit in DNA to a ribose.[18] While the RER pathway is the most efficient process for removing rNMPs from DNA, ribonucleotides can also be removed through the 3’-5’ exonuclease activity that is orchestrated by the proofreading activities of DNA polymerase during DNA replication.[19] Additional rNMP removal methods include incisions that are mediated by the isomerase enzyme, topoisomerase I.[20]

Consequences of ribonucleotides in DNA

If not removed and replaced, ribonucleotides in DNA can have structural, chemical, and functional impacts on cellular processes.[3] As first recognized by T.N. Jasihree et al in 1993, the presence of ribonucleotides can alter the helical shape of DNA from B form to A form.[21] This, in addition to a number of other chemical consequences,[22][23] can impact the binding of proteins to DNA and have implications for processes such as chromatin dynamics, DNA replication, transcription, repair, and meiosis.[3] Additionally, some DNA polymerases are unable to bypass embedded rNMPs during DNA replication, which is further exacerbated at sites with longer strands of rNMPs.[24] This can cause replication stress which can potentially interfere with transcription, chromosome segregation, and future replication.[3] Furthermore, certain byproducts or failures in the RER process can also potentially lead to DNA breaks, rearrangements, and genome instability.[25] [26]

Ribonucleotide mapping

By 2015, some of the mechanisms for, consequences of, and approximate rates for the incorporation of rNMPs into DNA were becoming increasingly clear.[1] However, a technique for identifying the location of rNMPs across the genome was still lacking. The growing accessibility of next-generation sequencing technologies[27] largely facilitated progress towards the goal of more precise identification of rNMPs throughout genomic DNA. In 2015, Kyung Duk Koh, Sathya Balachander, Jay R Hesselberth and Francesca Storici presented ribose-seq, a newly developed technique for the systematic profiling of rNMPs in DNA.[1] They applied ribose-seq to the Saccharomyces cerevisiae genome to document the genome-wide identity, frequency, and localization of rNMPs throughout mitochondrial and nuclear DNA. These findings ultimately provided new insight into patterns of rNMP incorporations across the genome. In parallel with ribose-seq, three other techniques were developed to capture sites of rNMPs embedded in DNA: emRibo-seq,[28] HydEn-seq,[29] and Pu-seq.[30] The ribose-seq procedure utilizes a number of unique features to capture and tag ribonucleotides in preparation for high-throughput sequencing.[1] To target the rNMP sites directly, the authors utilized an alkali treatment process which denatures and cleaves the DNA at the rNMP sites. The specificity for cleavage at these sites is driven by the relative sensitivity of the molecule at this region due to the additional hyrdroxyl group present on ribonucleotides. To allow for the stability and enrichment of these small fragments and protection from T5 exonuclease digestion ahead of PCR amplification, single stranded DNA fragments - containing the rNMP of interest and a unique molecular identification number - are self-ligated to form single stranded DNA (ssDNA) circles. In order to facilitate ligation of an rNMP and a DNA end, Koh et al. harness the distinctive ligation action of Arabidopsis thaliana tRNA ligase (AtRNL) to allow for the binding of the 2’-phosphate (rNMP side) to the 5′-phosphate (DNA side). Finally, high-throughput sequencing of the entire generated library from this technique identifies the genomic coordinates of all captured embedded rNMPs. In 2019, the Storici lab presented an updated and further optimized protocol for ribose-seq.[31] The modifications included using shorter molecular barcode adaptors, smaller fragments, altered PCR settings, and adjusted size ranges for cutting and purifying from the gel. Overall, it was suggested that the improved protocol may substantially increase the efficiency of the method for applications to large genomes.

Methods

Preparation

Genomic DNA extraction

To extract genomic DNA from cells of interest, cells must first be lysed using mechanical or chemical methods. This releases the DNA from each cell’s nucleus, allowing it to be purified and quantified. A total of 10uG of DNA is required for the ribose-seq protocol.[1]

Preparation of ribose-seq adaptor

The ribose-seq adaptor is made up of one dT (deoxyribonucleotide thymine) followed by an annealing sequence, two primer sequences, and a unique molecular identifier (UMI). To prepare the ribose-seq adaptor, the three oligonucleotides (oligos) are added to a solution and heated to 95-100°C. At this temperature, the oligo sequences anneal to each other, creating a suspension of ribose-seq adaptors. The amount of desalted double-stranded ribose-seq adaptor is then quantified. A concentration of 10 uM is needed for ribose-seq.[1]

Ribose-Seq

Fragmentation of Genomic DNA

Genomic DNA is fragmented using blunt-ended restriction enzyme digestion. Following digestion, fragmented DNA is purified and quantified. The restriction enzyme set and combination optimization tools for rNMP capture techniques (RESCOT) were developed to maximize rNMP capture rate.[32] The resulting solution contains a population of 500- to 3,000-bp genomic fragments with an average size of ~1.5 kb. 10ug of fragmented genomic DNA is required for subsequent dA-Tailing and Adaptor Ligation.[1]

dA-Tailing and adaptor ligation

dA molecules are bonded to the tails of each fragment, allowing the dT of the ribose-seq adaptor to bind it's corresponding dA tail in the subsequent step. Following this process, the dA-tailed DNA is purified and eluted with water. The ribose-seq adaptor, containing its randomized 8-base unique molecular identifier (UMI) and single dT overhang, is then ligated to the dA end of the target fragment by cooling the solution to 15C overnight. The sample is then purified and prepared for Alkali Treatment.[1]

Alkali treatment

During alkali treatment, double stranded DNA (dsDNA) is treated with sodium hydroxide, a strong base, to denature the dsDNA and cleave DNA at embedded ribonucleotide sites. This generates single-stranded DNA (ssDNA) fragments with exposed two ends: one 2′,3′-cyclic phosphate end and an opposite 5′-hydroxyl end. The alkaline solution is then neutralized, and the ssDNA is purified.[1]

Self-Ligation (Circularization)

Next, the ssDNA-containing solution is incubated with Arabidopsis thaliana tRNA ligase (AtRNL24). AtRNL aids in A. thalinana rRNA maturation by converting 2′,3′-cyclic phosphate ends of RNA to 2′-phosphate and ligates these to 5′-phosphate ends of RNA or DNA. Introducing AtRNL to the ssDNA solution catalyzes self-ligation of the ssDNA as the 2',3' cyclic phosphate ends ligate to the 5' phosphate end of the same ssDNA strand. The resulting cyclical ssDNA structures contain the targeted embedded ribonucleotide sequence, a UMI, and oligo sequences which make up the rest of the circular ssDNA structure.[1]

Removal of linear ssDNA

The circular ssDNA structures created by self-ligation are resistant to degradation by enzyme T5 exonuclease. However, unligated, linear ssDNA is not. Therefore, treating the products and remaining DNA fragments with T5 exonuclease selectively degrades unligated, linear ssDNA strands, while leaving self-ligated ssDNA circles intact. In this way, the unwanted linear ssDNA fragments are degraded and only successfully-ligated ssDNA circles remain.[1]

Removal of 2'Phosphate

Before PCR is possible, the 2' Phosphate groups at the ligated junctions of the ssDNA circles must be removed. In order to accomplish this, the circular ssDNA solution is treated with 2'phosphotransferase, Tpt1, removing the 2′ phosphate remaining at the ligation junction.[1]

Polymerase Chain Reaction (PCR)

Following ribose-seq, amplification is performed on each ssDNA circle using PCR amplification. In the first round of PCR, sequencing primers are introduced to the solution and the ssDNA circles are amplified. The sequencing primers selected must be compatible with the target library of choice. SeeApplications section for details. In the second round of PCR, more specific sequencing primers are introduced, and the ssDNA circles are amplified further. Polyacrylamide Gel Electrophoresis (PAGE) is performed to confirm the presence of the ribose-seq library, and to check the positive and negative controls.[1]

Size selection and gel purification

Gel purification is used to purify the ribose-seq library before sequencing. At this stage, fragments are selected by length to eliminate fragments that fall outside the target size range of interest.[1]

Use for downstream applications

After the sequences of interest are captured and purified, the samples can be sequenced and used for downstream applications. SeeApplications section for details.

Applications

Ribose-seq can be applied in numerous circumstances and for many cell types of various eukaryotic species (see: Advantages and Limitations).

Initial characterization

Ribose-seq was first utilized by Koh et al. (2015) to characterize the location of rNMPs in the mitochondrial and nuclear DNA of RNase-H2-deficient (rnh201∆) budding yeast.[1] They found that rNMPs were present at 0.449 rNMPs per kilobase (kb) in the nuclear genome and 19.5 rNMPs per kb in the mitochondrial genome. Additionally, Koh et al. (2015) were able to identify an apparent bias that exists for the incorporation of rCs and rGs into DNA. They found that in yeast cells, rCs and rGs are retained in genomic DNA at higher frequencies than rAs and rUs are, relative to the respective expected frequencies of all bases. Regardless of base type, rNMPs were often present at sites that were immediately downstream of a dA in both nuclear and mitochondrial DNA. Many additional "hotspots" for rNMP incorporation were also identified with respect to specific genomic features (e.g. at particular gene sites, strands, or chromosomal regions).

Multi-species comparisons

Expanding upon the initial characterization, Balachander et al. (2020) optimized their ribose-seq protocol and used it in conjunction with a bioinformatics tool (Ribose-Map,[33] available at https://github.com/agombolay/ribose-map) to map rNMPs throughout the mitochondrial and nuclear DNA of multiple strains of three yeast species: Saccharomyces cerevisiae, Saccharomyces paradoxus, and Schizosaccharomyces pombe.[34] In mitochondrial DNA (mtDNA) from each of these species and strains, it was demonstrated that the type of rNMP that is incorporated into a given genomic position is most impacted by the dNMP that is located at the immediate upstream position (-1). By using ribose-seq across multiple genomes, Balachander et al. were also able to elucidate various similarities and differences between yeast strains and species. The rNMP patterns in mtDNA of all strains showed low frequencies of rUs; however, S. paradoxus and S. cerevisiae indicate preferences for rC whereas S. pombe indicated preferences for rG. In nuclear DNA, the three yeast species were largely similar. Additionally, rNMPs were significantly enriched in short-nucleotide repeat sequences. This application of ribose-seq allowed for the characterization of rNMP incorporation patterns and the identification of potential drivers of this process within the sequence context.

Beyond the yeast genome

In the first ribonucleotide mapping study to be conducted in a non-yeast species, ribose-seq was used in conjunction with Ribose-Map to map rNMPs throughout the nuclear, mitochondrial, and chloroplast genomes of Chlamydomonas reinhardtii, a photosynthetic unicellular green alga.[10] In mitochondrial and chloroplast DNA of this species, rAs appeared to be preferentially incorporated over other bases, whereas nuclear DNA showed bias towards rG and rC incorporation. In all C. reinhardtii organelles, rU was consistently the least represented rNMP, which is consistent with findings from yeast mitochondria. Similar to yeast profiles, the incorporation of rNMPs in C. reinhardtii is most impacted by the dNMP immediately upstream (position -1). In this setting, the application of Ribose-seq allowed for the identification of rNMP incorporation distributions and patterns in this species for the first time.

Advantages and limitations

Advantages of ribose-seq

Because ribose-seq specifically targets embedded ribonucleotides in DNA, and does not capture RNA primers or Okazaki fragments, it is highly transferable and can effectively target embedded ribonucleotides in a wide variety of genomic DNA.[35] Therefore, ribose-seq allows for rNMP mapping in different types of genomes , including small genomic molecules or large nuclear genomes, without necessitating standardized procedures. Additionally, the specificity of ribose-seq allows for the mapping of rNMPs even in conditions in which the DNA is exposed to environmental stressors that damage the DNA by generating breaks and/or abasic sites.

Ribose-seq uses the UMI in the adaptors for high-throughput sequencing which allows to perform deduplication of the sequencing reads, and thus providing for a more accurate detection of rNMP hotspots.[34]

Ribose-seq is not specific to yeast cells; the technique could potentially be applied to any cell type from any species, making it a useful tool for elucidating the and location and pattern of rNMPs in most organisms.[36]

Ribose-seq uses common laboratory equipment and reagents that are fairly easy to source, which makes it an accessible technique for use in many molecular biology laboratories.[31]

Limitations of ribose-seq

Because it targets single rNMPs, ribose-seq cannot capture longer ribonucleotide tracts. Therefore, ribose-seq would not be an appropriate technique for mapping RNA primers or Okazaki fragments, which can be formed during DNA replication or breaks.[35]

Ribose-seq takes about 4 days to prepare the rNMP libraries from the genomic DNA.[31] The protocol is also relatively complex, requiring multiple steps.

Complementary methods

Ribose-Map

Ribose-Map is a bioinformatic tool for processing and analysis of large, complex rNMP raw sequencing datasets. [33] Ribose-Map is compatible with ribose-seq, emRiboSeq, Pu-seq, Alk-HydEn-seq, and RHII-HydEn-seq.[33] (see: Other Techniques for rNMP capture).

RESCOT (Restriction Enzyme Set and Combination Optimization Tools)

RESCOT is a computational method that calculates the genomic coverage of RNMP-captured regions for a given choice of restriction enzymes and then optimize the RE sets and select the RE set with highest coverage to maximize rNMP capture rate.[32]

Other techniques for rNMP capture

  • emRiboSeq: Uses human RNase H2 and T4 Quick ligase to capture the deoxyribonucleotide upstream of the embedded rNMP from the 5′ side.[28]
  • RHII-HydEn-seq: Uses Escherichia coli RNase HII and T4 RNA ligase to capture the embedded rNMP from its 5′ side. RHII-HydEn-seq can also capture the DNA sequence terminated with an rNMP, and Okazaki fragments.[29]
  • Pu-seq: Uses alkaline conditions and T4 RNA ligase to capture the deoxyribonucleotide downstream of the embedded rNMP from the 3′ side. Can capture Okazaki fragments.[30]
  • Alk-HydEn-seq: Uses alkaline conditions and T4 RNA ligase to capture the deoxyribonucleotide downstream of the embedded rNMP from the 3′ side. Can also capture Okazaki fragments.[37]
  • RiSQ-seq (Ribonucleotide Scanning Quantification sequencing): Uses absolute quantification to capture the upstream neighbor deoxyribonucleotide of the rNMP.[38]

References

  1. 1.00 1.01 1.02 1.03 1.04 1.05 1.06 1.07 1.08 1.09 1.10 1.11 1.12 1.13 1.14 1.15 1.16 1.17 Koh KD, Balachander S, Hesselberth JR, Storici F (2015). "Ribose-seq: global mapping of ribonucleotides embedded in genomic DNA.". Nature Methods 12 (3): 251-7, 3 p following 257. doi:10.1038/nmeth.3259. PMID 25622106. 
  2. Reijns MA, Rabe B, Rigby RE, Mill P, Astell KR, Lettice LA (2012). "Enzymatic removal of ribonucleotides from DNA is essential for mammalian genome integrity and development.". Cell 149 (5): 1008–22. doi:10.1016/j.cell.2012.04.011. PMID 22579044. 
  3. 3.0 3.1 3.2 3.3 Williams JS, Lujan SA, Kunkel TA (2016). "Processing ribonucleotides incorporated during eukaryotic DNA replication.". Nature Reviews Molecular Cell Biology 17 (6): 350–63. doi:10.1038/nrm.2016.37. PMID 27093943. 
  4. 4.0 4.1 Williams JS, Smith DJ, Marjavaara L, Lujan SA, Chabes A, Kunkel TA (2013). "Topoisomerase 1-mediated removal of ribonucleotides from nascent leading-strand DNA.". Molecular Cell 49 (5): 1010–5. doi:10.1016/j.molcel.2012.12.021. PMID 23375499. 
  5. Kim N, Huang SN, Williams JS, Li YC, Clark AB, Cho JE (2011). "Mutagenic processing of ribonucleotides in DNA by yeast topoisomerase I.". Science 332 (6037): 1561–4. doi:10.1126/science.1205016. PMID 21700875. Bibcode2011Sci...332.1561K. 
  6. "Ribose-seq identifies, locates ribonucleotides in genomic DNA". 26 January 2015. https://www.sciencedaily.com/releases/2015/01/150126112053.htm. "The findings resulted from collaboration between researchers in Storici's laboratory at the Georgia Institute of Technology -- with graduate students Kyung Duk Koh and Sathya Balachander -- and at the University of Colorado Anschutz Medical School with assistant professor Jay Hesselberth." 
  7. Kallen, Amanda (2020). "Chapter 1 - Basic genetics: mitosis, meiosis, chromosomes, DNA, RNA, and beyond". Human Reproductive Genetics. Academic Press. p. 3. ISBN 978-0-12-816561-4. https://www.sciencedirect.com/book/9780128165614/human-reproductive-genetics. 
  8. Grossman LI, Watson R, Vinograd J (1973). "The presence of ribonucleotides in mature closed-circular mitochondrial DNA.". Proceedings of the National Academy of Sciences USA 70 (12): 3339–43. doi:10.1073/pnas.70.12.3339. PMID 4202844. Bibcode1973PNAS...70.3339G. 
  9. Vengrova S, Dalgaard JZ (2006). "The wild-type Schizosaccharomyces pombe mat1 imprint consists of two ribonucleotides.". EMBO Reports 7 (1): 59–65. doi:10.1038/sj.embor.7400576. PMID 16299470. 
  10. 10.0 10.1 El-Sayed WMM, Gombolay AL, Xu P, Yang T, Jeon Y, Balachander S (2021). "Disproportionate presence of adenosine in mitochondrial and chloroplast DNA of Chlamydomonas reinhardtii.". iScience 24 (1): 102005. doi:10.1016/j.isci.2020.102005. PMID 33490913. Bibcode2021iSci...24j2005E. 
  11. Zhou ZX, Williams JS, Lujan SA, Kunkel TA (2021). "Ribonucleotide incorporation into DNA during DNA replication and its consequences.". Critical Reviews in Biochemistry and Molecular Biology 56 (1): 109–124. doi:10.1080/10409238.2020.1869175. PMID 33461360. 
  12. Champoux JJ, Gilboa E, Baltimore D (1984). "Mechanism of RNA primer removal by the RNase H activity of avian myeloblastosis virus reverse transcriptase.". Journal of Virology 49 (3): 686–91. doi:10.1128/JVI.49.3.686-691.1984. PMID 6199510. 
  13. Garg P, Stith CM, Sabouri N, Johansson E, Burgers PM (2004). "Idling by DNA polymerase delta maintains a ligatable nick during lagging-strand DNA replication.". Genes & Development 18 (22): 2764–73. doi:10.1101/gad.1252304. PMID 15520275. 
  14. 14.0 14.1 Clausen AR, Zhang S, Burgers PM, Lee MY, Kunkel TA (2013). "Ribonucleotide incorporation, proofreading and bypass by human DNA polymerase δ.". DNA Repair (Amst) 12 (2): 121–7. doi:10.1016/j.dnarep.2012.11.006. PMID 23245697. 
  15. Kasiviswanathan R, Copeland WC (2011). "Ribonucleotide discrimination and reverse transcription by the human mitochondrial DNA polymerase.". Journal of Biological Chemistry 286 (36): 31490–500. doi:10.1074/jbc.M111.252460. PMID 21778232. 
  16. McDonald JP, Vaisman A, Kuban W, Goodman MF, Woodgate R (2012). "Mechanisms employed by Escherichia coli to prevent ribonucleotide incorporation into genomic DNA by Pol V.". PLOS Genet 8 (11): e1003030. doi:10.1371/journal.pgen.1003030. PMID 23144626. 
  17. Sparks JL, Chon H, Cerritelli SM, Kunkel TA, Johansson E, Crouch RJ (2012). "RNase H2-initiated ribonucleotide excision repair.". Mol Cell 47 (6): 980–6. doi:10.1016/j.molcel.2012.06.035. PMID 22864116. 
  18. Randerath K, Reddy R, Danna TF, Watson WP, Crane AE, Randerath E (1992). "Formation of ribonucleotides in DNA modified by oxidative damage in vitro and in vivo. Characterization by 32P-postlabeling.". Mutation Research 275 (3–6): 355–66. doi:10.1016/0921-8734(92)90038-q. PMID 1383776. https://pubmed.ncbi.nlm.nih.gov/1383776. 
  19. Williams JS, Clausen AR, Nick McElhinny SA, Watts BE, Johansson E, Kunkel TA (2012). "Proofreading of ribonucleotides inserted into DNA by yeast DNA polymerase ɛ.". DNA Repair (Amst) 11 (8): 649–56. doi:10.1016/j.dnarep.2012.05.004. PMID 22682724. 
  20. Sekiguchi J, Shuman S (1997). "Site-specific ribonuclease activity of eukaryotic DNA topoisomerase I.". Molecular Cell 1 (1): 89–97. doi:10.1016/s1097-2765(00)80010-6. PMID 9659906. 
  21. Jaishree TN, van der Marel GA, van Boom JH, Wang AH (1993). "Structural influence of RNA incorporation in DNA: quantitative nuclear magnetic resonance refinement of d(CG)r(CG)d(CG) and d(CG)r(C)d(TAGCG).". Biochemistry 32 (18): 4903–11. doi:10.1021/bi00069a027. PMID 7683912. https://pubmed.ncbi.nlm.nih.gov/7683912. 
  22. Chiu HC, Koh KD, Evich M, Lesiak AL, Germann MW, Bongiorno A (2014). "RNA intrusions change DNA elastic properties and structure.". Nanoscale 6 (17): 10009–17. doi:10.1039/c4nr01794c. PMID 24992674. Bibcode2014Nanos...610009C. https://pubmed.ncbi.nlm.nih.gov/24992674. 
  23. Li Y, Breaker RR (1999). "Kinetics of RNA Degradation by Specific Base Catalysis of Transesterification Involving the 2'-Hydroxyl Group". Journal of the American Chemical Society 121 (123): 2364–5372. doi:10.1021/ja990592p. https://pubs.acs.org/doi/10.1021/ja990592p. 
  24. Clausen AR, Murray MS, Passer AR, Pedersen LC, Kunkel TA (2013). "Structure-function analysis of ribonucleotide bypass by B family DNA replicases.". Proceedings of the National Academy of Sciences USA 110 (42): 16802–7. doi:10.1073/pnas.1309119110. PMID 24082122. PMC 3801065. Bibcode2013PNAS..11016802C. https://www.ncbi.nlm.nih.gov/entre. 
  25. Wallace BD, Williams RS (2014). "Ribonucleotide triggered DNA damage and RNA-DNA damage responses.". RNA Biology 11 (11): 1340–6. doi:10.4161/15476286.2014.992283. PMID 25692233. 
  26. O'Connell K, Jinks-Robertson S, Petes TD (2015). "Elevated Genome-Wide Instability in Yeast Mutants Lacking RNase H Activity.". Genetics 201 (3): 963–75. doi:10.1534/genetics.115.182725. PMID 26400613. 
  27. Mardis ER (2008). "Next-generation DNA sequencing methods.". Annual Review of Genomics and Human Genetics 9 (1): 387–402. doi:10.1146/annurev.genom.9.081307.164359. PMID 18576944. https://pubmed.ncbi.nlm.nih.gov/18576944. 
  28. 28.0 28.1 Reijns MAM, Kemp H, Ding J, de Procé SM, Jackson AP, Taylor MS (2015). "Lagging-strand replication shapes the mutational landscape of the genome.". Nature 518 (7540): 502–506. doi:10.1038/nature14183. PMID 25624100. Bibcode2015Natur.518..502R. 
  29. 29.0 29.1 Zhou ZX, Lujan SA, Burkholder AB, Garbacz MA, Kunkel TA (2019). "Roles for DNA polymerase δ in initiating and terminating leading strand DNA replication.". Nature Communications 10 (1): 3992. doi:10.1038/s41467-019-11995-z. PMID 31488849. Bibcode2019NatCo..10.3992Z. 
  30. 30.0 30.1 Keszthelyi A, Daigaku Y, Ptasińska K, Miyabe I, Carr AM (2015). "Mapping ribonucleotides in genomic DNA and exploring replication dynamics by polymerase usage sequencing (Pu-seq).". Nature Protocols 10 (11): 1786–801. doi:10.1038/nprot.2015.116. PMID 26492137. https://pubmed.ncbi.nlm.nih.gov/26492137. 
  31. 31.0 31.1 31.2 Balachander S, Yang T, Newnam G, El-Sayed WMM, Koh KD, Storici F (2019). Capture of Ribonucleotides in Yeast Genomic DNA Using Ribose-Seq.. Methods in Molecular Biology. 2049. pp. 17–37. doi:10.1007/978-1-4939-9736-7_2. ISBN 978-1-4939-9735-0. 
  32. 32.0 32.1 Xu P, Storici F (2021). "RESCOT: Restriction Enzyme Set and Combination Optimization Tools for rNMP Capture Techniques.". Theoretical Computer Science 894: 203–213. doi:10.1016/j.tcs.2021.08.006. PMID 34924677. 
  33. 33.0 33.1 33.2 Gombolay AL, Vannberg FO, Storici F (2019). "Ribose-Map: a bioinformatics toolkit to map ribonucleotides embedded in genomic DNA.". Nucleic Acids Res 47 (1): e5. doi:10.1093/nar/gky874. PMID 30272244. 
  34. 34.0 34.1 Balachander S, Gombolay AL, Yang T, Xu P, Newnam G, Keskin H (2020). "Ribonucleotide incorporation in yeast genomic DNA shows preference for cytosine and guanosine preceded by deoxyadenosine.". Nature Communications 11 (1): 2447. doi:10.1038/s41467-020-16152-5. PMID 32415081. Bibcode2020NatCo..11.2447B. 
  35. 35.0 35.1 "Ribose-seq identifies, locates ribonucleotides in genomic DNA". 26 January 2015. https://www.sciencedaily.com/releases/2015/01/150126112053.htm. "Ribose-seq is specific to directly capturing ribonucleotides embedded in DNA and does not capture RNA primers or Okazaki fragments formed during DNA replication, breaks or abasic sites in DNA," Storici noted. "For this reason, ribose-seq has application for rNMP mapping in any genomic DNA, from large nuclear genomes to small genomic molecules such as plasmids and mitochondrial DNA, with no need of standardization procedures," she said. "It also allows mapping rNMPs even in conditions in which the DNA is exposed to environmental stressors that damage the DNA by generating breaks and/or abasic sites." 
  36. Jay Hesselberth, Francesca Storici, Kyung Duk Koh (19 May 2015). "Ribose-seq: ribonucleotides in DNA to Illumina library". Protocol Exchange (ProtocolExchange). doi:10.1038/protex.2015.044. https://protocolexchange.researchsquare.com/article/nprot-3903/v1. Retrieved 28 Feb 2022. "Ribose-seq could potentially be applied to any cell type of any organism to allow profiling of ribonucleotide incorporation into genomic DNA.". 
  37. Clausen AR, Lujan SA, Burkholder AB, Orebaugh CD, Williams JS, Clausen MF (2015). "Tracking replication enzymology in vivo by genome-wide mapping of ribonucleotide incorporation.". Nature Structural & Molecular Biology 22 (3): 185–91. doi:10.1038/nsmb.2957. PMID 25622295. 
  38. Iida T, Iida N, Sese J, Kobayashi T (2021). "Evaluation of repair activity by quantification of ribonucleotides in the genome.". Genes Cells 26 (8): 555–569. doi:10.1111/gtc.12871. PMID 33993586.