Benchmarking digital PCR partition classification methods with empirical and simulated duplex data.

Digital PCR (dPCR) is a highly accurate technique for the quantification of target nucleic acid(s). It has shown great potential in clinical applications, like tumor liquid biopsy and validation of biomarkers. Accurate classification of partitions based on end-point fluorescence intensities is crucial to avoid biased estimators of the concentration of the target molecules. We have evaluated many clustering methods, from general-purpose methods to specific methods for dPCR and flowcytometry, on both simulated and real-life data. Clustering method performance was evaluated by simulating various scenarios. Based on our extensive comparison of clustering methods, we describe the limits of these methods, and formulate guidelines for choosing an appropriate method. In addition, we have developed a novel method for simulating realistic dPCR data. The method is based on a mixture distribution of a Poisson point process and a skew-$t$ distribution, which enables the generation of irregularities of cluster shapes and randomness of partitions between clusters ('rain') as commonly observed in dPCR data. Users can fine-tune the model parameters and generate labeled datasets, using their own data as a template. Besides, the database of experimental dPCR data augmented with the labeled simulated data can serve as training and testing data for new clustering methods. The simulation method is available as an R Shiny app.

Digital PCR Partition Classification.

BACKGROUND: Partition classification is a critical step in the digital PCR data analysis pipeline. A range of partition classification methods have been developed, many motivated by specific experimental setups. An overview of these partition classification methods is lacking and their comparative properties are often unclear, likely impacting the proper application of these methods. CONTENT: This review provides a summary of all available digital PCR partition classification approaches and the challenges they aim to overcome, serving as a guide for the digital PCR practitioner wishing to apply them. We additionally discuss strengths and weaknesses of these methods, which can further guide practitioners in vigilant application of these existing methods. This review provides method developers with ideas for improving methods or designing new ones. The latter is further stimulated by our identification and discussion of application gaps in the literature, for which there are currently no or few methods available. SUMMARY: This review provides an overview of digital PCR partition classification methods, their properties, and potential applications. Ideas for further advances are presented and may bolster method development.

Identification of single target taxon-specific reference assays for the most commonly genetically transformed crops using digital droplet PCR.

Knowledge of the number of DNA sequences targeted by the taxon-specific reference assays is essential for correct GM quantification and is key to the harmonisation of measurement results. In the present study droplet digital PCR (ddPCR) was used to determine the number of DNA target copies of taxon-specific assays validated for real-time PCR for the four main genetically modified (GM) crops. The transferability of experimental conditions from real-time PCR to ddPCR was also explored, as well as the effect of DNA digestion. The results of this study indicate that for each crop at least one taxon-specific assay can be identified as having a single DNA target. A short list of taxon-specific reference assays is proposed as best candidates for the relative quantification of GM events for soybean, maize, cotton and oilseed rape. The investigated assays could be in most cases transferred to ddPCR without further optimisation. The use of DNA digestion did not improve ddPCR characteristics such as rain and resolution at the conditions tested.

Novel nuclear barcode regions for the identification of flatfish species.

The development of an efficient seafood traceability framework is crucial for the management of sustainable fisheries and the monitoring of potential substitution fraud across the food chain. Recent studies have shown the potential of DNA barcoding methods in this framework, with most of the efforts focusing on using mitochondrial targets such as the cytochrome oxidase 1 and cytochrome b genes. In this article, we show the identification of novel targets in the nuclear genome, and their associated primers, to be used for the efficient identification of flatfishes of the Pleuronectidae family. In addition, different in silico methods are described to generate a dataset of barcode reference sequences from the ever-growing wealth of publicly available sequence information, replacing, where possible, labour-intensive laboratory work. The short amplicon lengths render the analysis of these new barcode target regions ideally suited to next-generation sequencing techniques, allowing characterisation of multiple fish species in mixed and processed samples. Their location in the nucleus also improves currently used methods by allowing the identification of hybrid individuals.

Evaluation of qPCR curve analysis methods for reliable biomarker discovery: bias, resolution, precision, and implications.

RNA transcripts such as mRNA or microRNA are frequently used as biomarkers to determine disease state or response to therapy. Reverse transcription (RT) in combination with quantitative PCR (qPCR) has become the method of choice to quantify small amounts of such RNA molecules. In parallel with the democratization of RT-qPCR and its increasing use in biomedical research or biomarker discovery, we witnessed a growth in the number of gene expression data analysis methods. Most of these methods are based on the principle that the position of the amplification curve with respect to the cycle-axis is a measure for the initial target quantity: the later the curve, the lower the target quantity. However, most methods differ in the mathematical algorithms used to determine this position, as well as in the way the efficiency of the PCR reaction (the fold increase of product per cycle) is determined and applied in the calculations. Moreover, there is dispute about whether the PCR efficiency is constant or continuously decreasing. Together this has lead to the development of different methods to analyze amplification curves. In published comparisons of these methods, available algorithms were typically applied in a restricted or outdated way, which does not do them justice. Therefore, we aimed at development of a framework for robust and unbiased assessment of curve analysis performance whereby various publicly available curve analysis methods were thoroughly compared using a previously published large clinical data set (Vermeulen et al., 2009) [11]. The original developers of these methods applied their algorithms and are co-author on this study. We assessed the curve analysis methods' impact on transcriptional biomarker identification in terms of expression level, statistical significance, and patient-classification accuracy. The concentration series per gene, together with data sets from unpublished technical performance experiments, were analyzed in order to assess the algorithms' precision, bias, and resolution. While large differences exist between methods when considering the technical performance experiments, most methods perform relatively well on the biomarker data. The data and the analysis results per method are made available to serve as benchmark for further development and evaluation of qPCR curve analysis methods (http://qPCRDataMethods.hfrc.nl).

Enhanced analysis of real-time PCR data by using a variable efficiency model: FPK-PCR.

Current methodology in real-time Polymerase chain reaction (PCR) analysis performs well provided PCR efficiency remains constant over reactions. Yet, small changes in efficiency can lead to large quantification errors. Particularly in biological samples, the possible presence of inhibitors forms a challenge. We present a new approach to single reaction efficiency calculation, called Full Process Kinetics-PCR (FPK-PCR). It combines a kinetically more realistic model with flexible adaptation to the full range of data. By reconstructing the entire chain of cycle efficiencies, rather than restricting the focus on a 'window of application', one extracts additional information and loses a level of arbitrariness. The maximal efficiency estimates returned by the model are comparable in accuracy and precision to both the golden standard of serial dilution and other single reaction efficiency methods. The cycle-to-cycle changes in efficiency, as described by the FPK-PCR procedure, stay considerably closer to the data than those from other S-shaped models. The assessment of individual cycle efficiencies returns more information than other single efficiency methods. It allows in-depth interpretation of real-time PCR data and reconstruction of the fluorescence data, providing quality control. Finally, by implementing a global efficiency model, reproducibility is improved as the selection of a window of application is avoided.

Detection and Quantification of Botanical Impurities in Commercial Oregano (Origanum vulgare) Using Metabarcoding and Digital PCR.

DNA technology for food authentication is already well established, and with the advent of Next Generation Sequencing (NGS) and, more specifically, metabarcoding, compositional analysis of food at the molecular level has rapidly gained popularity. This has led to several reports in the media about the presence of foreign, non-declared species in several food commodities. As herbs and spices are attractive targets for fraudulent manipulation, a combination of digital PCR and metabarcoding by NGS was employed to check the purity of 285 oregano samples taken from the European market. By using novel primers and analytical approaches, it was possible to detect and quantify both adulterants and contaminants in these samples. The results highlight the high potential of NGS for compositional analysis, although its quantitative information (read count percentages) is unreliable, and other techniques are therefore needed to complement the sequencing information for assessing authenticity ('true to the name') of food ingredients.

DNA Accounting: Tallying Genomes to Detect Adulterated Saffron.

The EU General Food Law not only aims at ensuring food safety but also to 'prevent fraudulent or deceptive practices; the adulteration of food; and any other practices which may mislead the consumer'. Especially the partial or complete, deliberate, and intentional substitution of valuable ingredients (e.g., Saffron) for less valuable ones is of concern. Due to the variety of products on the market an approach to detect food adulteration that works well for one species may not be easily applicable to another. Here we present a broadly applicable approach for the detection of substitution of biological materials based on digital PCR. By simultaneously measuring and forecasting the number of genome copies in a sample, fraud is detectable as a discrepancy between these two values. Apart from the choice of target gene, the procedure is identical across all species. It is scalable, rapid, and has a high dynamic range. We provide proof of concept by presenting the analysis of 141 samples of Saffron (Crocus sativus) from across the European market by DNA accounting and the verification of these results by NGS analysis.

A theoretical introduction to "combinatory SYBRGreen qPCR screening", a matrix-based approach for the detection of materials derived from genetically modified plants.

The detection of genetically modified (GM) materials in food and feed products is a complex multi-step analytical process invoking screening, identification, and often quantification of the genetically modified organisms (GMO) present in a sample. "Combinatory qPCR SYBRGreen screening" (CoSYPS) is a matrix-based approach for determining the presence of GM plant materials in products. The CoSYPS decision-support system (DSS) interprets the analytical results of SYBRGREEN qPCR analysis based on four values: the C(t)- and T(m) values and the LOD and LOQ for each method. A theoretical explanation of the different concepts applied in CoSYPS analysis is given (GMO Universe, "Prime number tracing", matrix/combinatory approach) and documented using the RoundUp Ready soy GTS40-3-2 as an example. By applying a limited set of SYBRGREEN qPCR methods and through application of a newly developed "prime number"-based algorithm, the nature of subsets of corresponding GMO in a sample can be determined. Together, these analyses provide guidance for semi-quantitative estimation of GMO presence in a food and feed product.

The EU one-stop-shop collection of publicly available information on COVID-19 in vitro diagnostic medical devices.

The JRC COVID-19 In Vitro Diagnostic Devices and Test Methods Database, aimed to collect in a single place all publicly available information on performance of CE-marked in vitro diagnostic medical devices (IVDs) as well as in house laboratory-developed devices and related test methods for COVID-19, is here presented. The database, manually curated and regularly updated, has been developed as a follow-up to the Communication from the European Commission "Guidelines on in vitro diagnostic tests and their performance" of 15 April 2020 and is freely accessible at https://covid-19-diagnostics.jrc.ec.europa.eu/.

Nuclear DNA barcodes for cod identification in mildly-treated and processed food products.

Gadoids are a group of fish with historical importance in the fishing industry. The high demand for cod is one of the reasons why cod products are often mislabelled, and numerous observations have been made on the replacement of Atlantic cod (Gadus morhua) by cheaper species or its illegal capture in contravention of fish quotas. Fish species identification is traditionally based on morphological features, but this may be difficult in case of heat-treated or processed products, or where the species look similar, as in the Gadoid group. DNA-based approaches (using either nuclear or mitochondrial DNA) are most commonly used in this case, due to their high specificity and to the high resilience of the target molecules to food processing techniques. In this article, we identified, using an automated screening approach, novel barcode regions and their associated primers in the nuclear genome, to be used for the efficient identification of Gadoids. The barcode regions were tested on official and commercial samples, raw or mildly treated products, like frozen, or salted, as well as pre-cooked complex mixtures and processed samples, using next-generation sequencing (NGS) technique. The method proposed could complement existing fish identification strategies in establishing an efficient framework to detect and prevent frauds along the food chain.

Simulation of between repeat variability in real time PCR reactions.

While many decisions rely on real time quantitative PCR (qPCR) analysis few attempts have hitherto been made to quantify bounds of precision accounting for the various sources of variation involved in the measurement process. Besides influences of more obvious factors such as camera noise and pipetting variation, changing efficiencies within and between reactions affect PCR results to a degree which is not fully recognized. Here, we develop a statistical framework that models measurement error and other sources of variation as they contribute to fluorescence observations during the amplification process and to derived parameter estimates. Evaluation of reproducibility is then based on simulations capable of generating realistic variation patterns. To this end, we start from a relatively simple statistical model for the evolution of efficiency in a single PCR reaction and introduce additional error components, one at a time, to arrive at stochastic data generation capable of simulating the variation patterns witnessed in repeated reactions (technical repeats). Most of the variation in C(q) values was adequately captured by the statistical model in terms of foreseen components. To recreate the dispersion of the repeats' plateau levels while keeping the other aspects of the PCR curves within realistic bounds, additional sources of reagent consumption (side reactions) enter into the model. Once an adequate data generating model is available, simulations can serve to evaluate various aspects of PCR under the assumptions of the model and beyond.