A mechanism for ramified rolling circle amplification
© Beals et al; licensee BioMed Central Ltd. 2010
Received: 9 August 2010
Accepted: 7 December 2010
Published: 7 December 2010
Amplification of single-stranded DNA circles has wide utility for a variety of applications. The two-primer ramified rolling circle amplification (RAM) reaction provides exponential DNA amplification under isothermal conditions, creating a regular laddered series of double-stranded DNA products. However, the molecular mechanism of the RAM reaction remains unexplained.
A RAM reaction model predicts exponential accumulation of a double-stranded DNA product size series, and product-size ratios, that are consistent with observed RAM reaction products. The mechanism involves generation of a series of increasing size intermediate templates; those templates produce RAM products and recursively generate smaller intermediate templates. The model allows prediction of the number of rounds of circular template replication. Real-time RAM reaction data are consistent with the model. Analysis of RAM reaction products shows exponential growth limitation consistent with the model's predictions.
The model provides a rationale for the observed products of the RAM reaction, and the molecular yield among those products. Experimental results are consistent with the model.
A variety of isothermal nucleic acid amplification technologies (reviewed in ) are now available as complements or alternatives to the polymerase chain reaction (PCR). Detection of circularized single-stranded DNA (ssDNA) molecules (padlock probes, ) has been brought to bear on applications as diverse as SNP detection [3–5], miRNA detection , and molecular diagnostics [7, 8]. Isothermal ssDNA circle amplification has been performed by single-primer initiated rolling circle amplification (RCA) and by two-primer amplification methods variously called ramification amplification (RAM, ); hyperbranched RCA , cascade RCA , or exponential RCA .
In contrast to single primer RCA  and the PCR, no detailed and predictive model for the RAM reaction has yet appeared, although illustrative diagrams of two-primer ssDNA circle amplification have been published [11, 13]. Here, we present a model of the RAM reaction that accounts for the size and number of double-stranded DNA (dsDNA) product molecules.
Results and Discussion
RAM reaction products and components
Figure 2B depicts the primer extension and strand-displacement that initiates the RAM reaction. Extension of the forward primer on the circular template displaces the 5' end of that primer; further extension creates the linear ssDNA that we designate the primary transcript. Extension of the primary transcript creates complementary binding sites for reverse primers. An extended reverse primer copies, at its 3' end, the original forward primer creating a secondary transcript.
Binding and extension of a forward primer on an extended reverse primer creates a double-stranded product whose length depends on the relative positions of the forward and reverse primers on the circular template (Figure 2B, C). We assume here without loss of generality that primer placement produces a first product that is the length of the template circle circumference, referred to here as a unit length.
A RAM model derived from a reaction diagram
Primer extension and strand displacement produce a linear ssDNA primary transcript that is a concatamer of the complement of the circular template sequence. The first RAM-specific step is reverse primer binding to the primary transcript. (We assume that primers bind as soon as a complementary sequence is fully exposed.) As primary transcript elongation creates a second reverse-primer binding site, the first-bound reverse primer is extended to the end of the primary transcript (Figure 3A, step 2).
Figure 3A, step 3, shows the extended initial reverse primer displaced from the primary transcript by the extended second-bound reverse primer. The 3' end of the displaced strand is the complement of the forward primer; binding and extension of a forward primer results in a double-stranded product (Figure 3A, step 4). That product, and subsequently-formed double-stranded products that compose the RAM reaction product ladder, are considered inert in the sense that they do not serve further as template or primer (a test of this assumption is shown in Additional File 2).
When the second reverse primer extension product is displaced it is bound twice, successively (Figure 3A, step 4 and Figure 3B, step 5), by two forward primers. The extended product from the first-bound primer is displaced by the primer-extension product from the second-bound primer (Figure 3B, steps 5 and 6).
Figure 3B shows the extended second reverse primer completely displaced from the primary transcript at step 5. The strand orientation is such that (as on the primary transcript) strand displacement that exposes a 5' non-terminal primer-binding site allows a primer and its extension product to be displaced by extension of the adjacent primer. Reverse primer 2 will ultimately produce a double-stranded product with length equal to the circumference of the template circle plus the length of the first reverse primer product (Figure 3B, step 7).
A concise description of reverse primer 3's products can be made by noting that in Figure 3B step 10, an intermediate structure is formed (on an extended forward primer) that is similar in structure to the structure in Figure 3B, step 6 (reverse primer 2, marked by asterisks in Figure 3B) except that the positions of the forward and reverse primers are exchanged. Reverse primer 3 inert products are unit-lengths 1 (step 8) and 3 (step 10), plus products identical to reverse primer 2 products. The re-creation of functionally identical structures allows us to substitute a single description of a structure for each occurrence of that structure.
Initial formalization of the RAM reaction model
Generalizing, the essential steps of the RAM reaction are:
1) Rolling circle amplification successively exposes new primer binding sites on the primary template.
2) Primer extension on the primary transcript releases secondary templates that are one unit longer than their immediate predecessor secondary transcript.
3) Successive primer binding and primer extension on each secondary template recursively produces one-unit-shorter tertiary templates and inert double-stranded DNA RAM products. The process terminates when the secondary template contains a single primer-binding site.
where N is the number of dsDNA RAM product molecules of p units at the step at which reverse-primer number r has bound the primary transcript. (All formulae assume a single initial template molecule, and should scale linearly with template number.)
where ∑ N p is the sum of products of size p produced from the initial reverse primer up to reverse primer r. That sum ideally represents the number of length p molecules in the reaction products, and predicts that the number of molecules per reaction product size-class should decrease by half for each step-increase in size.
In words, we can estimate the number of reverse-primer binding sites on the primary transcript (and thus the number of rounds of circle replication) using the number of molecules observed in a given RAM product band.
The model predicts that each reaction step adds a larger product and that the number of products of any given length doubles at each step of the reaction; and therefore that product accumulation is exponential with respect to reaction step. We presume these properties to be the basis for the log2(molarity) dependence on product number shown in Figure 1C. This result suggests that if the reaction goes to completion the number of molecules in each DNA band (Figure 1A) should be half the number of molecules in the preceding band.
However, as an exponential process, the RAM reaction is not expected to go to completion; and the effects of reactant-limitation or changed reaction environment on RAM product equilibrium yield and size-distribution became a question for empirical investigation. To characterize the RAM reaction in light of this model we monitored the RAM reaction in real-time, and measured RAM-product mass-distribution at various reaction time-points.
Real-time analysis of the RAM reaction
RAM and PCR produce similar fluorescence-signal curves in real-time reactions. However, there are substantive differences between the isothermal RAM reaction and the thermocycling that is inherent to PCR. The RAM reaction is analyzed as a continuous time-dependent process rather than a discrete cycle-dependent process; the time to reach a threshold fluorescence is designated the response time (Rt)  corresponding to the cycle threshold (Ct) of real-time PCR. Further, in PCR the template is invariant, whereas the RAM reaction generates a dynamic collection of templates (Figure 3) as the reaction progresses.
To see whether template number affected RAM signal characteristics we measured the real-time kinetics of reactions initiated at various template levels. Because real-time instruments and analytical algorithms were developed for the PCR, we investigated whether mechanistic differences between PCR and RAM would affect real-time RAM signal analysis.
Model fit, response time, and analysis data for a real-time RAM reaction.
Figure 4B shows the log-linear slope of Rt vs. log2(template copy number). (Similar plots yield the efficiency of a PCR reaction .) The figure shows that response time is a log-linear function of template-dilution. Table 1 (efficiency) and diamond symbols in Figure 4B show a measure of the efficiency of the reaction as determined from the fitted models; roughly equal efficiency was obtained for all reactions. Together, these data suggest that the RAM reaction rate is uncorrelated to input template number.
These observations are consistent with the RAM reaction model. At low template number, reaction time differences may result from stochastic forward primer binding time. At very high template numbers it is conceivable that primary transcript elongation may be inhibited by reactant depletion. Between those extremes dsDNA products result from secondary and tertiary transcript processing and we see no compelling reason for distinct processing modes of these transcripts as a result of input template number. However the RAM reaction rate is a complex function of several variables (unpublished data) with no a priori expectation of maximum rate, as in the PCR, and some reaction conditions may lead to template-number-sensitive reactions. However, these data suggest no necessary dependence of real-time RAM reaction kinetics on input template number.
Detailed comparison of real-time data vs. the sigmoid model via plots of residuals (Additional File 3) show that the models' deviation from the data, although small with respect to the overall signal, has a distinct time-dependent pattern. Although Figure 4A represents RAM data collected as for PCR, isothermal reactions can be monitored essentially continuously. Model fitting in PCR is required to find the exponential amplification cycles at which the Ct is ideally measured, because real-time PCR data frequently have only a few cycles of exponential amplification . Continuous data collection from isothermal reactions allows simpler non-model based direct detection of exponential amplification .
Characterization of a RAM-reaction time-course
The RAM model predicts that if the reaction goes to completion each product peak should have half the number of molecules as its smaller predecessor. Initial capillary electrophoresis assessment of RAM products did not show the ideal log-linear -1 slope of log2(molecule number) vs. product peak number. We compared product ratio slopes from RAM reaction time-course samples to determine how reaction dynamics affected deviation from the ideal slope.
An explanation for deviation from the ideal slope but maintaining log-linearity follows from the reaction model and stoichiometric calculation. From RAM product yield and equation 3 we infer (Additional File 1) that at least 40 to 45 reverse primers have initiated secondary transcripts to yield the products shown in Figure 4A. Complete resolution of those transcripts, under realistic conditions, would require reactant amounts that would sufficiently reduce the reactant concentration to stop the reaction (calculation details are given in Additional File 1).
The proposed model predicts that, as secondary transcripts are recursively processed into smaller tertiary transcripts, larger dsDNA products are made before smaller dsDNA products. However, smaller secondary transcripts may be fully resolved into dsDNA RAM products before reactants become limiting. Our conception of the results shown in Figure 5 is that the observed RAM products are the combination of products from fully resolved smaller secondary transcripts with larger products from partially resolved larger secondary transcripts.
One motive for the time-course experiments described above was to look for evidence that would falsify the model. In the hundreds of RAM reactions we have characterized each has had the form of the data shown in Figure 5B - a log-linear relation with negative slopes always greater than -1.
This RAM model originated with an exhaustive diagram of the initial stages of the reaction, and is based on the RCA model and on the properties of a strand-displacing DNA polymerase. The utility of the Bst polymerase for the RAM reaction may lie in the balance of that enzyme's relative rate of primary transcript production and subsequent template processing. The orderly production of incrementally longer secondary transcripts may not occur if the primary transcript is produced faster than secondary transcripts can be initiated.
Our purpose is to provide a mechanistic working hypothesis for the RAM reaction. Here, we make quantitative predictions only for the theoretical reaction model; our description of the reaction products, although based on the model, is qualitative. An obvious extension of the model would be simulation of the relative contributions of completely processed vs. partially processed secondary and tertiary transcripts.
To our knowledge, this RAM reaction model is unique in its prediction of a molecular binary-power series and in invocation of a recursive molecular mechanism. These concepts are familiar from computational algorithms, and the RAM reaction may prove useful in molecular computation in addition to its current applications.
RAM reaction time-courses were done in RAM buffer 1: 1× New England Biolabs (NEB) Thermopol II buffer, 2 mM MgSO4, 0.2 mM each dNTP; 2% DMSO; 0.75 uM forward primer Cpr4PrFwd80_20; 1 uM Cpr4PrRvs09_19; 50 molecules circularized template Cpr4Stchy_ch154 per reaction; Bst DNA polymerase, large fragment (NEB), 0.13 units/ul.
Real-time RAM reactions were done in RAM buffer 2: [1× New England Biolabs (NEB) Thermopol II buffer, 2 mM MgSO4, 0.2 mM each dNTP; 5% DMSO; 40 mM NaCl; 0.6× SYBR-Green (Molecular Probes), 1 nM fluorescein (BioRad), 0.75 uM forward primer Cpr2Fwd47_20; 1 uM Cpr3Rvs50_19; circularized template Cpr3Stchy_ch381 in the range of 0 (zero) to 107 molecules per reaction; Bst DNA polymerase, large fragment (NEB), 0.13 units/ul.
RAM reaction time-course
A RAM reaction source-tube was assembled with all reaction components except primers and template. 40 microliters (ul) RAM reaction source was dispensed to each of 40 0.5 ml capacity reaction tubes. 10 ul of primer-template mix was added to each reaction tube, quickly mixed, then the reaction tube was held at 60°C. After incubation for the specified time (here, 33 to 72 minutes) 1 ul of 100 mM Na2EDTA was added and tubes were stored at 4°C.
Curve-fitting to real-time data
Fluorescence vs. time data were imported as an R data frame and analyzed using the qpcR package . Models were fit with the pcrbatch() function with arguments (...fct = l5(), retPar = TRUE). The parameters obtained for each reaction were used to generate the model curves at each time-point. cpD1, cpD2, experimental efficiency, and R2(adj) were taken from the table generated by pcrbatch().
C-probes and primers
Cpr4Stchy_ch154 (95 DNA bases (b.)) is a padlock probe  that contains, 5' to 3': 25 bases of NCBI-nucleotide record AF081468 (encoding Stachybotrys chartarum ribosomal RNA) bases 486-510; 60 bases of internal sequence; and 10 bases of AF081468, 476-485.
Cpr4Stchy_ch154: 5'AGTATTCTCT GAGTGGCAAA CGCAAATTAA ACTCCTCCTG TGATGTCTAA TATGGTACCG
TTCGTTAGAA TCAGCTGGCA AATTCTGTTT TTTTC.
Cpr3Stchy_ch381 (100 b) is a padlock probe that contains, 5' to 3': 24 bases of NCBI-nucleotide record AF081468 (encoding Stachybotrys chartarum ribosomal RNA) bases 397-420; 60 bases of internal sequence; and 16 bases of AF081468, 381-396.
Cpr3Stchy_ch381: TACCTATCGT TGCTTCGGCG GGAATGAAGC TTGTCCTAGT GTGTCAGTCG CACGCTTACC
AAGAGCAACT ACACGAACAG CTGTACCCTT ATGTGAACCG.
RAM amplification primers are designated by a descriptive text prefix followed by the Cprobe base-number that initiates synthesis of the 5' complement of a forward primer on the Cprobe template, or by the first Cprobe base-number of the reverse primer subsequence; followed by an underscore separator, followed by the oligo-length.
Cpr4Fwd80_20: 5' TGCCAGCTGATTCTAACGAA
Cpr4PrRvs09_19: 5' CTGAGTGGCAAACGCAAAT
Cpr2Fwd47_23: 5' CTGACACACTAGGACAAGCTTCA
Cpr3Rvs50_20: 5' GCACGCTTACCAAGAGCAAC
Agarose gel electrophoresis
A 300 ml, 1.2% agarose in 0.5 × TBE was cast and run in an Owl A2 electrophoresis chamber. 4 ul RAM reaction (58 minute sample from the time course described above) and 10 ul BRL 1 kb ladder were loaded to the wells of the gel lanes shown in Figure 1. The gel was run in 0.5× TBE running buffer for 200 minutes at 180 Vdc., stained in 1× SYBR-Gold (Molecular Probes) in 0.5 × TBE for 20 minutes, then imaged on a Spectroline model TI-312E UV-transilluminator through a Wratten No. 15 filter (sold by Molecular Probes as a SYBR Green/Gold photographic filter).
The gel image was captured with a Nikon 990 (E990V1.0) digital camera in aperture-priority mode; F6.0 at effective focal length f13.7 mm, 5.10 second exposure, "FULL FINE" mode. All digital manipulations were applied uniformly until separation of selected lanes for image assembly. The image was rotated to vertical orientation, reversed to create a black-on-white image; then contrast and brightness were adjusted. The chosen lanes were cropped and placed in the final figure.
Capillary electrophoresis (CE)
DNA mass and size measurements were done in an Agilent Bioanalyzer 2100 and analyzed with software version B02.06 SI 418 Patch level 01. DNA1000 (Series II) assays were performed according to the manufacturer's protocol. 1 ul RAM reaction mix was loaded to each CE chamber.
CE data processing: Electropherogram representations of the data were reviewed to correct occasional incorrect assignment of electrophoretic peak to internal markers. The data was exported as text files for further analysis. RAM products were assigned peak numbers, designating the lowest molecular weight product as peak 1. Slopes were determined by linear regression of log2(nMolar DNA) vs. peak number for the first 8 RAM product peaks.
Analysis of CE data: Figure 4B data panels were constructed in JMP (SAS Institute, Version 4.0.4) as a oneway analysis of variance of log2(nMolar concentration) vs. peak number for each time point shown. Figure 4C shows slope (determined as described above) for each reaction containing sufficient DNA for detection in the Agilent instrument.
We thank Dr. Harry McCoy and two anonymous reviewers for critical reading of the manuscript.
- Gill P, Ghaemi A: Nucleic acid isothermal amplification technologies: a review. Nucleosides Nucleotides Nucleic Acids. 2008, 27: 224-43. 10.1080/15257770701845204View ArticlePubMedGoogle Scholar
- Nilsson M, Malmgren H, Samiotaki M, Kwiatkowski M, Chowdhary BP, Landegren U: Padlock probes: circularizing oligonucleotides for localized DNA detection. Science. 1994, 265: 2085-8. 10.1126/science.7522346View ArticlePubMedGoogle Scholar
- Hardenbol P, Baner J, Jain M, Nilsson M, Namsaraev EA, Karlin-Neumann GA, Fakhrai-Rad H, Ronaghi M, Willis TD, Landegren U: Multiplexed genotyping with sequence-tagged molecular inversion probes. Nat Biotechnol. 2003, 21: 673-8. 10.1038/nbt821View ArticlePubMedGoogle Scholar
- Pickering J, Bamford A, Godbole V, Briggs J, Scozzafava G, Roe P, Wheeler C, Ghouze F, Cuss S: Integration of DNA ligation and rolling circle amplification for the homogeneous, end-point detection of single nucleotide polymorphisms. Nucleic Acids Res. 2002, 30: e60- 10.1093/nar/gnf060View ArticlePubMedPubMed CentralGoogle Scholar
- Alsmadi OA, Bornarth CJ, Song W, Wisniewski M, Du J, Brockman JP, Faruqi AF, Hosono S, Sun Z, Du Y: High accuracy genotyping directly from genomic DNA using a rolling circle amplification based assay. BMC Genomics. 2003, 4: 21- 10.1186/1471-2164-4-21View ArticlePubMedPubMed CentralGoogle Scholar
- Yao B, Li J, Huang H, Sun C, Wang Z, Fan Y, Chang Q, Li S, Xi J: Quantitative analysis of zeptomole microRNAs based on isothermal ramification amplification. RNA. 2009, 15: 1787-94. 10.1261/rna.1555209View ArticlePubMedPubMed CentralGoogle Scholar
- Zhang W, Cohenford M, Lentrichia B, Isenberg HD, Simson E, Li H, Yi J, Zhang DY: Detection of Chlamydia trachomatis by isothermal ramification amplification method: a feasibility study. J Clin Microbiol. 2002, 40: 128-32. 10.1128/JCM.40.1.128-132.2002View ArticlePubMedPubMed CentralGoogle Scholar
- Li F, Zhao C, Zhang W, Cui S, Meng J, Wu J, Zhang DY: Use of ramification amplification assay for detection of Escherichia coli O157:H7 and other E. coli Shiga toxin-producing strains. J Clin Microbiol. 2005, 43: 6086-90. 10.1128/JCM.43.12.6086-6090.2005View ArticlePubMedPubMed CentralGoogle Scholar
- Zhang DY, Brandwein M, Hsuih TC, Li H: Amplification of target-specific, ligation-dependent circular probe. Gene. 1998, 211: 277-85. 10.1016/S0378-1119(98)00113-9View ArticlePubMedGoogle Scholar
- Lizardi PM, Huang X, Zhu Z, Bray-Ward P, Thomas DC, Ward DC: Mutation detection and single-molecule counting using isothermal rolling-circle amplification. Nat Genet. 1998, 19: 225-32. 10.1038/898View ArticlePubMedGoogle Scholar
- Thomas DC, Nardone GA, Randall SK: Amplification of padlock probes for DNA diagnostics by cascade rolling circle amplification or the polymerase chain reaction. Arch Pathol Lab Med. 1999, 123: 1170-6.PubMedGoogle Scholar
- Fire A, Xu SQ: Rolling replication of short DNA circles. Proc Natl Acad Sci USA. 1995, 92: 4641-5. 10.1073/pnas.92.10.4641View ArticlePubMedPubMed CentralGoogle Scholar
- Isaksson A, Landegren U: Accessing genomic information: alternatives to PCR. Curr Opin Biotechnol. 1999, 10: 11-5. 10.1016/S0958-1669(99)80003-2View ArticlePubMedGoogle Scholar
- Burg JL, Cahill PB, Kutter M, Stefano JE, Mahan DE: Real-time fluorescence detection of RNA amplified by Q beta replicase. Anal Biochem. 1995, 230: 263-72. 10.1006/abio.1995.1473View ArticlePubMedGoogle Scholar
- Spiess AN, Feig C, Ritz C: Highly accurate sigmoidal fitting of real-time PCR data by introducing a parameter for asymmetry. BMC Bioinformatics. 2008, 9: 221- 10.1186/1471-2105-9-221View ArticlePubMedPubMed CentralGoogle Scholar
- Zhao S, Fernald RD: Comprehensive algorithm for quantitative real-time polymerase chain reaction. J Comput Biol. 2005, 12: 1047-64. 10.1089/cmb.2005.12.1047View ArticlePubMedPubMed CentralGoogle Scholar
- Ritz C, Spiess AN: qpcR: an R package for sigmoidal model selection in quantitative real-time polymerase chain reaction analysis. Bioinformatics. 2008, 24: 1549-51. 10.1093/bioinformatics/btn227View ArticlePubMedGoogle Scholar
- R Development Core Team: R: A language and environment for statistical computing. ISBN 3-900051-07-0, R Foundation for Statistical Computing, Vienna, Austria, http://www.R-project.org
- Pfaffl MW: A new mathematical model for relative quantification in real-time RT-PCR. Nucleic Acids Res. 2001, 29: e45- 10.1093/nar/29.9.e45View ArticlePubMedPubMed CentralGoogle Scholar
- Beals TP: Application to RAM Amplification of Real-time Analysis Methodologies Developed for PCR. qPCR 2009 Event. Diagnostics & Molecular Markers. 2009, P021-Freising, GermanyGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.