Incorporation of measurement of DNA integrity into qPCR assays

Optimal accuracy of quantitative PCR (qPCR) requires correction for integrity of the target sequence. Here we combine the mathematics of the Poisson distribution and exponential amplification to show that the frequency of lesions per base (which prevent PCR amplification) can be derived from the slope of the regression line between cycle threshold (Ct) and amplicon length. We found that the amplifiable fraction (AF) of a target can be determined from this frequency and the target length. Experimental results from this method correlated with both the magnitude of a damaging agent and with other measures of DNA damage. Applying the method to a reference sequence, we determined the values for lesions/base in control samples, as well as in the AFs of the target sequence in qPCR samples collected from leukemic patients. The AFs used to calculate the final qPCR result were generally >0.5, but were <0.2 in a few samples, indicating significant degradation. We conclude that DNA damage is not always predictable; quantifying the DNA integrity of a sample and determining the AF of a specific qPCR target will improve the accuracy of qPCR and aid in the interpretation of negative results.     

Chemically modified primers for improved multiplex polymerase chain reaction

Multiplex polymerase chain reaction (PCR), the amplification of multiple targets in a single reaction, presents a new set of challenges that further complicate more traditional PCR setups. These complications include a greater probability for nonspecific amplicon formation and for imbalanced amplification of different targets, each of which can compromise quantification and detection of multiple targets. Despite these difficulties, multiplex PCR is frequently used in applications such as pathogen detection, RNA quantification, mutation analysis, and (recently) next generation DNA sequencing. Here we investigated the utility of primers with one or two thermolabile 4-oxo-1-pentyl phosphotriester modifications in improving multiplex PCR performance. Initial endpoint and real-time analyses revealed a decrease in off-target amplification and a subsequent increase in amplicon yield. Furthermore, the use of modified primers in multiplex setups revealed a greater limit of detection and more uniform amplification of each target as compared with unmodified primers. Overall, the thermolabile modified primers present a novel and exciting avenue for improving multiplex PCR performance.     

Quantification of siRNA using competitive qPCR

We have developed a PCR-based short interfering RNA (siRNA) quantification method based on competition between siRNA and a homologous DNA primer for annealing to template DNA, avoiding the requirement for prior conversion of RNA to cDNA. Primers and probe were designed to amplify regions of the human papillomavirus E6 or enhanced green fluorescent protein genes. Having confirmed siRNA could not act as primer for amplicon generation, the lowest competing primer concentration yielding a linear relationship between template DNA amount (0.1–50 ng) and cycle of threshold (Ct) was determined (6.25 nM). Under these conditions addition of sequence-specific siRNA to the competitive quantitative PCR (cqPCR), resulted in a dose-dependent linear increase in Ct value. 2′-O-methyl ribose-modified siRNA retained an ability to inhibit template amplification in serum, unlike unmodified siRNAs that were susceptible to endonucleases. Mismatch-bearing or truncated siRNAs failed to inhibit template amplification confirming sequence specificity and an ability to discriminate between degraded and non-degraded siRNA sequences. Following delivery of E6 siRNA to C33-A cells using oligofectamine or His6 reducible polymers, siRNA uptake was quantified by cqPCR, revealing dose-dependent uptake. We anticipate that cqPCR will allow accurate determination of siRNA pharmacokinetics following in vivo delivery, greatly facilitating development of therapeutic siRNA delivery strategies.     

A method for quantification of absolute amounts of nucleic acids by (RT)–PCR and a new mathematical model for data analysis

Accurate quantification of nucleic acids by competitive (RT)–PCR requires a valid internal standard, a reference for data normalization and an adequate mathematical model for data analysis. We report here an effective procedure for the generation of homologous RNA internal standards and a strategy for synthesizing and using a reference target RNA in quantification of absolute amounts of nucleic acids. Further, a new mathematical model describing the general kinetic features of competitive PCR was developed. The model extends the validity of quantitative competitive (RT)–PCR beyond the exponential phase. The new method eliminates the errors arising from different amplification efficiencies of the co-amplified sequences and from heteroduplex formation in the system. The high accuracy (relative error <2%) is comparable to the recently developed real time detection 5′-nuclease PCR. Also, corresponding computer software has been devised for practical data analysis.     

Mathematics of quantitative kinetic PCR and the application of standard curves

Fluorescent monitoring of DNA amplification is the basis of real-time PCR, from which target DNA concentration can be determined from the fractional cycle at which a threshold amount of amplicon DNA is produced. Absolute quantification can be achieved using a standard curve constructed by amplifying known amounts of target DNA. In this study, the mathematics of quantitative PCR are examined in detail, from which several fundamental aspects of the threshold method and the application of standard curves are illustrated. The construction of five replicate standard curves for two pairs of nested primers was used to examine the reproducibility and degree of quantitative variation using SYBER® Green I fluorescence. Based upon this analysis the application of a single, well- constructed standard curve could provide an estimated precision of ±6–21%, depending on the number of cycles required to reach threshold. A simplified method for absolute quantification is also proposed, in which quantitative scale is determined by DNA mass at threshold.     

Antisense PCR: A simple and robust method for performing nested single-tube PCR

To overcome the disadvantages of two-round nested PCR, we developed a simple and robust closed single-tube nested PCR method (antisense PCR). The method uses antisense oligonucleotides that carry a 5′ tag and that can potentially hybridize to the 3′ ends of the outer primers, depending on the annealing temperature. During initial cycles, which are performed at a high annealing temperature, the antisense oligonucleotides do not hybridize and amplification is directed by the outer primers. During later cycles, for which the annealing temperature is decreased, the outer primers hybridize to the antisense oligonucleotides, extend to produce sequences that are mismatched to the amplicon templates, and consequently become inactivated, whereas the inner primers hybridize to the amplicon templates and continue amplification. Antisense quantitative PCR (qPCR) was compared with one-round qPCR for real-time amplification of four PCR targets (BCR, APC, N-RAS, and a rearranged IGH gene). It had equal amplification efficiency but produced much less nonspecific amplification. Antisense PCR enables both endpoint detection and real-time quantification. It can substitute for two-round nested PCRs but may also be applicable to instances of one-round PCR in which nonspecificity is a problem.     

Assessing the Performance Capabilities of LRE-Based Assays for Absolute Quantitative Real-Time PCR

Background

Linear regression of efficiency or LRE introduced a new paradigm for conducting absolute quantification, which does not require standard curves, can generate absolute accuracies of ±25% and has single molecule sensitivity. Derived from adapting the classic Boltzmann sigmoidal function to PCR, target quantity is calculated directly from the fluorescence readings within the central region of an amplification profile, generating 4–8 determinations from each amplification reaction.

Findings

Based on generating a linear representation of PCR amplification, the highly visual nature of LRE analysis is illustrated by varying reaction volume and amplification efficiency, which also demonstrates how LRE can be used to model PCR. Examining the dynamic range of LRE further demonstrates that quantitative accuracy can be maintained down to a single target molecule, and that target quantification below ten molecules conforms to that predicted by Poisson distribution. Essential to the universality of optical calibration, the fluorescence intensity generated by SYBR Green I (FU/bp) is shown to be independent of GC content and amplicon size, further verifying that absolute scale can be established using a single quantitative standard. Two high-performance lambda amplicons are also introduced that in addition to producing highly precise optical calibrations, can be used as benchmarks for performance testing. The utility of limiting dilution assay for conducting platform-independent absolute quantification is also discussed, along with the utility of defining assay performance in terms of absolute accuracy.

Conclusions

Founded on the ability to exploit lambda gDNA as a universal quantitative standard, LRE provides the ability to conduct absolute quantification using few resources beyond those needed for sample preparation and amplification. Combined with the quantitative and quality control capabilities of LRE, this kinetic-based approach has the potential to fundamentally transform how real-time qPCR is conducted.     

SYBR GreenⅠ实时荧光定量RT-PCR测定肠道病毒71型（EV71）RNA拷贝数方法的建立

目的建立SYBR GreenⅠ荧光染料实时定量RT-PCR方法,测定实验动物等来源的EV71病毒RNA。方法运用EV71VP1保守区引物,优化real time RT-PCR条件,运用NASBA方法扩增EV71病毒RNA,计算拷贝数,经10倍系列稀释做出标准曲线,作为EV71病毒RNA定量检测的外标准品。结果应用Qiagen公司QuantiTect SYBR Green RT-PCR Kit,该标准品可精确定量到100copies/μL,PCR扩增效率达到99.5%。结论 SYBRGreenⅠ荧光染料实时定量PCR法测定EV71病毒RNA拷贝数的方法敏感性高、稳定性好,可用于EV71病毒RNA载量的定量测定。     

The effect of low-level activation on the mechanical properties of isolated frog muscle fibers

The mechanical properties, as revealed by minute length changes, of isolated twitch fibers of the frog have been studied at rest and during low-level activation. Resting tension is 77 ± 23 mN/cm² (mean ± SD) at 2.2 µm sarcomere length.¹ The slope of the tension curve (ΔP/ΔL) recorded during a constant-speed length change of a resting fiber is initially large. At length changes exceeding about 0.18 % of the initial length of the fiber ΔP/ΔL falls abruptly and remains close to zero during the rest of the length change. The amplitude of the tension response is reduced after a length change and returns to normal in about 3 min. Hypertonic sucrose-Ringer solutions cause a small, maintained rise in tension up to 1.4–1.6 times normal osmotic strength. Higher sucrose concentrations cause relatively large, transient tension responses. The initial ΔP/ΔL is increased in moderately hypertonic solutions; it may be reduced in more strongly hypertonic solutions. Elevated [K]_o (range 10–17.5 mM) causes a marked reduction in ΔP/ΔL. In this range of [K]_o the reduction is not accompanied by changes in resting tension. Addition of 1–1.5 mM caffeine to the Ringer solution affects the resting tension very little but also reduces ΔP/ΔL. The results suggest that stiffness and tension development are not related in a simple way.     

Telomerase-dependent repeat divergence at the 3' ends of yeast telomeres

Yeast telomeres consist of ~300 nt of degenerate repeats with the consensus sequence G_2–3(TG)_1–6. We developed a method for the amplification of a genetically marked telomere by PCR, allowing precise length and sequence determination of the G-rich strand including the 3′ terminus. We examined wild-type cells, telomerase RNA deficient cells and a strain deleted for YKU70, which encodes for a protein involved in telomere maintenance and DNA double strand break repair. The 3′ end of the G-rich strand was found to be at a variable position within the telomeric repeat. No preference for either thymine or guanine as the 3′ base was detected. Comparison of telomere sequences from clonal populations revealed that telomeres consist of a centromere-proximal region of stable sequence and a distal region with differing degenerate repeats. In wild-type as well as yku70-Δ cells, variation in the degenerate telomeric repeats was detected starting 40–100 nt from the 3′ end. Sequence divergence was abolished after deletion of the telomerase RNA gene. Thus, this region defines the domain where telomere shortening and telomerase-mediated extension occurs. Since this domain is much larger than the number of nucleo­tides lost per generation in the absence of telomerase, we propose that telomerase does not extend a given telomere in every cell cycle.     

Predicting loop–helix tertiary structural contacts in RNA pseudoknots

Tertiary interactions between loops and helical stems play critical roles in the biological function of many RNA pseudoknots. However, quantitative predictions for RNA tertiary interactions remain elusive. Here we report a statistical mechanical model for the prediction of noncanonical loop–stem base-pairing interactions in RNA pseudoknots. Central to the model is the evaluation of the conformational entropy for the pseudoknotted folds with defined loop–stem tertiary structural contacts. We develop an RNA virtual bond-based conformational model (Vfold model), which permits a rigorous computation of the conformational entropy for a given fold that contains loop–stem tertiary contacts. With the entropy parameters predicted from the Vfold model and the energy parameters for the tertiary contacts as inserted parameters, we can then predict the RNA folding thermodynamics, from which we can extract the tertiary contact thermodynamic parameters from theory–experimental comparisons. These comparisons reveal a contact enthalpy (ΔH) of −14 kcal/mol and a contact entropy (ΔS) of −38 cal/mol/K for a protonated C⁺•(G–C) base triple at pH 7.0, and (ΔH = −7 kcal/mol, ΔS = −19 cal/mol/K) for an unprotonated base triple. Tests of the model for a series of pseudoknots show good theory–experiment agreement. Based on the extracted energy parameters for the tertiary structural contacts, the model enables predictions for the structure, stability, and folding pathways for RNA pseudoknots with known or postulated loop–stem tertiary contacts from the nucleotide sequence alone.     

PCR Bias in Ecological Analysis: a Case Study for Quantitative Taq Nuclease Assays in Analyses of Microbial Communities

Succession of ecotypes, physiologically diverse strains with negligible rRNA sequence divergence, may explain the dominance of small, red-pigmented (phycoerythrin-rich) cyanobacteria in the autotrophic picoplankton of deep lakes (C. Postius and A. Ernst, Arch. Microbiol. 172:69–75, 1999). In order to test this hypothesis, it is necessary to determine the abundance of specific ecotypes or genotypes in a mixed background of phylogenetically similar organisms. In this study, we examined the performance of Taq nuclease assays (TNAs), PCR-based assays in which the amount of an amplicon is monitored by hydrolysis of a labeled oligonucleotide (TaqMan probe) when hybridized to the amplicon. High accuracy and a 7-order detection range made the real-time TNA superior to the corresponding end point technique. However, in samples containing mixtures of homologous target sequences, quantification can be biased due to limited specificity of PCR primers and probe oligonucleotides and due to accumulation of amplicons that are not detected by the TaqMan probe. A decrease in reaction efficiency, which can be recognized by direct monitoring of amplification, provides experimental evidence for the presence of such a problem and emphasizes the need for real-time technology in quantitative PCR. Use of specific primers and probes and control of amplification efficiency allow correct quantification of target DNA in the presence of an up to 10⁴-fold excess of phylogenetically similar DNA and of an up to 10⁷-fold excess of dissimilar DNA.     

Performance of ultra-deep pyrosequencing in analysis of HIV-1 pol gene variation

Introduction

Ultra-deep pyrosequencing (UDPS) has been used to detect minority variants within HIV-1 populations. Some aspects of the quality and reproducibility of UDPS have been previously evaluated, but comprehensive studies are still needed.

Principal Finding

In this study the UDPS technology (FLX platform) was evaluated by analyzing a 120 base pair fragment of the HIV-1 pol gene from plasma samples from two patients and artificial mixtures of molecular clones. UDPS was performed using an optimized experimental protocol and an in-house data cleaning strategy. Nine samples and mixtures were analyzed and the average number of reads per sample was 19,404 (range 8,858–26,846). The two patient plasma samples were analyzed twice and quantification of viral variants was found to be highly repeatable for variants representing >0.27% of the virus population, whereas some variants representing 0.11–0.27% were detected in only one of the two UDPS runs. Bland-Altman analysis showed that a repeated measurement would have a 95% likelihood to lie approximately within ±0.5 log₁₀ of the initial estimate. A similar level of agreement was observed for variant frequency estimates in forward vs. reverse sequencing direction, but here the agreement was higher for common variants than for rare variants. UDPS following PCR amplification with alternative primers indicated that some variants may be incorrectly quantified due to primer-related selective amplification. Finally, the in vitro recombination rate during PCR was evaluated using artificial mixtures of clones and was found to be low. The most abundant in vitro recombinant represented 0.25% of all UDPS reads.

Conclusion

This study demonstrates that this UDPS protocol results in low experimental noise and high repeatability, which is relevant for future research and clinical use of the UDPS technology. The low rate of in vitro recombination suggests that this UDPS system can be used to study genetic variants and mutational linkage.     

Evaluation of PCR Amplification Bias by Terminal Restriction Fragment Length Polymorphism Analysis of Small-Subunit rRNA and mcrA Genes by Using Defined Template Mixtures of Methanogenic Pure Cultures and Soil DNA Extracts

Terminal restriction fragment length polymorphism (T-RFLP) analysis is a widely used method for profiling microbial community structure in different habitats by targeting small-subunit (SSU) rRNA and also functional marker genes. It is not known, however, whether relative gene frequencies of individual community members are adequately represented in post-PCR amplicon frequencies as shown by T-RFLP. In this study, precisely defined artificial template mixtures containing genomic DNA of four different methanogens in various ratios were prepared for subsequent T-RFLP analysis. PCR amplicons were generated from defined mixtures targeting not only the SSU rRNA but also the methyl-coenzyme M reductase (mcrA/mrtA) genes of methanogens. Relative amplicon frequencies of microorganisms were quantified by comparing fluorescence intensities of characteristic terminal restriction fragments. SSU ribosomal DNA (rDNA) template ratios in defined template mixtures of the four-membered community were recovered absolutely by PCR-T-RFLP analysis, which demonstrates that the T-RFLP analysis evaluated can give a quantitative view of the template pool. SSU rDNA-targeted T-RFLP analysis of a natural community was found to be highly reproducible, independent of PCR annealing temperature, and unaffected by increasing PCR cycle numbers. Ratios of mcrA-targeted T-RFLP analysis were biased, most likely by PCR selection due to the degeneracy of the primers used. Consequently, for microbial community analyses, each primer system used should be evaluated carefully for possible PCR bias. In fact, such bias can be detected by using T-RFLP analysis as a tool for the precise quantification of the PCR product pool.     

The effects of low levels of dystrophin on mouse muscle function and pathology

Duchenne muscular dystrophy (DMD) is a severe progressive muscular disorder caused by reading frame disrupting mutations in the DMD gene, preventing the synthesis of functional dystrophin. As dystrophin provides muscle fiber stability during contractions, dystrophin negative fibers are prone to exercise-induced damage. Upon exhaustion of the regenerative capacity, fibers will be replaced by fibrotic and fat tissue resulting in a progressive loss of function eventually leading to death in the early thirties. With several promising approaches for the treatment of DMD aiming at dystrophin restoration in clinical trials, there is an increasing need to determine more precisely which dystrophin levels are sufficient to restore muscle fiber integrity, protect against muscle damage and improve muscle function.To address this we generated a new mouse model (mdx-Xist ^Δhs) with varying, low dystrophin levels (3–47%, mean 22.7%, stdev 12.1, n = 24) due to skewed X-inactivation. Longitudinal sections revealed that within individual fibers, some nuclei did and some did not express dystrophin, resulting in a random, mosaic pattern of dystrophin expression within fibers. Mdx-Xist ^Δhs, mdx and wild type females underwent a 12 week functional test regime consisting of different tests to assess muscle function at base line, or after chronic treadmill running exercise. Overall, mdx-Xist ^Δhs mice with 3–14% dystrophin outperformed mdx mice in the functional tests. Improved histopathology was observed in mice with 15–29% dystrophin and these levels also resulted in normalized expression of pro-inflammatory biomarker genes, while for other parameters >30% of dystrophin was needed. Chronic exercise clearly worsened pathology, which needed dystrophin levels >20% for protection. Based on these findings, we conclude that while even dystrophin levels below 15% can improve pathology and performance, levels of >20% are needed to fully protect muscle fibers from exercise-induced damage.     

Highly parallel genome-wide expression analysis of single mammalian cells

Background

We have developed a high-throughput amplification method for generating robust gene expression profiles using single cell or low RNA inputs.

Methodology/Principal Findings

The method uses tagged priming and template-switching, resulting in the incorporation of universal PCR priming sites at both ends of the synthesized cDNA for global PCR amplification. Coupled with a whole-genome gene expression microarray platform, we routinely obtain expression correlation values of R²∼0.76–0.80 between individual cells and R²∼0.69 between 50 pg total RNA replicates. Expression profiles generated from single cells or 50 pg total RNA correlate well with that generated with higher input (1 ng total RNA) (R²∼0.80). Also, the assay is sufficiently sensitive to detect, in a single cell, approximately 63% of the number of genes detected with 1 ng input, with approximately 97% of the genes detected in the single-cell input also detected in the higher input.

Conclusions/Significance

In summary, our method facilitates whole-genome gene expression profiling in contexts where starting material is extremely limiting, particularly in areas such as the study of progenitor cells in early development and tumor stem cell biology.