Assessment of Digital PCR as a Primary Reference Measurement Procedure to Support Advances in Precision Medicine.

BACKGROUND: Genetic testing of tumor tissue and circulating cell-free DNA for somatic variants guides patient treatment of many cancers. Such measurements will be fundamental in the future support of precision medicine. However, there are currently no primary reference measurement procedures available for nucleic acid quantification that would support translation of tests for circulating tumor DNA into routine use. METHODS: We assessed the accuracy of digital PCR (dPCR) for copy number quantification of a frequently occurring single-nucleotide variant in colorectal cancer (KRAS c.35G>A, p.Gly12Asp, from hereon termed G12D) by evaluating potential sources of uncertainty that influence dPCR measurement. RESULTS: Concentration values for samples of KRAS G12D and wild-type plasmid templates varied by <1.2-fold when measured using 5 different assays with varying detection chemistry (hydrolysis, scorpion probes, and intercalating dyes) and <1.3-fold with 4 commercial dPCR platforms. Measurement trueness of a selected dPCR assay and platform was validated by comparison with an orthogonal method (inductively coupled plasma mass spectrometry). The candidate dPCR reference measurement procedure showed linear quantification over a wide range of copies per reaction and high repeatability and interlaboratory reproducibility (CV, 2%-8% and 5%-10%, respectively). CONCLUSIONS: This work validates dPCR as an SI-traceable reference measurement procedure based on enumeration and demonstrates how it can be applied for assignment of copy number concentration and fractional abundance values to DNA reference materials in an aqueous solution. High-accuracy measurements using dPCR will support the implementation and traceable standardization of molecular diagnostic procedures needed for advancements in precision medicine.

International Interlaboratory Digital PCR Study Demonstrating High Reproducibility for the Measurement of a Rare Sequence Variant.

This study tested the claim that digital PCR (dPCR) can offer highly reproducible quantitative measurements in disparate laboratories. Twenty-one laboratories measured four blinded samples containing different quantities of a KRAS fragment encoding G12D, an important genetic marker for guiding therapy of certain cancers. This marker is challenging to quantify reproducibly using quantitative PCR (qPCR) or next generation sequencing (NGS) due to the presence of competing wild type sequences and the need for calibration. Using dPCR, 18 laboratories were able to quantify the G12D marker within 12% of each other in all samples. Three laboratories appeared to measure consistently outlying results; however, proper application of a follow-up analysis recommendation rectified their data. Our findings show that dPCR has demonstrable reproducibility across a large number of laboratories without calibration. This could enable the reproducible application of molecular stratification to guide therapy and, potentially, for molecular diagnostics.

Inter-laboratory assessment of different digital PCR platforms for quantification of human cytomegalovirus DNA.

Quantitative PCR (qPCR) is an important tool in pathogen detection. However, the use of different qPCR components, calibration materials and DNA extraction methods reduces comparability between laboratories, which can result in false diagnosis and discrepancies in patient care. The wider establishment of a metrological framework for nucleic acid tests could improve the degree of standardisation of pathogen detection and the quantification methods applied in the clinical context. To achieve this, accurate methods need to be developed and implemented as reference measurement procedures, and to facilitate characterisation of suitable certified reference materials. Digital PCR (dPCR) has already been used for pathogen quantification by analysing nucleic acids. Although dPCR has the potential to provide robust and accurate quantification of nucleic acids, further assessment of its actual performance characteristics is needed before it can be implemented in a metrological framework, and to allow adequate estimation of measurement uncertainties. Here, four laboratories demonstrated reproducibility (expanded measurement uncertainties below 15%) of dPCR for quantification of DNA from human cytomegalovirus, with no calibration to a common reference material. Using whole-virus material and extracted DNA, an intermediate precision (coefficients of variation below 25%) between three consecutive experiments was noted. Furthermore, discrepancies in estimated mean DNA copy number concentrations between laboratories were less than twofold, with DNA extraction as the main source of variability. These data demonstrate that dPCR offers a repeatable and reproducible method for quantification of viral DNA, and due to its satisfactory performance should be considered as candidate for reference methods for implementation in a metrological framework.

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.

Assessing the accuracy of quantitative molecular microbial profiling.

The application of high-throughput sequencing in profiling microbial communities is providing an unprecedented ability to investigate microbiomes. Such studies typically apply one of two methods: amplicon sequencing using PCR to target a conserved orthologous sequence (typically the 16S ribosomal RNA gene) or whole (meta)genome sequencing (WGS). Both methods have been used to catalog the microbial taxa present in a sample and quantify their respective abundances. However, a comparison of the inherent precision or bias of the different sequencing approaches has not been performed. We previously developed a metagenomic control material (MCM) to investigate error when performing different sequencing strategies. Amplicon sequencing using four different primer strategies and two 16S rRNA regions was examined (Roche 454 Junior) and compared to WGS (Illumina HiSeq). All sequencing methods generally performed comparably and in good agreement with organism specific digital PCR (dPCR); WGS notably demonstrated very high precision. Where discrepancies between relative abundances occurred they tended to differ by less than twofold. Our findings suggest that when alternative sequencing approaches are used for microbial molecular profiling they can perform with good reproducibility, but care should be taken when comparing small differences between distinct methods. This work provides a foundation for future work comparing relative differences between samples and the impact of extraction methods. We also highlight the value of control materials when conducting microbial profiling studies to benchmark methods and set appropriate thresholds.

A standardised methodology for the extraction and quantification of cell-free DNA in cerebrospinal fluid and application to evaluation of Alzheimer's disease and brain cancers.

Cerebrospinal fluid (CSF) is a source of diagnostic biomarkers for a range of neurological conditions. Cell-free DNA (cfDNA) is detected in CSF and differences in the concentration of cell-free mitochondrial DNA have been reported in studies of neurodegenerative disorders including Alzheimer's disease (AD). However, the influence of pre-analytical steps has not been investigated for cfDNA in CSF and there is no standardised approach for quantification of total cfDNA (copies of nuclear genome or mitochondria-derived gene targets). In this study, the suitability of four extraction methods was evaluated: QIAamp Circulating Nucleic Acid (Qiagen), Quick-cfDNA Serum & Plasma (Zymo), NucleoSnap® DNA Plasma (Macherey-Nagel) and Plasma/Serum Circulating DNA Purification Mini (Norgen) kits, for cfDNA extraction from CSF of controls and AD dementia patients, utilising a spike-in control for extraction efficiency and fragment size. One of the optimal extraction methods was applied to a comparison of cfDNA concentrations in CSF from control subjects, AD dementia and primary and secondary brain tumour patients. Extraction efficiency based on spike-in recovery was similar in all three groups whilst both endogenous mitochondrial and nucleus-derived cfDNA was significantly higher in CSF from cancer patients compared to control and AD groups, which typically contained < 100 genome copies/mL. This study shows that it is feasible to measure low concentration nuclear and mitochondrial gene targets in CSF and that normalisation of extraction yield can help control pre-analytical variability influencing biomarker measurements.

Defining the early stages of intestinal colonisation by whipworms.

Whipworms are large metazoan parasites that inhabit multi-intracellular epithelial tunnels in the large intestine of their hosts, causing chronic disease in humans and other mammals. How first-stage larvae invade host epithelia and establish infection remains unclear. Here we investigate early infection events using both Trichuris muris infections of mice and murine caecaloids, the first in-vitro system for whipworm infection and organoid model for live helminths. We show that larvae degrade mucus layers to access epithelial cells. In early syncytial tunnels, larvae are completely intracellular, woven through multiple live dividing cells. Using single-cell RNA sequencing of infected mouse caecum, we reveal that progression of infection results in cell damage and an expansion of enterocytes expressing of Isg15, potentially instigating the host immune response to the whipworm and tissue repair. Our results unravel intestinal epithelium invasion by whipworms and reveal specific host-parasite interactions that allow the whipworm to establish its multi-intracellular niche.

An inter-laboratory study to investigate the impact of the bioinformatics component on microbiome analysis using mock communities.

Despite the advent of whole genome metagenomics, targeted approaches (such as 16S rRNA gene amplicon sequencing) continue to be valuable for determining the microbial composition of samples. Amplicon microbiome sequencing can be performed on clinical samples from a normally sterile site to determine the aetiology of an infection (usually single pathogen identification) or samples from more complex niches such as human mucosa or environmental samples where multiple microorganisms need to be identified. The methodologies are frequently applied to determine both presence of micro-organisms and their quantity or relative abundance. There are a number of technical steps required to perform microbial community profiling, many of which may have appreciable precision and bias that impacts final results. In order for these methods to be applied with the greatest accuracy, comparative studies across different laboratories are warranted. In this study we explored the impact of the bioinformatic approaches taken in different laboratories on microbiome assessment using 16S rRNA gene amplicon sequencing results. Data were generated from two mock microbial community samples which were amplified using primer sets spanning five different variable regions of 16S rRNA genes. The PCR-sequencing analysis included three technical repeats of the process to determine the repeatability of their methods. Thirteen laboratories participated in the study, and each analysed the same FASTQ files using their choice of pipeline. This study captured the methods used and the resulting sequence annotation and relative abundance output from bioinformatic analyses. Results were compared to digital PCR assessment of the absolute abundance of each target representing each organism in the mock microbial community samples and also to analyses of shotgun metagenome sequence data. This ring trial demonstrates that the choice of bioinformatic analysis pipeline alone can result in different estimations of the composition of the microbiome when using 16S rRNA gene amplicon sequencing data. The study observed differences in terms of both presence and abundance of organisms and provides a resource for ensuring reproducible pipeline development and application. The observed differences were especially prevalent when using custom databases and applying high stringency operational taxonomic unit (OTU) cut-off limits. In order to apply sequencing approaches with greater accuracy, the impact of different analytical steps needs to be clearly delineated and solutions devised to harmonise microbiome analysis results.

A comparison of miRNA isolation and RT-qPCR technologies and their effects on quantification accuracy and repeatability.

MicroRNAs (miRNAs) are short (~22 nucleotides), non-coding RNA molecules that post-transcriptionally regulate gene expression. As the miRNA field is still in its relative infancy, there is currently a lack of consensus regarding optimal methodologies for miRNA quantification, data analysis and data standardization. To investigate miRNA measurement we selected a panel of both synthetic miRNA spikes and endogenous miRNAs to evaluate assay performance, copy number estimation, and relative quantification. We compared two different miRNA quantification methodologies and also assessed the impact of short RNA enrichment on the miRNA measurement. We found that both short RNA enrichment and quantification strategy used had a significant impact on miRNA measurement. Our findings illustrate that miRNA quantification can be influenced by the choice of methodology and this must be considered when interpreting miRNA analyses. Furthermore, we show that synthetic miRNA spikes can be used as effective experimental controls for the short RNA enrichment procedure.

TGF-ß/Smad2/3 signaling directly regulates several miRNAs in mouse ES cells and early embryos.

The Transforming Growth Factor-ß (TGF-ß) signaling pathway is one of the major pathways essential for normal embryonic development and tissue homeostasis, with anti-tumor but also pro-metastatic properties in cancer. This pathway directly regulates several target genes that mediate its downstream functions, however very few microRNAs (miRNAs) have been identified as targets. miRNAs are modulators of gene expression with essential roles in development and a clear association with diseases including cancer. Little is known about the transcriptional regulation of the primary transcripts (pri-miRNA, pri-miR) from which several mature miRNAs are often derived. Here we present the identification of miRNAs regulated by TGF-ß signaling in mouse embryonic stem (ES) cells and early embryos. We used an inducible ES cell system to maintain high levels of the TGF-ß activated/phosphorylated Smad2/3 effectors, which are the transcription factors of the pathway, and a specific inhibitor that blocks their activation. By performing short RNA deep-sequencing after 12 hours Smad2/3 activation and after 16 hours inhibition, we generated a database of responsive miRNAs. Promoter/enhancer analysis of a subset of these miRNAs revealed that the transcription of pri-miR-181c/d and the pri-miR-341∼3072 cluster were found to depend on activated Smad2/3. Several of these miRNAs are expressed in early mouse embryos, when the pathway is known to play an essential role. Treatment of embryos with TGF-ß inhibitor caused a reduction of their levels confirming that they are targets of this pathway in vivo. Furthermore, we showed that pri-miR-341∼3072 transcription also depends on FoxH1, a known Smad2/3 transcription partner during early development. Together, our data show that miRNAs are regulated directly by the TGF-ß/Smad2/3 pathway in ES cells and early embryos. As somatic abnormalities in functions known to be regulated by the TGF-ß/Smad2/3 pathway underlie tumor suppression and metastasis, this research also provides a resource for miRNAs involved in cancer.

microRNA-449 is a putative regulator of choroid plexus development and function.

microRNAs are short RNA molecules that are often expressed in specific tissues and regulate a variety of developmental processes. We used locked nucleic acid probes in in situ hybridisation reactions to study the distribution of microRNA-449 (miR449) during mouse embryonic development in order to obtain clues about its function/s. Between E9.75 and E11.5, miR449 was found to be expressed specifically in the developing roof plate of the fourth ventricle within the domain of roof plate marker, Lmx1a. From E12.5 onwards, this expression became restricted to the epithelial cell layer of the fourth ventricle choroid plexus. MiR449 also became detectable specifically in the choroid plexuses of the lateral and 3rd ventricles at E13.5 and E15.5, respectively. Northern blot analysis of adult brain also showed a selective and enriched expression in the choroid plexus tissue. Potential target genes regulated by miR449 were selected for experimental validation in luciferase-reporter assays and the transcription factor E2f5, which regulates CSF production, was verified as a miR449 target gene. Taken together, these findings suggest that miR449 has a specific role in the development and functioning of choroid plexuses.