A kinetic-based sigmoidal model for the polymerase chain reaction and its application to high-capacity absolute quantitative real-time PCR

Background

Based upon defining a common reference point, current real-time quantitative PCR technologies compare relative differences in amplification profile position. As such, absolute quantification requires construction of target-specific standard curves that are highly resource intensive and prone to introducing quantitative errors. Sigmoidal modeling using nonlinear regression has previously demonstrated that absolute quantification can be accomplished without standard curves; however, quantitative errors caused by distortions within the plateau phase have impeded effective implementation of this alternative approach.

Results

Recognition that amplification rate is linearly correlated to amplicon quantity led to the derivation of two sigmoid functions that allow target quantification via linear regression analysis. In addition to circumventing quantitative errors produced by plateau distortions, this approach allows the amplification efficiency within individual amplification reactions to be determined. Absolute quantification is accomplished by first converting individual fluorescence readings into target quantity expressed in fluorescence units, followed by conversion into the number of target molecules via optical calibration. Founded upon expressing reaction fluorescence in relation to amplicon DNA mass, a seminal element of this study was to implement optical calibration using lambda gDNA as a universal quantitative standard. Not only does this eliminate the need to prepare target-specific quantitative standards, it relegates establishment of quantitative scale to a single, highly defined entity. The quantitative competency of this approach was assessed by exploiting "limiting dilution assay" for absolute quantification, which provided an independent gold standard from which to verify quantitative accuracy. This yielded substantive corroborating evidence that absolute accuracies of ± 25% can be routinely achieved. Comparison with the LinReg and Miner automated qPCR data processing packages further demonstrated the superior performance of this kinetic-based methodology.

Conclusion

Called "linear regression of efficiency" or LRE, this novel kinetic approach confers the ability to conduct high-capacity absolute quantification with unprecedented quality control capabilities. The computational simplicity and recursive nature of LRE quantification also makes it amenable to software implementation, as demonstrated by a prototypic Java program that automates data analysis. This in turn introduces the prospect of conducting absolute quantification with little additional effort beyond that required for the preparation of the amplification reactions.     

Real-Time Relative qPCR without Reference to Control Samples and Estimation of Run-Specific PCR Parameters from Run-Internal Mini-Standard Curves

Background

Real-Time quantitative PCR is an important tool in research and clinical settings. Here, we describe two new approaches that broaden the scope of real-time quantitative PCR; namely, run-internal mini standard curves (RIMS) and direct real-time relative quantitative PCR (drqPCR). RIMS are an efficient alternative to traditional standard curves and provide both run-specific and target-specific estimates of PCR parameters. The drqPCR enables direct estimation of target ratios without reference to conventional control samples.

Methodology/Principal Findings

In this study, we compared RIMS-based drqPCR with classical quantifications based on external standard curves and the “comparative Ct method”. Specifically, we used a raw real-time PCR dataset as the basis for more than two-and-a-half million simulated quantifications with various user-defined conditions. Compared with classical approaches, we found that RIMS-based drqPCR provided superior precision and comparable accuracy.

Conclusions/Significance

The obviation of referencing to control samples is attractive whenever unpaired samples are quantified. This may be in clinical and research settings; for instance, studies on chimerism, TREC quantifications, copy number variations etc. Also, lab-to-lab comparability can be greatly simplified.     

A Mechanistic Model of PCR for Accurate Quantification of Quantitative PCR Data

Background

Quantitative PCR (qPCR) is a workhorse laboratory technique for measuring the concentration of a target DNA sequence with high accuracy over a wide dynamic range. The gold standard method for estimating DNA concentrations via qPCR is quantification cycle () standard curve quantification, which requires the time- and labor-intensive construction of a standard curve. In theory, the shape of a qPCR data curve can be used to directly quantify DNA concentration by fitting a model to data; however, current empirical model-based quantification methods are not as reliable as standard curve quantification.

Principal Findings

We have developed a two-parameter mass action kinetic model of PCR (MAK2) that can be fitted to qPCR data in order to quantify target concentration from a single qPCR assay. To compare the accuracy of MAK2-fitting to other qPCR quantification methods, we have applied quantification methods to qPCR dilution series data generated in three independent laboratories using different target sequences. Quantification accuracy was assessed by analyzing the reliability of concentration predictions for targets at known concentrations. Our results indicate that quantification by MAK2-fitting is as reliable as standard curve quantification for a variety of DNA targets and a wide range of concentrations.

Significance

We anticipate that MAK2 quantification will have a profound effect on the way qPCR experiments are designed and analyzed. In particular, MAK2 enables accurate quantification of portable qPCR assays with limited sample throughput, where construction of a standard curve is impractical.     

A Novel Multiplex Tetra-Primer ARMS-PCR for the Simultaneous Genotyping of Six Single Nucleotide Polymorphisms Associated with Female Cancers

Background

The tetra-primer amplification refractory mutation system PCR (T-ARMS-PCR) is a fast and economical means of assaying SNP''s, requiring only PCR amplification and subsequent electrophoresis for the determination of genotypes. To improve the throughput and efficiency of T-ARMS-PCR, we combined T-ARMS-PCR with a chimeric primer-based temperature switch PCR (TSP) strategy, and used capillary electrophoresis (CE) for amplicon separation and identification. We assessed this process in the simultaneous genotyping of four breast cancer–and two cervical cancer risk–related SNPs.

Methods

A total of 24 T-ARMS-PCR primers, each 5′-tagged with a universal sequence and a pair of universal primers, were pooled together to amplify the 12 target alleles of 6 SNPs in 186 control female blood samples. Direct sequencing of all samples was also performed to assess the accuracy of this method.

Results

Of the 186 samples, as many as 11 amplicons can be produced in one single PCR and separated by CE. Genotyping results of the multiplex T-ARMS-PCR were in complete agreement with direct sequencing of all samples.

Conclusions

This novel multiplex T-ARMS-PCR method is the first reported method allowing one to genotype six SNPs in a single reaction with no post-PCR treatment other than electrophoresis. This method is reliable, fast, and easy to perform.     

Clinical Benefit of Liver Stiffness Measurement at 3 Months after Kasai Hepatoportoenterostomy to Predict the Liver Related Events in Biliary Atresia

Background

The progression of hepatic fibrosis may result in decompensated hepatic failure with cirrhosis, liver related events (LRE) such as ascites, variceal bleeding, and death after successful and timely Kasai hepatoportoenterostomy (HPE) in biliary atresia. The aim of this study is to suggest clinical benefit of the liver stiffness measurement (LSM) using transient elastography at 3 months after the Kasai operation to predict LRE.

Methods

Between January 2007 and December 2011, 69 eligible biliary atresia patients who underwent Kasai HPE and performed transient elastography before and 3 months after HPE were included. The occurrences of LRE were analyzed for all patients. All patients were divided into 2 groups (with and without LRE) for comparison. Multivariate analysis was used to detect the independent risk factors of LRE. The area under the receiver operation characteristics curve (AUROC) was used to establish the LSM optimal cutoff value of 3 months after Kasai operation in predicting LRE.

Results

LSM value, aminotransferase, albumin, bilirubin, and PT-INR significantly differed among the two groups. Multivariate analysis demonstrated LSM value as the most powerful independent factor of the development of LRE. The cut-off value of 19.9 kPa was calculated to be optimal for predicting LRE development with total sensitivity and specificity of 1.804. AUROC resulted in 0.943, with sensitivity of 85.3% and specificity of 95.2%.

Conclusions

The LSM value of 3 months after Kasai HPE can be a useful predictor of LRE development.     

Reducing amplification artifacts in high multiplex amplicon sequencing by using molecular barcodes

Background

PCR amplicon sequencing has been widely used as a targeted approach for both DNA and RNA sequence analysis. High multiplex PCR has further enabled the enrichment of hundreds of amplicons in one simple reaction. At the same time, the performance of PCR amplicon sequencing can be negatively affected by issues such as high duplicate reads, polymerase artifacts and PCR amplification bias. Recently researchers have made some good progress in addressing these shortcomings by incorporating molecular barcodes into PCR primer design. So far, most work has been demonstrated using one to a few pairs of primers, which limits the size of the region one can analyze.

Results

We developed a simple protocol, which enables the use of molecular barcodes in high multiplex PCR with hundreds of amplicons. Using this protocol and reference materials, we demonstrated the applications in accurate variant calling at very low fraction over a large region and in targeted RNA quantification. We also evaluated the protocol’s utility in profiling FFPE samples.

Conclusions

We demonstrated the successful implementation of molecular barcodes in high multiplex PCR, with multiplex scale many times higher than earlier work. We showed that the new protocol combines the benefits of both high multiplex PCR and molecular barcodes, i.e. the analysis of a very large region, low DNA input requirement, very good reproducibility and the ability to detect as low as 1 % mutations with minimal false positives (FP).

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1806-8) contains supplementary material, which is available to authorized users.     

A statistical approach to quantification of genetically modified organisms (GMO) using frequency distributions

Background

According to Regulation (EU) No 619/2011, trace amounts of non-authorised genetically modified organisms (GMO) in feed are tolerated within the EU if certain prerequisites are met. Tolerable traces must not exceed the so-called ‘minimum required performance limit’ (MRPL), which was defined according to the mentioned regulation to correspond to 0.1% mass fraction per ingredient. Therefore, not yet authorised GMO (and some GMO whose approvals have expired) have to be quantified at very low level following the qualitative detection in genomic DNA extracted from feed samples. As the results of quantitative analysis can imply severe legal and financial consequences for producers or distributors of feed, the quantification results need to be utterly reliable.

Results

We developed a statistical approach to investigate the experimental measurement variability within one 96-well PCR plate. This approach visualises the frequency distribution as zygosity-corrected relative content of genetically modified material resulting from different combinations of transgene and reference gene Cq values. One application of it is the simulation of the consequences of varying parameters on measurement results. Parameters could be for example replicate numbers or baseline and threshold settings, measurement results could be for example median (class) and relative standard deviation (RSD). All calculations can be done using the built-in functions of Excel without any need for programming. The developed Excel spreadsheets are available (see section ‘Availability of supporting data’ for details). In most cases, the combination of four PCR replicates for each of the two DNA isolations already resulted in a relative standard deviation of 15% or less.

Conclusions

The aims of the study are scientifically based suggestions for minimisation of uncertainty of measurement especially in —but not limited to— the field of GMO quantification at low concentration levels. Four PCR replicates for each of the two DNA isolations seem to be a reasonable minimum number to narrow down the possible spread of results.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-014-0407-x) contains supplementary material, which is available to authorized users.     

Simultaneous Detection and Differentiation of Human Papillomavirus Genotypes 6, 11, 16 and 18 by AllGlo Quadruplex Quantitative PCR

Background

Human papillomaviruses (HPV) are classified into high-risk HPV and low-risk HPV. The most common high-risk HPV types in cervical cancer are HPV 16 and 18, and the most common low-risk types causing genital warts are HPV 6 and HPV 11. In this study, applying novel AllGlo fluorescent probes, we established a quadruplex quantitative PCR method to simultaneously detect and differentiate HPV 6, 11, 16 and 18 in a single tube.

Methods

The specificity, the sensitivity, the detection limit, the reproducibility and the standard curve of this method were examined. Finally, clinical samples that had been tested previously by TaqMan PCR and HPV GenoArray (GA) test were used to verify the accuracy and sensitivity of the method.

Results

The assay has a sensitivity of 10¹ to 10² copies/test and a linear detection range from 10¹ to 10⁸ copies/test. The mean amplification efficiencies for HPV 6, 11, 16, and 18 were 0.97, 1.10, 0.93 and 1.20, respectively, and the mean correlation coefficient (r²) of each standard curve was above 0.99 for plasmid templates ranging from 10³ to 10⁷ copies/test. There was 100% agreement between the AllGlo quadruplex quantitative PCR, HPV GA test and TaqMan uniplex qPCR methods.

Conclusions

AllGlo quadruplex quantitative PCR in a single tube has the advantages of relatively high throughput, good reproducibility, high sensitivity, high specificity, and a wide linear range of detection. The convenient single tube format makes this assay a powerful tool for the studies of mixed infections by multiple pathogens, viral typing and viral load quantification.     

Quantification of epigenetic biomarkers: an evaluation of established and emerging methods for DNA methylation analysis

Background

DNA methylation is an important epigenetic mechanism in several human diseases, most notably cancer. The quantitative analysis of DNA methylation patterns has the potential to serve as diagnostic and prognostic biomarkers, however, there is currently a lack of consensus regarding the optimal methodologies to quantify methylation status. To address this issue we compared five analytical methods: (i) MethyLight qPCR, (ii) MethyLight digital PCR (dPCR), methylation-sensitive and -dependent restriction enzyme (MSRE/MDRE) digestion followed by (iii) qPCR or (iv) dPCR, and (v) bisulfite amplicon next generation sequencing (NGS). The techniques were evaluated for linearity, accuracy and precision.

Results

MethyLight qPCR displayed the best linearity across the range of tested samples. Observed methylation measured by MethyLight- and MSRE/MDRE-qPCR and -dPCR were not significantly different to expected values whilst bisulfite amplicon NGS analysis over-estimated methylation content. Bisulfite amplicon NGS showed good precision, whilst the lower precision of qPCR and dPCR analysis precluded discrimination of differences of < 25% in methylation status. A novel dPCR MethyLight assay is also described as a potential method for absolute quantification that simultaneously measures both sense and antisense DNA strands following bisulfite treatment.

Conclusions

Our findings comprise a comprehensive benchmark for the quantitative accuracy of key methods for methylation analysis and demonstrate their applicability to the quantification of circulating tumour DNA biomarkers by using sample concentrations that are representative of typical clinical isolates.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1174) contains supplementary material, which is available to authorized users.     

Application of quantitative proteomics to the integrated analysis of the ubiquitylated and global proteomes of xenograft tumor tissues

Background

Post-translational modification by ubiquitin is a fundamental regulatory mechanism that is implicated in many cellular processes including the cell cycle, apoptosis, cell adhesion, angiogenesis, and tumor growth. The low stoichiometry of ubiquitylation presents an analytical challenge for the detection of endogenously modified proteins in the absence of enrichment strategies. The recent availability of antibodies recognizing peptides with Lys residues containing a di-Gly ubiquitin remnant (K-ε-GG) has greatly improved the ability to enrich and identify ubiquitylation sites from complex protein lysates via mass spectrometry. To date, there have not been any published studies that quantitatively assess the changes in endogenous ubiquitin-modification protein stoichiometry status at the proteome level from different tissues.

Results

In this study, we applied an integrated quantitative mass spectrometry based approach using isobaric tags for relative and absolute quantitation (iTRAQ) to interrogate the ubiquitin-modified proteome and the cognate global proteome levels from luminal and basal breast cancer patient-derived xenograft tissues. Among the proteins with quantitative global and ubiquitylation data, 91 % had unchanged levels of total protein relative abundance, and less than 5 % of these proteins had up- or down-regulated ubiquitylation levels. Of particular note, greater than half of the proteins with observed changes in their total protein level also had up- or down-regulated changes in their ubiquitylation level.

Conclusions

This is the first report of the application of iTRAQ-based quantification to the integrated analysis of the ubiquitylated and global proteomes at the tissue level. Our results underscore the importance of conducting integrated analyses of the global and ubiquitylated proteomes toward elucidating the specific functional significance of ubiquitylation.

Electronic supplementary material

The online version of this article (doi:10.1186/s12014-015-9086-5) contains supplementary material, which is available to authorized users.     

An Efficient Multistrategy DNA Decontamination Procedure of PCR Reagents for Hypersensitive PCR Applications

Background

PCR amplification of minute quantities of degraded DNA for ancient DNA research, forensic analyses, wildlife studies and ultrasensitive diagnostics is often hampered by contamination problems. The extent of these problems is inversely related to DNA concentration and target fragment size and concern (i) sample contamination, (ii) laboratory surface contamination, (iii) carry-over contamination, and (iv) contamination of reagents.

Methodology/Principal Findings

Here we performed a quantitative evaluation of current decontamination methods for these last three sources of contamination, and developed a new procedure to eliminate contaminating DNA contained in PCR reagents. We observed that most current decontamination methods are either not efficient enough to degrade short contaminating DNA molecules, rendered inefficient by the reagents themselves, or interfere with the PCR when used at doses high enough to eliminate these molecules. We also show that efficient reagent decontamination can be achieved by using a combination of treatments adapted to different reagent categories. Our procedure involves γ- and UV-irradiation and treatment with a mutant recombinant heat-labile double-strand specific DNase from the Antarctic shrimp Pandalus borealis. Optimal performance of these treatments is achieved in narrow experimental conditions that have been precisely analyzed and defined herein.

Conclusions/Significance

There is not a single decontamination method valid for all possible contamination sources occurring in PCR reagents and in the molecular biology laboratory and most common decontamination methods are not efficient enough to decontaminate short DNA fragments of low concentration. We developed a versatile multistrategy decontamination procedure for PCR reagents. We demonstrate that this procedure allows efficient reagent decontamination while preserving the efficiency of PCR amplification of minute quantities of DNA.     

Optimizing Taq Polymerase Concentration for Improved Signal-to-Noise in the Broad Range Detection of Low Abundance Bacteria

Background

PCR in principle can detect a single target molecule in a reaction mixture. Contaminating bacterial DNA in reagents creates a practical limit on the use of PCR to detect dilute bacterial DNA in environmental or public health samples. The most pernicious source of contamination is microbial DNA in DNA polymerase preparations. Importantly, all commercial Taq polymerase preparations inevitably contain contaminating microbial DNA. Removal of DNA from an enzyme preparation is problematical.

Methodology/Principal Findings

This report demonstrates that the background of contaminating DNA detected by quantitative PCR with broad host range primers can be decreased greater than 10-fold through the simple expedient of Taq enzyme dilution, without altering detection of target microbes in samples. The general method is: For any thermostable polymerase used for high-sensitivity detection, do a dilution series of the polymerase crossed with a dilution series of DNA or bacteria that work well with the test primers. For further work use the concentration of polymerase that gave the least signal in its negative control (H₂O) while also not changing the threshold cycle for dilutions of spiked DNA or bacteria compared to higher concentrations of Taq polymerase.

Conclusions/Significance

It is clear from the studies shown in this report that a straightforward procedure of optimizing the Taq polymerase concentration achieved “treatment-free” attenuation of interference by contaminating bacterial DNA in Taq polymerase preparations. This procedure should facilitate detection and quantification with broad host range primers of a small number of bona fide bacteria (as few as one) in a sample.     

DAG Expression: High-Throughput Gene Expression Analysis of Real-Time PCR Data Using Standard Curves for Relative Quantification

Background

Real-time quantitative PCR (qPCR) is still the gold-standard technique for gene-expression quantification. Recent technological advances of this method allow for the high-throughput gene-expression analysis, without the limitations of sample space and reagent used. However, non-commercial and user-friendly software for the management and analysis of these data is not available.

Results

The recently developed commercial microarrays allow for the drawing of standard curves of multiple assays using the same n-fold diluted samples. Data Analysis Gene (DAG) Expression software has been developed to perform high-throughput gene-expression data analysis using standard curves for relative quantification and one or multiple reference genes for sample normalization. We discuss the application of DAG Expression in the analysis of data from an experiment performed with Fluidigm technology, in which 48 genes and 115 samples were measured. Furthermore, the quality of our analysis was tested and compared with other available methods.

Conclusions

DAG Expression is a freely available software that permits the automated analysis and visualization of high-throughput qPCR. A detailed manual and a demo-experiment are provided within the DAG Expression software at http://www.dagexpression.com/dage.zip.     

Quantification and Dynamic Monitoring of EGFR T790M in Plasma Cell-Free DNA by Digital PCR for Prognosis of EGFR-TKI Treatment in Advanced NSCLC

Background

Among advanced non-small cell lung cancer (NSCLC) patients with an acquired resistance to epidermal growth factor receptor-tyrosine kinase inhibitors (EGFR-TKI), about 50% carry the T790M mutation, but this frequency in EGFR-TKI-naïve patients and dynamic change during therapy remains unclear. This study investigated the quantification and dynamic change of T790M mutation in plasma cell-free DNA (cf-DNA) of advanced NSCLC patients to assess the clinical outcomes of EGFR-TKI therapy.

Materials and Methods

We retrospectively investigated 135 patients with advanced NSCLC who obtained progression-free survival (PFS) after EGFR-TKI for >6 months for their EGFR sensitive mutations and T790M mutation in matched pre- and post-TKI plasma samples, using denaturing high-performance liquid chromatography (DHPLC), amplification refractory mutation system (ARMS), and digital-PCR (D-PCR). Real-time PCR was performed to measure c-MET amplification.

Results

Detection limit of D-PCR in assessing the T790M mutation was approximately 0.03%. D-PCR identified higher frequency of T790M than ARMS in pre-TKI (31.3% vs. 5.5%) and post-TKI (43.0% vs. 25.2%) plasma samples. Patients with pre-TKI T790M showed inferior PFS (8.9 vs. 12.1 months, p = 0.007) and overall survival (OS, 19.3 vs. 31.9 months, p = 0.001) compared with those without T790M. In patients harboring EGFR sensitive mutation, high quantities of pre-TKI T790M predicted poorer PFS (p = 0.001) on EGFR-TKI than low ones. Moreover, patients who experienced increased quantity of T790M during EGFR-TKI treatment showed superior PFS and OS compared with those with decreased changes (p = 0.044 and p = 0.015, respectively).

Conclusion

Qualitative and quantitative T790M in plasma cf-DNA by D-PCR provided a non-invasive and sensitive assay to predict EGFR-TKI prognosis.     

Quality of DNA Extracted from Mouthwashes

Background

A cost effective, safe and efficient method of obtaining DNA samples is essential in large scale genetic analyses. Buccal cells are an attractive source of DNA, as their collection is non-invasive and can be carried out by mail. However, little attention has been given to the quality of DNA extracted from mouthwashes.

Methodology

Mouthwash-derived DNA was extracted from 500 subjects participating in a genetic study of high myopia. DNA quality was investigated using two standard techniques: agarose gel electrophoresis and quantitative polymerase chain reaction (qPCR).

Principal Findings

Whereas the majority of mouthwash-derived DNA samples showed a single band of high molecular weight DNA by gel electrophoresis, 8.9% (95% CI: 7.1–10.7%) of samples contained only a smear of low-to-medium molecular weight, degraded DNA. The odds of DNA degradation in a subject''s second mouthwash sample, given degradation of the first, was significantly greater than one (OR = 3.13; 95% CI: 1.22–7.39; Fisher''s test P = 0.009), suggesting that DNA degradation was at least partially a subject-specific phenomenon. Approximately 12.4% (95% CI: 10.4–14.4%) of mouthwash-derived DNA failed to PCR amplify efficiently (using an ∼200 bp microsatellite marker). However, we found there was no significant difference in amplification success rate between DNA samples judged to be degraded or non-degraded by gel electrophoresis (Fisher''s test P = 0.5).

Conclusions

This study demonstrated that DNA degradation affects a significant minority of saline mouthwashes, and that the phenomenon is partially subject-specific. Whilst the level of degradation did not significantly prevent successful amplification of short PCR fragments, previous studies suggest that such DNA degradation would compromise more demanding applications.