Efficient Genotyping of KRAS Mutant Non-Small Cell Lung Cancer Using a Multiplexed Droplet Digital PCR Approach

© 2015 Pender et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.Droplet digital PCR (ddPCR) can be used to detect low frequency mutations in oncogenedriven lung cancer. The range of KRAS point mutations observed in NSCLC necessitates a multiplex approach to efficient mutation detection in circulating DNA. Here we report the design and optimisation of three discriminatory ddPCR multiplex assays investigating nine different KRAS mutations using PrimePCRddPCRMutation Assays and the Bio-Rad QX100 system. Together these mutations account for 95% of the nucleotide changes found in KRAS in human cancer. Multiplex reactions were optimised on genomic DNA extracted from KRAS mutant cell lines and tested on DNA extracted from fixed tumour tissue from a cohort of lung cancer patients without prior knowledge of the specific KRAS genotype. The multiplex ddPCR assays had a limit of detection of better than 1 mutant KRAS molecule in 2,000 wild-type KRAS molecules, which compared favourably with a limit of detection of 1 in 50 for next generation sequencing and 1 in 10 for Sanger sequencing. Multiplex ddPCR assays thus provide a highly efficient methodology to identify KRAS mutations in lung adenocarcinoma

Plasmodium knowlesi Genome Sequences from Clinical Isolates Reveal Extensive Genomic Dimorphism.

Plasmodium knowlesi is a newly described zoonosis that causes malaria in the human population that can be severe and fatal. The study of P. knowlesi parasites from human clinical isolates is relatively new and, in order to obtain maximum information from patient sample collections, we explored the possibility of generating P. knowlesi genome sequences from archived clinical isolates. Our patient sample collection consisted of frozen whole blood samples that contained excessive human DNA contamination and, in that form, were not suitable for parasite genome sequencing. We developed a method to reduce the amount of human DNA in the thawed blood samples in preparation for high throughput parasite genome sequencing using Illumina HiSeq and MiSeq sequencing platforms. Seven of fifteen samples processed had sufficiently pure P. knowlesi DNA for whole genome sequencing. The reads were mapped to the P. knowlesi H strain reference genome and an average mapping of 90% was obtained. Genes with low coverage were removed leaving 4623 genes for subsequent analyses. Previously we identified a DNA sequence dimorphism on a small fragment of the P. knowlesi normocyte binding protein xa gene on chromosome 14. We used the genome data to assemble full-length Pknbpxa sequences and discovered that the dimorphism extended along the gene. An in-house algorithm was developed to detect SNP sites co-associating with the dimorphism. More than half of the P. knowlesi genome was dimorphic, involving genes on all chromosomes and suggesting that two distinct types of P. knowlesi infect the human population in Sarawak, Malaysian Borneo. We use P. knowlesi clinical samples to demonstrate that Plasmodium DNA from archived patient samples can produce high quality genome data. We show that analyses, of even small numbers of difficult clinical malaria isolates, can generate comprehensive genomic information that will improve our understanding of malaria parasite diversity and pathobiology

Early ctDNA dynamics as a surrogate for progression-free survival in advanced breast cancer in the BEECH trial.

Background Dynamic changes in circulating tumour DNA (ctDNA) levels may predict long-term outcome. We utilised samples from a phase I/II randomised trial (BEECH) to assess ctDNA dynamics as a surrogate for progression-free survival (PFS) and early predictor of drug efficacy.Patients and methods Patients with estrogen receptor-positive advanced metastatic breast cancer (ER+ mBC) in the BEECH study, paclitaxel plus placebo versus paclitaxel plus AKT inhibitor capivasertib, had plasma samples collected for ctDNA analysis at baseline and at multiple time points in the development cohort (safety run-in, part A) and validation cohort (randomised, part B). Baseline sample ctDNA sequencing identified mutations for longitudinal analysis and mutation-specific digital droplet PCR (ddPCR) assays were utilised to assess change in ctDNA abundance (allele fraction) between baseline and 872 on-treatment samples. Primary objective was to assess whether early suppression of ctDNA, based on pre-defined criteria from the development cohort, independently predicted outcome in the validation cohort.Results In the development cohort, suppression of ctDNA was apparent after 8 days of treatment (P = 0.014), with cycle 2 day 1 (4 weeks) identified as the optimal time point to predict PFS from early ctDNA dynamics. In the validation cohort, median PFS was 11.1 months in patients with suppressed ctDNA at 4 weeks and 6.4 months in patients with high ctDNA (hazard ratio = 0.20, 95% confidence interval 0.083-0.50, P < 0.0001). There was no difference in the level of ctDNA suppression between patients randomised to capivasertib or placebo overall (P = 0.904) nor in the PIK3CA mutant subpopulation (P = 0.071). Clonal haematopoiesis of indeterminate potential (CHIP) was evident in 30% (18/59) baseline samples, although CHIP had no effect on tolerance of chemotherapy nor on PFS.Conclusion Early on-treatment ctDNA dynamics are a surrogate for PFS. Dynamic ctDNA assessment has the potential to substantially enhance early drug development.Clinical registration number NCT01625286

Drug-resistant genotypes and multi-clonality in Plasmodium falciparum analysed by direct genome sequencing from peripheral blood of malaria patients.

Naturally acquired blood-stage infections of the malaria parasite Plasmodium falciparum typically harbour multiple haploid clones. The apparent number of clones observed in any single infection depends on the diversity of the polymorphic markers used for the analysis, and the relative abundance of rare clones, which frequently fail to be detected among PCR products derived from numerically dominant clones. However, minority clones are of clinical interest as they may harbour genes conferring drug resistance, leading to enhanced survival after treatment and the possibility of subsequent therapeutic failure. We deployed new generation sequencing to derive genome data for five non-propagated parasite isolates taken directly from 4 different patients treated for clinical malaria in a UK hospital. Analysis of depth of coverage and length of sequence intervals between paired reads identified both previously described and novel gene deletions and amplifications. Full-length sequence data was extracted for 6 loci considered to be under selection by antimalarial drugs, and both known and previously unknown amino acid substitutions were identified. Full mitochondrial genomes were extracted from the sequencing data for each isolate, and these are compared against a panel of polymorphic sites derived from published or unpublished but publicly available data. Finally, genome-wide analysis of clone multiplicity was performed, and the number of infecting parasite clones estimated for each isolate. Each patient harboured at least 3 clones of P. falciparum by this analysis, consistent with results obtained with conventional PCR analysis of polymorphic merozoite antigen loci. We conclude that genome sequencing of peripheral blood P. falciparum taken directly from malaria patients provides high quality data useful for drug resistance studies, genomic structural analyses and population genetics, and also robustly represents clonal multiplicity

Whole Genome Sequencing Highlights Genetic Changes Associated with Laboratory Domestication of C. elegans

Defining the mutational landscape when individuals of a species grow separately and diverge over many generations can provide insights into trait evolution. A specific example of this involves studying changes associated with domestication where different lines of the same wild stock have been cultivated independently in different standard environments. Whole genome sequence comparison of such lines permits estimation of mutation rates, inference of genes' ancestral states and ancestry of existing strains, and correction of sequencing errors in genome databases. Here we study domestication of the C. elegans Bristol strain as a model, and report the genome sequence of LSJ1 (Bristol), a sibling of the standard C. elegans reference wild type N2 (Bristol). The LSJ1 and N2 lines were cultivated separately from shortly after the Bristol strain was isolated until methods to freeze C. elegans were developed. We find that during this time the two strains have accumulated 1208 genetic differences. We describe phenotypic variation between N2 and LSJ1 in the rate at which embryos develop, the rate of production of eggs, the maturity of eggs at laying, and feeding behavior, all the result of post-isolation changes. We infer the ancestral alleles in the original Bristol isolate and highlight 2038 likely sequencing errors in the original N2 reference genome sequence. Many of these changes modify genome annotation. Our study provides a starting point to further investigate genotype-phenotype association and offers insights into the process of selection as a result of laboratory domestication

ConDeTri - A Content Dependent Read Trimmer for Illumina Data

During the last few years, DNA and RNA sequencing have started to play an increasingly important role in biological and medical applications, especially due to the greater amount of sequencing data yielded from the new sequencing machines and the enormous decrease in sequencing costs. Particularly, Illumina/Solexa sequencing has had an increasing impact on gathering data from model and non-model organisms. However, accurate and easy to use tools for quality filtering have not yet been established. We present ConDeTri, a method for content dependent read trimming for next generation sequencing data using quality scores of each individual base. The main focus of the method is to remove sequencing errors from reads so that sequencing reads can be standardized. Another aspect of the method is to incorporate read trimming in next-generation sequencing data processing and analysis pipelines. It can process single-end and paired-end sequence data of arbitrary length and it is independent from sequencing coverage and user interaction. ConDeTri is able to trim and remove reads with low quality scores to save computational time and memory usage during de novo assemblies. Low coverage or large genome sequencing projects will especially gain from trimming reads. The method can easily be incorporated into preprocessing and analysis pipelines for Illumina data

Q&A: ChIP-seq technologies and the study of gene regulation

10.1186/1741-7007-8-56BMC Biology85

Mutation tracking in circulating tumor DNA predicts relapse in early breast cancer.

The identification of early-stage breast cancer patients at high risk of relapse would allow tailoring of adjuvant therapy approaches. We assessed whether analysis of circulating tumor DNA (ctDNA) in plasma can be used to monitor for minimal residual disease (MRD) in breast cancer. In a prospective cohort of 55 early breast cancer patients receiving neoadjuvant chemotherapy, detection of ctDNA in plasma after completion of apparently curative treatment-either at a single postsurgical time point or with serial follow-up plasma samples-predicted metastatic relapse with high accuracy [hazard ratio, 25.1 (confidence interval, 4.08 to 130.5; log-rank P < 0.0001) or 12.0 (confidence interval, 3.36 to 43.07; log-rank P < 0.0001), respectively]. Mutation tracking in serial samples increased sensitivity for the prediction of relapse, with a median lead time of 7.9 months over clinical relapse. We further demonstrated that targeted capture sequencing analysis of ctDNA could define the genetic events of MRD, and that MRD sequencing predicted the genetic events of the subsequent metastatic relapse more accurately than sequencing of the primary cancer. Mutation tracking can therefore identify early breast cancer patients at high risk of relapse. Subsequent adjuvant therapeutic interventions could be tailored to the genetic events present in the MRD, a therapeutic approach that could in part combat the challenge posed by intratumor genetic heterogeneity

Investigation into the annotation of protocol sequencing steps in the sequence read archive

BACKGROUND: The workflow for the production of high-throughput sequencing data from nucleic acid samples is complex. There are a series of protocol steps to be followed in the preparation of samples for next-generation sequencing. The quantification of bias in a number of protocol steps, namely DNA fractionation, blunting, phosphorylation, adapter ligation and library enrichment, remains to be determined. RESULTS: We examined the experimental metadata of the public repository Sequence Read Archive (SRA) in order to ascertain the level of annotation of important sequencing steps in submissions to the database. Using SQL relational database queries (using the SRAdb SQLite database generated by the Bioconductor consortium) to search for keywords commonly occurring in key preparatory protocol steps partitioned over studies, we found that 7.10%, 5.84% and 7.57% of all records (fragmentation, ligation and enrichment, respectively), had at least one keyword corresponding to one of the three protocol steps. Only 4.06% of all records, partitioned over studies, had keywords for all three steps in the protocol (5.58% of all SRA records). CONCLUSIONS: The current level of annotation in the SRA inhibits systematic studies of bias due to these protocol steps. Downstream from this, meta-analyses and comparative studies based on these data will have a source of bias that cannot be quantified at present