Real-time PCR quantification using a variable reaction efficiency model

The quantitative real-time polymerase chain reaction (PCR) remains a cornerstone technique in gene expression analysis and sequence characterization. Despite the importance of the approach to experimental biology, the confident assignment of reaction efficiency to the early cycles of real-time PCR reactions remains problematic. Considerable noise may be generated when few cycles in the amplification are available to estimate peak efficiency. An alternate approach that uses data from beyond the log–linear amplification phase is explored in this article with the aim of reducing noise and adding confidence to efficiency estimates. PCR reaction efficiency is regressed to estimate the per-cycle profile of an asymptotically departed peak efficiency even when this is not closely approximated in the measurable cycles. The process can be repeated over replicates to develop a robust estimate of peak reaction efficiency. This leads to an estimate of the maximum reaction efficiency that may be considered primer design specific. Using a series of biological scenarios, we demonstrate that this approach can provide an accurate estimate of initial template concentration.     

长片段融合PCR在构建烟曲霉rho1基因回补株中的应用

目的融合PCR是一种常用的构建重组片段或重组质粒的手段,但长片段融合PCR的难度较大。文中将探讨长片段融合PCR过程中引物设计及扩增条件对产物的影响。方法以构建烟曲霉rho 1基因回补株为例,采用融合PCR的方法扩增重组片段(长达6.5 kb),在引物设计时引入不同大小的同源区,并设置不同的扩增体系。结果当设计引物的同源区为35 bp,选用具有高扩增效率、高保真性的DNA聚合酶,以及各片段在融合PCR反应体系中的浓度为15 ng/μL时,实现了长达6.5 kb的片段扩增并完成了烟曲霉rho 1回补株的构建。结论在合适的PCR引物设计、片段浓度配比及聚合酶条件下,长片段融合PCR在丝状真菌的基因敲除及回补株的构建中是一种非常有效的工具。     

美国黑核桃SSR反应体系优化

优化SSR-PCR反应体系是黑核桃（Juglans nigra L.）SSR基因鉴定和群体遗传等研究的基础。本研究通过对PCR反应中Mg²⁺浓度、牛血清白蛋白（Bovine Serum Albumin,BSA）浓度、Taq聚合酶用量、dNTPs浓度、引物浓度和模板DNA量的组合以及PCR程序组合试验,确定了黑核桃SSR的最佳反应体系,即在10 μL的PCR反应体系中,含10 ng模板DNA,0.1 mg·mL^-1牛血清蛋白（BSA）,0.25 mmol·L^-1 dNTPs,1.5 mmol·L^-1 Mg²⁺ 1 μL 10X Taq DNA聚合酶反应缓冲液,0.5 U Taq聚合酶,1.0 mmol·L^-1单对引物（0.5 mmol·L^-13对引物）。SSR-PCR反应扩增程序为：94℃变性3 min;93℃变性15 s,50℃或者53.5℃退火1 min,72℃延伸30 s,32个循环;72℃后延伸10 min,置4℃保存。利用此反应体系对黑核桃进行PCR扩增并电泳检测,其结果清晰、稳定、可靠,适合进一步对黑核桃群体遗传、基因型鉴定和分子生态研究。     

High GC Contents of Primer 5′-End Increases Reaction Efficiency in Polymerase Chain Reaction

Many studies have suggested that regulation of the polymerase chain reaction (PCR) is influenced by several factors. However, the understanding of reaction efficiency factors is not sufficient. Here we propose that high GC contents of primer 5′-end increases reaction efficiency in PCR. Using 71 primers (45 pairs), we analyzed factors that affect reaction efficiency, and statistically tested the correlation between the amplification signals and several factors. As a result, there were significant correlations between the amplification signals and the GC contents in the first 1～3 bps of primer 5′-end.     

葡萄ISSR-PCR反应体系的建立与优化

采用正交设计L9(34)对影响葡萄ISSR-PCR反应体系的4个因素(dNTP、TaqDNA聚合酶、引物、模板DNA)在3个浓度水平上进行试验,并通过直观分析初步确定其反应体系;在此基础上,通过单因素试验探讨了dNTP、TaqDNA聚合酶、引物、模板DNA、退火温度及循环次数等因素或条件对葡萄ISSR-PCR扩增结果的影响,确定最佳反应水平。最终建立了葡萄ISSR-PCR扩增的最佳反应体系:在25μL的反应体系中,dNTP浓度0.2 mmol/L,TaqDNA聚合酶的用量0.5 U,引物浓度0.4mmol/L,DNA模板用量40 ng。反应程序:94℃预变性5 min;94℃变性1 min,52℃退火1 min,72℃延伸1 min 30 s,40次循环;最后72℃延伸10 min,10℃保存。     

A general model of error-prone PCR

In this paper, we generalize a previously-described model of the error-prone polymerase chain reaction (PCR) reaction to conditions of arbitrarily variable amplification efficiency and initial population size. Generalisation of the model to these conditions improves the correspondence to observed and expected behaviours of PCR, and restricts the extent to which the model may explore sequence space for a prescribed set of parameters. Error-prone PCR in realistic reaction conditions is predicted to be less effective at generating grossly divergent sequences than the original model. The estimate of mutation rate per cycle by sampling sequences from an in vitro PCR experiment is correspondingly affected by the choice of model and parameters.     

Sequence-Dependent Biophysical Modeling of DNA Amplification

A theoretical framework for prediction of the dynamic evolution of chemical species in DNA amplification reactions, for any specified sequence and operating conditions, is reported. Using the polymerase chain reaction (PCR) as an example, we developed a sequence- and temperature-dependent kinetic model for DNA amplification using first-principles biophysical modeling of DNA hybridization and polymerization. We compare this kinetic model with prior PCR models and discuss the features of our model that are essential for quantitative prediction of DNA amplification efficiency for arbitrary sequences and operating conditions. Using this model, the kinetics of PCR is analyzed. The ability of the model to distinguish between the dynamic evolution of distinct DNA sequences in DNA amplification reactions is demonstrated. The kinetic model is solved for a typical PCR temperature protocol to motivate the need for optimization of the dynamic operating conditions of DNA amplification reactions. It is shown that amplification efficiency is affected by dynamic processes that are not accurately represented in the simplified models of DNA amplification that form the basis of conventional temperature cycling protocols. Based on this analysis, a modified temperature protocol that improves PCR efficiency is suggested. Use of this sequence-dependent kinetic model in a control theoretic framework to determine the optimal dynamic operating conditions of DNA amplification reactions, for any specified amplification objective, is discussed.     

Sequence-Dependent Biophysical Modeling of DNA Amplification

A theoretical framework for prediction of the dynamic evolution of chemical species in DNA amplification reactions, for any specified sequence and operating conditions, is reported. Using the polymerase chain reaction (PCR) as an example, we developed a sequence- and temperature-dependent kinetic model for DNA amplification using first-principles biophysical modeling of DNA hybridization and polymerization. We compare this kinetic model with prior PCR models and discuss the features of our model that are essential for quantitative prediction of DNA amplification efficiency for arbitrary sequences and operating conditions. Using this model, the kinetics of PCR is analyzed. The ability of the model to distinguish between the dynamic evolution of distinct DNA sequences in DNA amplification reactions is demonstrated. The kinetic model is solved for a typical PCR temperature protocol to motivate the need for optimization of the dynamic operating conditions of DNA amplification reactions. It is shown that amplification efficiency is affected by dynamic processes that are not accurately represented in the simplified models of DNA amplification that form the basis of conventional temperature cycling protocols. Based on this analysis, a modified temperature protocol that improves PCR efficiency is suggested. Use of this sequence-dependent kinetic model in a control theoretic framework to determine the optimal dynamic operating conditions of DNA amplification reactions, for any specified amplification objective, is discussed.     

多倍体罗汉果PCR-RFLP参数的优化与应用

以多倍体罗汉果DNA为材料,采用L16(4~5)正交组合试验和单因素梯度试验,研究Mg~(2+)、dNTP、引物、Taq DNA聚合酶、模板DNA浓度和退火温度、循环次数等对PCR扩增结果以及内切酶量、酶切时间对酶切反应的影响。结果表明,多倍体罗汉果RFLP最优PCR反应体系和扩增参数为:在25μL扩增反应体系中,10×Buffer 2.5μL,MgCl_2 1.5 mmol/L,dNTP 0.2 mmol/L,引物0.1μmol/L,Taq DNA聚合酶2.0 U,模板DNA 60 ng;退火温度为56℃,循环次数为35次。酶切反应体系:内切酶10×Buffer 2.0μL,内切酶5.0U,PCR产物15μL,超纯水补至20μL;酶切时间2 h。     

The choice of universal primers and the characteristics of the species mixture determine when DNA metabarcoding can be quantitative

DNA metabarcoding is a technique used to survey biodiversity in many ecological settings, but there are doubts about whether it can provide quantitative results, that is, the proportions of each species in the mixture as opposed to a species list. While there are several experimental studies that report quantitative metabarcoding results, there are a similar number that fail to do so. Here, we provide the rationale to understand under what circumstances the technique can be quantitative. In essence, we simulate a mixture of DNA of S species with a defined initial abundance distribution. In the simulated PCR, each species increases its concentration following a certain amplification efficiency. The final DNA concentration will reflect the initial one when the efficiency is similar for all species; otherwise, the initial and final DNA concentrations would be poorly related. Although there are many known factors that modulate amplification efficiency, we focused on the number of primer–template mismatches, arguably the most important one. We used 15 common primers pairs targeting the mitochondrial COI region and the mitogenomes of ca. 1,200 insect species. The results showed that some primers pairs produced quantitative results under most circumstances, whereas some other primers failed to do so. In conclusion, depending on the primer pair used in the PCR amplification and on the characteristics of the mixture analysed (i.e., high species richness, low evenness), DNA metabarcoding can provide a quantitative estimate of the relative abundances of different species.     

Increased efficiency of genetic profiling through quantity and quality assessment of fluorescently labeled oligonucleotide primers

Optimizing the amount of primer to use in PCR amplification is one of the most important steps when developing protocols for genetic profiling, where subtle changes in primer concentration result in major impacts on the amount of desired product that is amplified. However; there are frequently discrepancies between the reported and actual quantity of primers delivered by suppliers, resulting in a need for re-optimization of conditions between primer orders and limiting the ability to standardize conditions between laboratories. To increase the consistency of genetic profiling protocols, we have developed a simple method to assess the quantity and quality of fluorescently labeled primers and therefore standardize reaction conditions through time and across laboratories. The method is based on analysis by electrophoresis with an automated fluorescent DNA analyzer.     

Stepwise Kinetic Equilibrium Models of Quantitative Polymerase Chain Reaction

ABSTRACT: BACKGROUND: Numerous models for use in interpreting quantitative PCR (qPCR) data are present in recent literature. The most commonly used models assume the amplification in qPCR is exponential and fit an exponential model with a constant rate of increase to a select part of the curve. Kinetic theory may be used to model the annealing phase and does not assume constant efficiency of amplification. Mechanistic models describing the annealing phase with kinetic theory offer the most potential for accurate interpretation of qPCR data. Even so, they have not been thoroughly investigated and are rarely used for interpretation of qPCR data. New results for kinetic modeling of qPCR are presented. RESULTS: Two models are presented in which the efficiency of amplification is based on equilibrium solutions for the annealing phase of the qPCR process. Model 1 assumes annealing of complementary targets strands and annealing of target and primers are both reversible reactions and reach a dynamic equilibrium. Model 2 assumes all annealing reactions are nonreversible and equilibrium is static. Both models include the effect of primer concentration during the annealing phase. Analytic formulae are given for the equilibrium values of all single and double stranded molecules at the end of the annealing step. The equilibrium values are then used in a stepwise method to describe the whole qPCR process. Rate constants of kinetic models are the same for solutions that are identical except for possibly having different initial target concentrations. Analysis of qPCR curves from such solutions are thus analyzed by simultaneous non-linear curve fitting with the same rate constant values applying to all curves and each curve having a unique value for initial target concentration. The models were fit to two data sets for which the true initial target concentrations are known. Both models give better fit to observed qPCR data than other kinetic models present in the literature. They also give better estimates of initial target concentration. Model 1 was found to be slightly more robust than model 2 giving better estimates of initial target concentration when estimation of parameters was done for qPCR curves with very different initial target concentration. Both models may be used to estimate the initial absolute concentration of target sequence when a standard curve is not available. CONCLUSIONS: It is argued that the kinetic approach to modeling and interpreting quantitative PCR data has the potential to give more precise estimates of the true initial target concentrations than other methods currently used for analysis of qPCR data. The two models presented here give a unified model of the qPCR process in that they explain the shape of the qPCR curve for a wide variety of initial target concentrations.     

Universal and blocking primer mismatches limit the use of high‐throughput DNA sequencing for the quantitative metabarcoding of arthropods

The quantification of the biological diversity in environmental samples using high‐throughput DNA sequencing is hindered by the PCR bias caused by variable primer–template mismatches of the individual species. In some dietary studies, there is the added problem that samples are enriched with predator DNA, so often a predator‐specific blocking oligonucleotide is used to alleviate the problem. However, specific blocking oligonucleotides could coblock nontarget species to some degree. Here, we accurately estimate the extent of the PCR biases induced by universal and blocking primers on a mock community prepared with DNA of twelve species of terrestrial arthropods. We also compare universal and blocking primer biases with those induced by variable annealing temperature and number of PCR cycles. The results show that reads of all species were recovered after PCR enrichment at our control conditions (no blocking oligonucleotide, 45 °C annealing temperature and 40 cycles) and high‐throughput sequencing. They also show that the four factors considered biased the final proportions of the species to some degree. Among these factors, the number of primer–template mismatches of each species had a disproportionate effect (up to five orders of magnitude) on the amplification efficiency. In particular, the number of primer–template mismatches explained most of the variation (~3/4) in the amplification efficiency of the species. The effect of blocking oligonucleotide concentration on nontarget species relative abundance was also significant, but less important (below one order of magnitude). Considering the results reported here, the quantitative potential of the technique is limited, and only qualitative results (the species list) are reliable, at least when targeting the barcoding COI region.     

Optimization and design of oligonucleotide setup for strand displacement amplification

Several advantages of strand displacement amplification (SDA) as an all-purpose DNA amplification reaction are due to it isothermal mechanism. The major problem of isothermal amplification mechanism is the accumulation of non-predictable byproduct especially for longer incubation time and low concentrations of initial template DNA. New theoretical strategies to tackle the difficulties regarding the specificity of the reaction are experimentally verified. Besides improving the reaction conditions, the stringency of primer hybridization can be distinctly improved by computer based sequence prediction algorithms based on the thermodynamic stability of DNA hybrid a described by the partition function of the hybridization reaction. An alternative SDA mechanism, with sequences developed by this means is also investigated.     

DNA analysis with multiplex microarray-enhanced PCR

We have developed a highly sensitive method for DNA analysis on 3D gel element microarrays, a technique we call multiplex microarray-enhanced PCR (MME-PCR). Two amplification strategies are carried out simultaneously in the reaction chamber: on or within gel elements, and in bulk solution over the gel element array. MME-PCR is initiated by multiple complex primers containing gene-specific, forward and reverse, sequences appended to the 3′ end of a universal amplification primer. The complex primer pair is covalently tethered through its 5′ end to the polyacryl- amide backbone. In the bulk solution above the gel element array, a single pair of unattached universal primers simultaneously directs pseudo-monoplex PCR of all targets according to normal solution-phase PCR. The presence of a single universal PCR primer pair in solution accelerates amplification within gel elements and eliminates the problem of primer interference that is common to conventional multiplex PCR. We show 10⁶-fold amplification of targeted DNA after 50 cycles with average amplification efficiency 1.34 per cycle, and demonstrate specific on-chip amplification of six genes in Bacillus subtilis. All six genes were detected at 4.5 pg of bacterial genomic DNA (equivalent to 10³ genomes) in 60 independent amplification reactions performed simultaneously in single reaction chamber.     

重要肠道病原菌多重PCR-基因芯片检测方法研究

目的:建立并初步评价一种针对重要肠道病原菌的多重PCR 基因芯片检测方法。方法：对筛选出的特异引物进行多重PCR优化,将引物分别按种属内混合和种属间混合的方案排查引物间的竞争性抑制现象,再将不同菌属的模板混合,用相对应的混合引物扩增,探寻高效特异的引物组合。分别掺入和不掺入荧光素,验证其对混合PCR反应的影响,并与芯片杂交,探寻多重PCR扩增效率对芯片杂交的影响。分析不同数量引物组合产生的杂交结果,筛选出无交叉反应的最优引物组合。结果:种属内引物混合均得到特异性扩增结果。种属间混合霍乱弧菌和空肠弯曲菌得到部分预期条带,随着混合引物数量的增加,交叉抑制现象也增多。杂交信号强度随多重PCR扩增效率的增加而增强。反应中掺入荧光素的扩增条带产量低于无荧光素的产物。可将35对混合引物拆成3个体系分别标记样品,以避免假阴性结果。结论:PCR反应中掺入荧光素降低扩增效率和杂交效率,但并不影响对杂交结果的判读和数据分析。基因芯片杂交信号强度取决于多重PCR的扩增效率。肠道病原菌多重PCR 基因芯片检测方法具有较高的特异性,混合PCR可以分别按照种属内和种属间的引物组合方案用于多病原的筛检。该基因芯片检测可以采用3个引物体系完成样品标记。     

A quantitative assay for assessing allelic proportions by iterative gap ligation.

A variety of techniques are currently available for detecting point mutations in DNA. These techniques are frequently not sensitive enough to be applied as quantitative assays in evaluation of relative occurrence of alleles in cases of polymorphism or when variations in allelic gene expression are being evaluated at the level of RNA. We report here the establishment of an iterative gap ligation (IGL) assay that is both quantitative and sensitive. The design of the assay is such that ligation of an upstream to a downstream primer across a single nucleotide gap will only occur if the gap is filled with a deoxynucleotide complementary to the wild-type or mutant sequence. Under conditions in which excess upstream primer saturates the template concurrently with limiting amounts of downstream primer quantitative ligation is absolutely dependent on provision of the appropriate gap filling nucleotide. When gap ligation occurs in a single incubation, or cycle, the amount of ligated product is a linear function of the relative amount of mutant sequence, with a sensitivity and detection limit of approximately 3% over a range of relative concentrations of 0-100%. When the reaction occurs over multiple cycles, or iterations, gap ligation becomes a non-linear function such that small changes in the relative proportions of alleles produce a disproportionately large amount of ligation. As a consequence, the sensitivity and limits of detection of the assay improve to 0.2% after only 8 cycles. The development of this assay provides a unique means of quantifying allelic polymorphisms in both DNA and RNA (after initial amplification by PCR or RT-PCR) and should be applicable to any experimental settings in which nucleic acids from tissues or mixed populations of cells are being evaluated.     

Simultaneous measurement of multiple mRNAs with a single control by quantitative competitive reverse transcriptase-polymerase chain reaction: glucose transporters Glut1 and Glut4

Theoretical considerations for extending the application of quantitative competitive polymerase chain reaction (qc-PCR) to include the simultaneous measurement of multiple mRNAs, specifically the mammalian glucose transporters Glut1 and Glut4, are presented with experimental data in which the accuracy and flexibility of the system are examined. This method reliably measures changes in the initial concentration for each of three target DNA sequences. The reaction is not acutely sensitive to variations in either the primer sites or internal sequence, and although the initial concentrations of the target DNAs did affect the relative amplification efficiencies, the effect was limited and did not prohibit quantification. This PCR system was able to reliably detect differences as little as 50% in the initial concentration of the Glut1 target DNA sequence. Therefore, with the appropriate controls, PCR can be extended to include the simultaneous quantification of more than one target DNA with a single internal control.     

Estimation of the reaction efficiency in polymerase chain reaction

Polymerase chain reaction (PCR) is largely used in molecular biology for increasing the copy number of a specific DNA fragment. The succession of 20 replication cycles makes it possible to multiply the quantity of the fragment of interest by a factor of 1 million. The PCR technique has revolutionized genomics research. Several quantification methodologies are available to determine the DNA replication efficiency of the reaction which is the probability of replication of a DNA molecule at a replication cycle. We elaborate a quantification procedure based on the exponential phase and the early saturation phase of PCR. The reaction efficiency is supposed to be constant in the exponential phase, and decreasing in the saturation phase. We propose to model the PCR amplification process by a branching process which starts as a Galton-Watson branching process followed by a size-dependent process. Using this stochastic modelling and the conditional least-squares estimation method, we infer the reaction efficiency from a single PCR trajectory.