Impact of variance components on reliability of absolute quantification using digital PCR

Background

Digital polymerase chain reaction (dPCR) is an increasingly popular technology for detecting and quantifying target nucleic acids. Its advertised strength is high precision absolute quantification without needing reference curves. The standard data analytic approach follows a seemingly straightforward theoretical framework but ignores sources of variation in the data generating process. These stem from both technical and biological factors, where we distinguish features that are 1) hard-wired in the equipment, 2) user-dependent and 3) provided by manufacturers but may be adapted by the user. The impact of the corresponding variance components on the accuracy and precision of target concentration estimators presented in the literature is studied through simulation.

Results

We reveal how system-specific technical factors influence accuracy as well as precision of concentration estimates. We find that a well-chosen sample dilution level and modifiable settings such as the fluorescence cut-off for target copy detection have a substantial impact on reliability and can be adapted to the sample analysed in ways that matter. User-dependent technical variation, including pipette inaccuracy and specific sources of sample heterogeneity, leads to a steep increase in uncertainty of estimated concentrations. Users can discover this through replicate experiments and derived variance estimation. Finally, the detection performance can be improved by optimizing the fluorescence intensity cut point as suboptimal thresholds reduce the accuracy of concentration estimates considerably.

Conclusions

Like any other technology, dPCR is subject to variation induced by natural perturbations, systematic settings as well as user-dependent protocols. Corresponding uncertainty may be controlled with an adapted experimental design. Our findings point to modifiable key sources of uncertainty that form an important starting point for the development of guidelines on dPCR design and data analysis with correct precision bounds. Besides clever choices of sample dilution levels, experiment-specific tuning of machine settings can greatly improve results. Well-chosen data-driven fluorescence intensity thresholds in particular result in major improvements in target presence detection. We call on manufacturers to provide sufficiently detailed output data that allows users to maximize the potential of the method in their setting and obtain high precision and accuracy for their experiments.

Electronic supplementary material

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

Docosahexaenoic Acid Supplementation Does Not Improve Western Diet-Induced Cardiomyopathy in Rats

Obesity increases risk for cardiomyopathy in the absence of hypertension, diabetes or ischemia. The fatty acid milieu, modulated by diet, may modify myocardial structure and function, lending partial explanation for the array of cardiomyopathic phenotypy. We sought to identify gross, cellular and ultrastructural myocardial changes associated with Western diet intake, and subsequent modification with docosahexaenoic acid (DHA) supplementation. Wistar and Sprague-Dawley (SD) rats received 1 of 3 diets: control (CON); Western (WES); Western + DHA (WES+DHA). After 12 weeks of treatment, echocardiography was performed and myocardial adiponectin, fatty acids, collagen, area occupied by lipid and myocytes, and ultrastructure were determined. Strain effects included higher serum adiponectin in Wistar rats, and differences in myocardial fatty acid composition. Diet effects were evident in that both WES and WES+DHA feeding were associated with similarly increased left ventricular (LV) diastolic cranial wall thickness (LVW_cr/d) and decreased diastolic internal diameter (LVID_d), compared to CON. Unexpectedly, WES+DHA feeding was associated additionally with increased thickness of the LV cranial wall during systole (LVW_cr/s) and the caudal wall during diastole (LVW_ca/d) compared to CON; this was observed concomitantly with increased serum and myocardial adiponectin. Diastolic dysfunction was present in WES+DHA rats compared to both WES and CON. Myocyte cross sectional area (CSA) was greater in WES compared to CON rats. In both fat-fed groups, transmission electron microscopy (TEM) revealed myofibril degeneration, disorganized mitochondrial cristae, lipid inclusions and vacuolation. In the absence of hypertension and whole body insulin resistance, WES+DHA intake was associated with more global LV thickening and with diastolic dysfunction, compared to WES feeding alone. Myocyte hypertrophy, possibly related to subcellular injury, is an early change that may contribute to gross hypertrophy. Strain differences in adipokines and myocardial fatty acid accretion may underlie heterogeneous data from rodent studies.     

Molecular evolution of mammalian lactate dehydrogenase-A genes and pseudogenes: association of a mouse processed pseudogene with a B1 repetitive sequence

A mouse genomic clone containing a lactate dehydrogenase-A (LDH-A) processed pseudogene and a B1 repetitive element was isolated, and a nucleotide sequence of approximately 3 kb was determined. The pseudogene and B1 element are flanked by perfect 13-bp repeats, and the B1 sequence starts at 14 nucleotides 3' to the presumptive polyadenylation signal of the pseudogene. The nucleotide sequences of the LDH-A genes and processed pseudogenes from mouse, rat, and human were compared, and a phylogenetic tree was constructed. The rate and pattern of nucleotide substitutions in the LDH-A pseudogenes are similar to previously reported results (Li et al. 1984). The average rate of nucleotide substitutions in the LDH-A pseudogenes is 4.3 X 10(- 9)/site/year. The substitutions of C----T and G----A are most frequent, and A----G substitutions are relatively high. The rate of synonymous substitutions in the LDH-A genes is 5.3 X 10(-9), which is not significantly higher than the average rate of 4.7 X 10(-9) for 35 mammalian genes. The rate of nonsynonymous substitutions in the LDH-A genes is 0.20 X 10(-9), which is considerably lower than the average rate of 0.88 X 10(-9) for 35 mammalian genes. Thus, the mammalian LDH-A gene appears to be highly conserved in evolution.      

The glucose transport system of the hyperthermophilic anaerobic bacterium Thermotoga neapolitana

The glucose transport system of the extremely thermophilic anaerobic bacterium Thermotoga neapolitana was studied with the nonmetabolizable glucose analog 2-deoxy-D-glucose (2-DOG). T. neapolitana accumulated 2-DOG against a concentration gradient in an intracellular free sugar pool that was exchangeable with external source of energy, such as pyruvate, and was inhibited by arsenate and gramicidin D. There was no phosphoenolpyruvate-dependent phosphorylation of glucose, 2-DOG, or fructose by cell extracts or toluene-treated cells, indicating the absence of a phosphoenolpyruvate:sugar phosphotransferase system. These data indicate that D-glucose is taken up by T. neapolitana via an active transport system that is energized by an ion gradient generated by ATP, derived from substrate-level phosphorylation.     

Ustilago maydis Mating Hyphae Orient Their Growth toward Pheromone Sources

Snetselaar, K. M., Bolker, M., and Kahmann, R. 1996. Ustilago maydis mating hyphae orient their growth toward pheromone sources. Fungal Genetics and Biology 20, 299-312. When small drops of Ustilago maydis sporidia were placed 100-200 μm apart on agar surfaces and covered with paraffin oil, sporidia from one drop formed thin hyphae that grew in a zig-zag fashion toward the other drop if it contained sporidia making the appropriate pheromone. For example, a2b2 mating hyphae grew toward a1b1 and a1b2 mating hyphae, and the filaments eventually fused tip to tip. Time-lapse photography indicated that the mating hyphae can rapidly change orientation in response to nearby compatible sporidia. When exposed to pheromone produced by cells in an adjacent drop, haploid sporidia with the a2 allele began elongating before sporidia with the a1 allele. Sporidia without functional pheromone genes responded to pheromone although they did not induce a response, and sporidia without pheromone receptors induced formation of mating hyphae although they did not form mating hyphae. Diploid sporidia heterozygous at b but not at a formed straight, rigid, aerial filaments when exposed to pheromone produced by the appropriate haploid sporidia. Again, the a2a2b1b2 strain formed filaments more quickly than the a1a1b1b2 strain. Taken together, these results suggest that the a2 pheromone diffuses less readily or is degraded more quickly than the a1 pheromone.