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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

Bifurcation of droplet flows within capillaries.

Flows of droplets through networks of microchannels differ significantly from the flow of simple fluids. Our report focuses on the paths of individual droplets through the simplest possible network: a channel that splits into two arms that subsequently recombine. This simple system exhibits complex patterns of flow: both periodic and irregular, depending on the frequency at which the drops are fed into the "loop." A numerical model explains these results and shows regions of regular patterns separated by regions of high complexity. Our results elicit new questions regarding the dynamics of flow of discrete elements of fluids through networks, and point to potential opportunities and difficulties in the design of integrated mini-laboratories operating on droplets.

Effect of high spontaneous polarization on defect structures and orientational dynamics of tilted chiral smectic freely suspended films.

The director structure around topological defects and in 2pi walls in the two-dimensional orientation field of thin freely suspended films of tilted chiral smectic liquid crystal is observed to minimize splay of the spontaneous polarization. Concentric ring patterns in the director field unwind more slowly in higher polarization films. These experiments confirm that polarization space charge increases the effective elasticity of static polarization-splay distortions and that it attracts ionic charge, leading to an increase in the effective orientational viscosity of the director field.

Electric-field-induced transition between the polarization-modulated and ferroelectric smectic- C(S) P(F)* liquid crystalline states studied using microbeam x-ray diffraction.

We report x-ray microbeam studies of a bent-core liquid crystalline material with chiral citronellyl tails. This material has an equilibrium polarization-modulated smectic- CP (PM-SmCP) state exhibiting the B7 texture upon slow cooing from the isotropic while a metastable chiral synclinic ferroelectric Sm-CP state (Sm- C(S) P(*)(F) ) is obtained on quenching from the isotropic. The polarization modulated phase PM-Sm C(S) P(*)(F) shows typical x-ray patterns having multiple satellite peaks around the first-order layer reflection, indicating undulated layers, while the metastable Sm- C(S) P(*)(F) state exhibits a single layering peak indicating flat layers. The Sm- C(S) P(*)(F) state is also induced by the application of an electric field larger than the threshold field ( E(th) ) and thermally returns to the polarization modulated PM-Sm C(S) P(*)(F) structure.

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.

Flexoelectric polarization in hybrid nematic films.

Hybrid films of the nematic liquid crystal penthylcyanobiphenyl (5CB) on the surface of glycerol were studied in the presence of an in-plane electric field. In this geometry there is a macroscopic flexoelectric polarization P(f) with a component in the film plane that couples linearly to the field. Here we report the electric field and film thickness dependence of the width of the 2pi reorientations (2pi walls) in the c director that are trapped by the field and interpret these observations based on continuum elastic theory. We find the difference between the flexoelectric coefficients to be e(1)-e(3)=+11 pC/m.

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 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.

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.