Circulating tumor DNA is a prognostic marker of tumor recurrence in stage II and III colorectal cancer: multicentric, prospective cohort study (ALGECOLS).

BACKGROUND: In non-metastatic colorectal cancer (CRC), we evaluated prospectively the pertinence of longitudinal detection and quantification of circulating tumor DNA (ctDNA) as a prognostic marker of recurrence. METHOD: The presence of ctDNA was assessed from plasma collected before and after surgery for 184 patients classified as stage II or III and at each visit during 3-4 years of follow-up. The ctDNA analysis was performed by droplet-based digital polymerase chain reaction, targeting mutation and methylation markers, blindly from the clinical outcomes. Multivariate analyses were adjusted on age, gender, stage, and adjuvant chemotherapy. RESULTS: Before surgery, 27.5% of patients were positive for ctDNA detection. The rate of recurrence was 32.7% and 11.6% in patients with or without detectable ctDNA respectively (P = 0.001). Time to recurrence (TTR) was significantly shorter in patients with detectable ctDNA before (adjusted hazard ratio [HR] = 3.58, 95% confidence interval [CI] 1.71-7.47) or immediately after surgery (adjusted HR = 3.22, 95% CI 1.32-7.89). The TTR was significantly shorter in patients with detectable ctDNA during the early postoperative follow-up (1-6 months) (adjusted HR = 5, 95% CI 1.9-12.9). Beyond this period, ctDNA remained a prognostic marker with a median anticipated diagnosis of recurrence of 13.1 weeks (interquartile range 28 weeks) when compared to imaging follow-up. The rate of ctDNA+ might be underestimated knowing that consensus pre-analytical conditions were not described at initiation of the study. CONCLUSION: This prospective study confirms the relevance of ctDNA as a recurrence risk factor in stage II and III CRC before surgery and as a marker of minimal residual disease after surgery that may predict recurrence several months before imaging techniques.

A Study of Hypermethylated Circulating Tumor DNA as a Universal Colorectal Cancer Biomarker.

BACKGROUND: Circulating tumor DNA (ctDNA) has emerged as a good candidate for tracking tumor dynamics in different cancer types, potentially avoiding repeated tumor biopsies. Many different genes can be mutated within a tumor, complicating procedures for tumor monitoring, even with highly sensitive next-generation sequencing (NGS) strategies. Droplet-based digital PCR (dPCR) is a highly sensitive and quantitative procedure, allowing detection of very low amounts of circulating tumor genetic material, but can be limited in the total number of target loci monitored. METHODS: We analyzed hypermethylation of 3 genes, by use of droplet-based dPCR in different stages of colorectal cancer (CRC), to identify universal markers for tumor follow-up. RESULTS: Hypermethylation of WIF1 (WNT inhibitory factor 1) and NPY (neuropeptide Y) genes was significantly higher in tumor tissue compared to normal tissue, independently of tumor stage. All tumor tissues appeared positive for one of the 2 markers. Methylated ctDNA (MetctDNA) was detected in 80% of metastatic CRC and 45% of localized CRC. For samples with detectable mutations in ctDNA, MetctDNA and mutant ctDNA (MutctDNA) fractions were correlated. During follow-up of different stage CRC patients, MetctDNA changes allowed monitoring of tumor evolution. CONCLUSIONS: These results indicate that MetctDNA could be used as a universal surrogate marker for tumor follow-up in CRC patients, and monitoring MetctDNA by droplet-based dPCR could avoid the need for monitoring mutations.

Clinical relevance of KRAS-mutated subclones detected with picodroplet digital PCR in advanced colorectal cancer treated with anti-EGFR therapy.

PURPOSE: KRAS mutations are predictive of nonresponse to anti-EGFR therapies in metastatic colorectal cancer (mCRC). However, only 50% of nonmutated patients benefit from them. KRAS-mutated subclonal populations nondetectable by conventional methods have been suggested as the cause of early progression. Molecular analysis technology with high sensitivity and precision is required to test this hypothesis. EXPERIMENTAL DESIGN: From two cohorts of patients with mCRC, 136 KRAS, NRAS, and BRAF wild-type tumors with sufficient tumor material to perform highly sensitive picodroplet digital PCR (dPCR) and 41 KRAS-mutated tumors were selected. All these patients were treated by anti-EGFR therapy. dPCR was used for KRAS or BRAF mutation screening and compared with qPCR. Progression-free survival (PFS) and overall survival (OS) were analyzed according to the KRAS-mutated allele fraction. RESULTS: In addition to the confirmation of the 41 patients with KRAS-mutated tumors, dPCR also identified KRAS mutations in 22 samples considered as KRAS wild-type by qPCR. The fraction of KRAS-mutated allele quantified by dPCR was inversely correlated with anti-EGFR therapy response rate (P < 0.001). In a Cox model, the fraction of KRAS-mutated allele was associated with worse PFS and OS. Patients with less than 1% of mutant KRAS allele have similar PFS and OS than those with wild-type KRAS tumors. CONCLUSIONS: This study suggests that patients with mCRC with KRAS-mutated subclones (at least those with a KRAS-mutated subclones fraction lower or equal to 1%) had a benefit from anti-EGFR therapies.

Determining lower limits of detection of digital PCR assays for cancer-related gene mutations.

Digital PCR offers very high sensitivity compared to many other technologies for processing molecular detection assays. Herein, a process is outlined for determining the lower limit of detection (LoD) of two droplet-based digital PCR assays for point mutations of the epidermal growth factor receptor (EGFR) gene. Hydrolysis probe mutation-detection assays for EGFR p.L858R and p.T790M mutations were characterized in detail. Furthermore, sixteen additional cancer-related mutation assays were explored by the same approach. For the EGFR L8585R assay, the assay sensitivity is extremely good, and thus, the LoD is limited by the amount of amplifiable DNA that is analyzed. With 95% confidence limits, the LoD is one mutant in 180,000 wild-type molecules for the evaluation of 3.3 µg of genomic DNA, and detection of one mutant molecule in over 4 million wild-type molecules was achieved when 70 million copies of DNA were processed. The measured false-positive rate for the EGFR L8585R assay is one in 14 million, which indicates the theoretical LoD if an unlimited amount of DNA is evaluated. For the EFGR T790M assay, the LoD is one mutant in 13,000 for analysis of a 3.3 µg sample of genomic DNA, and the dPCR assay limit sensitivity approaches one mutant in 22,000 wild-type molecules.

Multiplex picodroplet digital PCR to detect KRAS mutations in circulating DNA from the plasma of colorectal cancer patients.

BACKGROUND: Multiplex digital PCR (dPCR) enables noninvasive and sensitive detection of circulating tumor DNA with performance unachievable by current molecular-detection approaches. Furthermore, picodroplet dPCR facilitates simultaneous screening for multiple mutations from the same sample. METHODS: We investigated the utility of multiplex dPCR to screen for the 7 most common mutations in codons 12 and 13 of the KRAS (Kirsten rat sarcoma viral oncogene homolog) oncogene from plasma samples of patients with metastatic colorectal cancer. Fifty plasma samples were tested from patients for whom the primary tumor biopsy tissue DNA had been characterized by quantitative PCR. RESULTS: Tumor characterization revealed that 19 patient tumors had KRAS mutations. Multiplex dPCR analysis of the plasma DNA prepared from these samples identified 14 samples that matched the mutation identified in the tumor, 1 sample contained a different KRAS mutation, and 4 samples had no detectable mutation. Among the tumor samples that were wild type for KRAS, 2 KRAS mutations were identified in the corresponding plasma samples. Duplex dPCR (i.e., wild-type and single-mutation assay) was also used to analyze plasma samples from patients with KRAS-mutated tumors and 5 samples expected to contain the BRAF (v-raf murine sarcoma viral oncogene homolog B) V600E mutation. The results for the duplex analysis matched those for the multiplex analysis for KRAS-mutated samples and, owing to its higher sensitivity, enabled detection of 2 additional samples with low levels of KRAS-mutated DNA. All 5 samples with BRAF mutations were detected. CONCLUSIONS: This work demonstrates the clinical utility of multiplex dPCR to screen for multiple mutations simultaneously with a sensitivity sufficient to detect mutations in circulating DNA obtained by noninvasive blood collection.

Multiplex picoliter-droplet digital PCR for quantitative assessment of DNA integrity in clinical samples.

BACKGROUND: Assessment of DNA integrity and quantity remains a bottleneck for high-throughput molecular genotyping technologies, including next-generation sequencing. In particular, DNA extracted from paraffin-embedded tissues, a major potential source of tumor DNA, varies widely in quality, leading to unpredictable sequencing data. We describe a picoliter droplet-based digital PCR method that enables simultaneous detection of DNA integrity and the quantity of amplifiable DNA. METHODS: Using a multiplex assay, we detected 4 different target lengths (78, 159, 197, and 550 bp). Assays were validated with human genomic DNA fragmented to sizes of 170 bp to 3000 bp. The technique was validated with DNA quantities as low as 1 ng. We evaluated 12 DNA samples extracted from paraffin-embedded lung adenocarcinoma tissues. RESULTS: One sample contained no amplifiable DNA. The fractions of amplifiable DNA for the 11 other samples were between 0.05% and 10.1% for 78-bp fragments and ≤1% for longer fragments. Four samples were chosen for enrichment and next-generation sequencing. The quality of the sequencing data was in agreement with the results of the DNA-integrity test. Specifically, DNA with low integrity yielded sequencing results with lower levels of coverage and uniformity and had higher levels of false-positive variants. CONCLUSIONS: The development of DNA-quality assays will enable researchers to downselect samples or process more DNA to achieve reliable genome sequencing with the highest possible efficiency of cost and effort, as well as minimize the waste of precious samples.

Accurate variant detection across non-amplified and whole genome amplified DNA using targeted next generation sequencing.

BACKGROUND: Many hypothesis-driven genetic studies require the ability to comprehensively and efficiently target specific regions of the genome to detect sequence variations. Often, sample availability is limited requiring the use of whole genome amplification (WGA). We evaluated a high-throughput microdroplet-based PCR approach in combination with next generation sequencing (NGS) to target 384 discrete exons from 373 genes involved in cancer. In our evaluation, we compared the performance of six non-amplified gDNA samples from two HapMap family trios. Three of these samples were also preamplified by WGA and evaluated. We tested sample pooling or multiplexing strategies at different stages of the tested targeted NGS (T-NGS) workflow. RESULTS: The results demonstrated comparable sequence performance between non-amplified and preamplified samples and between different indexing strategies [sequence specificity of 66.0% ± 3.4%, uniformity (coverage at 0.2× of the mean) of 85.6% ± 0.6%]. The average genotype concordance maintained across all the samples was 99.5% ± 0.4%, regardless of sample type or pooling strategy. We did not detect any errors in the Mendelian patterns of inheritance of genotypes between the parents and offspring within each trio. We also demonstrated the ability to detect minor allele frequencies within the pooled samples that conform to predicted models. CONCLUSION: Our described PCR-based sample multiplex approach and the ability to use WGA material for NGS may enable researchers to perform deep resequencing studies and explore variants at very low frequencies and cost.

Dynamic regulation of HIV-1 mRNA populations analyzed by single-molecule enrichment and long-read sequencing.

Alternative RNA splicing greatly expands the repertoire of proteins encoded by genomes. Next-generation sequencing (NGS) is attractive for studying alternative splicing because of the efficiency and low cost per base, but short reads typical of NGS only report mRNA fragments containing one or few splice junctions. Here, we used single-molecule amplification and long-read sequencing to study the HIV-1 provirus, which is only 9700 bp in length, but encodes nine major proteins via alternative splicing. Our data showed that the clinical isolate HIV-1(89.6) produces at least 109 different spliced RNAs, including a previously unappreciated ∼1 kb class of messages, two of which encode new proteins. HIV-1 message populations differed between cell types, longitudinally during infection, and among T cells from different human donors. These findings open a new window on a little studied aspect of HIV-1 replication, suggest therapeutic opportunities and provide advanced tools for the study of alternative splicing.

High-resolution dose-response screening using droplet-based microfluidics.

A critical early step in drug discovery is the screening of a chemical library. Typically, promising compounds are identified in a primary screen and then more fully characterized in a dose-response analysis with 7-10 data points per compound. Here, we describe a robust microfluidic approach that increases the number of data points to approximately 10,000 per compound. The system exploits Taylor-Aris dispersion to create concentration gradients, which are then segmented into picoliter microreactors by droplet-based microfluidics. The large number of data points results in IC(50) values that are highly precise (± 2.40% at 95% confidence) and highly reproducible (CV = 2.45%, n = 16). In addition, the high resolution of the data reveals complex dose-response relationships unambiguously. We used this system to screen a chemical library of 704 compounds against protein tyrosine phosphatase 1B, a diabetes, obesity, and cancer target. We identified a number of novel inhibitors, the most potent being sodium cefsulodine, which has an IC(50) of 27 ± 0.83 µM.

Detection of low prevalence somatic mutations in solid tumors with ultra-deep targeted sequencing.

Ultra-deep targeted sequencing (UDT-Seq) can identify subclonal somatic mutations in tumor samples. Early assays' limited breadth and depth restrict their clinical utility. Here, we target 71 kb of mutational hotspots in 42 cancer genes. We present novel methods enhancing both laboratory workflow and mutation detection. We evaluate UDT-Seq true sensitivity and specificity (> 94% and > 99%, respectively) for low prevalence mutations in a mixing experiment and demonstrate its utility using six tumor samples. With an improved performance when run on the Illumina Miseq, the UDT-Seq assay is well suited for clinical applications to guide therapy and study clonal selection in heterogeneous samples.

Pipeline for large-scale microdroplet bisulfite PCR-based sequencing allows the tracking of hepitype evolution in tumors.

Cytosine methylation provides an epigenetic level of cellular plasticity that is important for development, differentiation and cancerogenesis. We adopted microdroplet PCR to bisulfite treated target DNA in combination with second generation sequencing to simultaneously assess DNA sequence and methylation. We show measurement of methylation status in a wide range of target sequences (total 34 kb) with an average coverage of 95% (median 100%) and good correlation to the opposite strand (rho = 0.96) and to pyrosequencing (rho = 0.87). Data from lymphoma and colorectal cancer samples for SNRPN (imprinted gene), FGF6 (demethylated in the cancer samples) and HS3ST2 (methylated in the cancer samples) serve as a proof of principle showing the integration of SNP data and phased DNA-methylation information into "hepitypes" and thus the analysis of DNA methylation phylogeny in the somatic evolution of cancer.

Application of microdroplet PCR for large-scale targeted bisulfite sequencing.

Cytosine methylation of DNA CpG dinucleotides in gene promoters is an epigenetic modification that regulates gene transcription. While many methods exist to interrogate methylation states, few current methods offer large-scale, targeted, single CpG resolution. We report an approach combining bisulfite treatment followed by microdroplet PCR with next-generation sequencing to assay the methylation state of 50 genes in the regions 1 kb upstream of and downstream from their transcription start sites. This method yielded 96% coverage of the targeted CpGs and demonstrated high correlation between CpG island (CGI) DNA methylation and transcriptional regulation. The method was scaled to interrogate the methylation status of 77,674 CpGs in the promoter regions of 2100 genes in primary CD4 T cells. The 2100 gene library yielded 97% coverage of all targeted CpGs and 99% of the target amplicons.

Quantitative and sensitive detection of rare mutations using droplet-based microfluidics.

Somatic mutations within tumoral DNA can be used as highly specific biomarkers to distinguish cancer cells from their normal counterparts. These DNA biomarkers are potentially useful for the diagnosis, prognosis, treatment and follow-up of patients. In order to have the required sensitivity and specificity to detect rare tumoral DNA in stool, blood, lymph and other patient samples, a simple, sensitive and quantitative procedure to measure the ratio of mutant to wild-type genes is required. However, techniques such as dual probe TaqMan(®) assays and pyrosequencing, while quantitative, cannot detect less than ∼1% mutant genes in a background of non-mutated DNA from normal cells. Here we describe a procedure allowing the highly sensitive detection of mutated DNA in a quantitative manner within complex mixtures of DNA. The method is based on using a droplet-based microfluidic system to perform digital PCR in millions of picolitre droplets. Genomic DNA (gDNA) is compartmentalized in droplets at a concentration of less than one genome equivalent per droplet together with two TaqMan(®) probes, one specific for the mutant and the other for the wild-type DNA, which generate green and red fluorescent signals, respectively. After thermocycling, the ratio of mutant to wild-type genes is determined by counting the ratio of green to red droplets. We demonstrate the accurate and sensitive quantification of mutated KRAS oncogene in gDNA. The technique enabled the determination of mutant allelic specific imbalance (MASI) in several cancer cell-lines and the precise quantification of a mutated KRAS gene in the presence of a 200,000-fold excess of unmutated KRAS genes. The sensitivity is only limited by the number of droplets analyzed. Furthermore, by one-to-one fusion of drops containing gDNA with any one of seven different types of droplets, each containing a TaqMan(®) probe specific for a different KRAS mutation, or wild-type KRAS, and an optical code, it was possible to screen the six common mutations in KRAS codon 12 in parallel in a single experiment.

Multiplex digital PCR: breaking the one target per color barrier of quantitative PCR.

Quantitative polymerase chain reactions (qPCR) based on real-time PCR constitute a powerful and sensitive method for the analysis of nucleic acids. However, in qPCR, the ability to multiplex targets using differently colored fluorescent probes is typically limited to 4-fold by the spectral overlap of the fluorophores. Furthermore, multiplexing qPCR assays requires expensive instrumentation and most often lengthy assay development cycles. Digital PCR (dPCR), which is based on the amplification of single target DNA molecules in many separate reactions, is an attractive alternative to qPCR. Here we report a novel and easy method for multiplexing dPCR in picolitre droplets within emulsions-generated and read out in microfluidic devices-that takes advantage of both the very high numbers of reactions possible within emulsions (>10(6)) as well as the high likelihood that the amplification of only a single target DNA molecule will initiate within each droplet. By varying the concentration of different fluorogenic probes of the same color, it is possible to identify the different probes on the basis of fluorescence intensity. Adding multiple colors increases the number of possible reactions geometrically, rather than linearly as with qPCR. Accurate and precise copy numbers of up to sixteen per cell were measured using a model system. A 5-plex assay for spinal muscular atrophy was demonstrated with just two fluorophores to simultaneously measure the copy number of two genes (SMN1 and SMN2) and to genotype a single nucleotide polymorphism (c.815A>G, SMN1). Results of a pilot study with SMA patients are presented.

Microdroplet-based PCR enrichment for large-scale targeted sequencing.

Targeted enrichment of specific loci of the human genome is a promising approach to enable sequencing-based studies of genetic variation in large populations. Here we describe an enrichment approach based on microdroplet PCR, which enables 1.5 million amplifications in parallel. We sequenced six samples enriched by microdroplet or traditional singleplex PCR using primers targeting 435 exons of 47 genes. Both methods generated similarly high-quality data: 84% of the uniquely mapping reads fell within the targeted sequences; coverage was uniform across approximately 90% of targeted bases; sequence variants were called with >99% accuracy; and reproducibility between samples was high (r(2) = 0.9). We scaled the microdroplet PCR to 3,976 amplicons totaling 1.49 Mb of sequence, sequenced the resulting sample with both Illumina GAII and Roche 454, and obtained data with equally high specificity and sensitivity. Our results demonstrate that microdroplet technology is well suited for processing DNA for massively parallel enrichment of specific subsets of the human genome for targeted sequencing.

Droplet microfluidic technology for single-cell high-throughput screening.

We present a droplet-based microfluidic technology that enables high-throughput screening of single mammalian cells. This integrated platform allows for the encapsulation of single cells and reagents in independent aqueous microdroplets (1 pL to 10 nL volumes) dispersed in an immiscible carrier oil and enables the digital manipulation of these reactors at a very high-throughput. Here, we validate a full droplet screening workflow by conducting a droplet-based cytotoxicity screen. To perform this screen, we first developed a droplet viability assay that permits the quantitative scoring of cell viability and growth within intact droplets. Next, we demonstrated the high viability of encapsulated human monocytic U937 cells over a period of 4 days. Finally, we developed an optically-coded droplet library enabling the identification of the droplets composition during the assay read-out. Using the integrated droplet technology, we screened a drug library for its cytotoxic effect against U937 cells. Taken together our droplet microfluidic platform is modular, robust, uses no moving parts, and has a wide range of potential applications including high-throughput single-cell analyses, combinatorial screening, and facilitating small sample analyses.

Fluorescence-activated droplet sorting (FADS): efficient microfluidic cell sorting based on enzymatic activity.

We describe a highly efficient microfluidic fluorescence-activated droplet sorter (FADS) combining many of the advantages of microtitre-plate screening and traditional fluorescence-activated cell sorting (FACS). Single cells are compartmentalized in emulsion droplets, which can be sorted using dielectrophoresis in a fluorescence-activated manner (as in FACS) at rates up to 2000 droplets s(-1). To validate the system, mixtures of E. coli cells, expressing either the reporter enzyme beta-galactosidase or an inactive variant, were compartmentalized with a fluorogenic substrate and sorted at rates of approximately 300 droplets s(-1). The false positive error rate of the sorter at this throughput was <1 in 10(4) droplets. Analysis of the sorted cells revealed that the primary limit to enrichment was the co-encapsulation of E. coli cells, not sorting errors: a theoretical model based on the Poisson distribution accurately predicted the observed enrichment values using the starting cell density (cells per droplet) and the ratio of active to inactive cells. When the cells were encapsulated at low density ( approximately 1 cell for every 50 droplets), sorting was very efficient and all of the recovered cells were the active strain. In addition, single active droplets were sorted and cells were successfully recovered.

Detection and analysis of low-abundance cell-surface biomarkers using enzymatic amplification in microfluidic droplets.

Finding the few: Cell-surface proteins are useful disease biomarkers, but current high-throughput methods are limited to detecting cells expressing more than several hundred proteins. Enzymatic amplification in microfluidic droplets (see picture) is a high-throughput method for detection and analysis of cell-surface biomarkers expressed at very low levels on individual human cells. Droplet optical labels allow concurrent analysis of several samples.

High-throughput quantitative polymerase chain reaction in picoliter droplets.

Limiting dilution PCR has become an increasingly useful technique for the detection and quantification of rare species in a population, but the limit of detection and accuracy of quantification are largely determined by the number of reactions that can be analyzed. Increased throughput may be achieved by reducing the reaction volume and increasing processivity. We have designed a high-throughput microfluidic chip that encapsulates PCR reagents in millions of picoliter droplets in a continuous oil flow. The oil stream conducts the droplets through alternating denaturation and annealing zones, resulting in rapid (55-s cycles) and efficient PCR amplification. Inclusion of fluorescent probes in the PCR reaction mix permits the amplification process to be monitored within individual droplets at specific locations within the microfluidic chip. We show that amplification of a 245-bp adenovirus product can be detected and quantified in 35 min at starting template concentrations as low as 1 template molecule/167 droplets (0.003 pg/microL). The frequencies of positive reactions over a range of template concentrations agree closely with the frequencies predicted by Poisson statistics, demonstrating both the accuracy and sensitivity of this platform for limiting dilution and digital PCR applications.

Control and measurement of the phase behavior of aqueous solutions using microfluidics.

A microfluidic device denoted the Phase Chip has been designed to measure and manipulate the phase diagram of multicomponent fluid mixtures. The Phase Chip exploits the permeation of water through poly(dimethylsiloxane) (PDMS) in order to controllably vary the concentration of solutes in aqueous nanoliter volume microdrops stored in wells. The permeation of water in the Phase Chip is modeled using the diffusion equation, and good agreement between experiment and theory is obtained. The Phase Chip operates by first creating drops of the water/solute mixture whose composition varies sequentially. Next, drops are transported down channels and guided into storage wells using surface tension forces. Finally, the solute concentration of each stored drop is simultaneously varied and measured. Two applications of the Phase Chip are presented. First, the phase diagram of a polymer/salt mixture is measured on-chip and validated off-chip, and second, protein crystallization rates are enhanced through the manipulation of the kinetics of nucleation and growth.