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

Validation and implementation of targeted capture and sequencing for the detection of actionable mutation, copy number variation, and gene rearrangement in clinical cancer specimens

Recent years have seen development and implementation of anticancer therapies targeted to particular gene mutations, but methods to assay clinical cancer specimens in a comprehensive way for the critical mutations remain underdeveloped. We have developed UW-OncoPlex, a clinical molecular diagnostic assay to provide simultaneous deep-sequencing information, based on >500× average coverage, for all classes of mutations in 194 clinically relevant genes. To validate UW-OncoPlex, we tested 98 previously characterized clinical tumor specimens from 10 different cancer types, including 41 formalin-fixed paraffin-embedded tissue samples. Mixing studies indicated reliable mutation detection in samples with ≥ 10% tumor cells. In clinical samples with ≥ 10% tumor cells, UW-OncoPlex correctly identified 129 of 130 known mutations [sensitivity 99.2%, (95% CI, 95.8%-99.9%)], including single nucleotide variants, small insertions and deletions, internal tandem duplications, gene copy number gains and amplifications, gene copy losses, chromosomal gains and losses, and actionable genomic rearrangements, including ALK-EML4, ROS1, PML-RARA, and BCR-ABL. In the same samples, the assay also identified actionable point mutations in genes not previously analyzed and novel gene rearrangements of MLL and GRIK4 in melanoma, and of ASXL1, PIK3R1, and SGCZ in acute myeloid leukemia. To best guide existing and emerging treatment regimens and facilitate integration of genomic testing with patient care, we developed a framework for data analysis, decision support, and reporting clinically actionable results.

Next-generation sequencing-based multigene mutational screening for acute myeloid leukemia using MiSeq: applicability for diagnostics and disease monitoring

Routine molecular testing in acute myeloid leukemia involves screening several genes of therapeutic and prognostic significance for mutations. A comprehensive analysis using single-gene assays requires large amounts of DNA, is cumbersome and timely consolidation of results for clinical reporting is challenging. High throughput, next-generation sequencing platforms widely used in research have not been tested vigorously for clinical application. Here we describe the clinical application of MiSeq, a next-generation sequencing platform to screen mutational hotspots in 54 cancer-related genes including genes relevant in acute myeloid leukemia (NRAS, KRAS, FLT3, NPM1, DNMT3A, IDH1/2, JAK2, KIT and EZH2). We sequenced 63 samples from patients with acute myeloid leukemia/myelodysplastic syndrome using MiSeq and compared the results with those obtained using another next-generation sequencing platform, Ion-Torrent Personal Genome Machine and other conventional testing platforms. MiSeq detected a total of 100 single nucleotide variants and 23 NPM1 insertions that were confirmed by Ion Torrent or conventional platforms, indicating complete concordance. FLT3-internal tandem duplications (n=10) were not detected; however, re-analysis of the MiSeq output by Pindel, an indel detection algorithm, did detect them. Dilution studies of cancer cell-line DNA showed that the quantitative accuracy of mutation detection was up to an allelic frequency of 1.5% with a high level of inter- and intra-run assay reproducibility, suggesting potential utility for monitoring response to therapy, clonal heterogeneity and evolution. Examples demonstrating the advantages of MiSeq over conventional platforms for disease monitoring are provided. Easy work-flow, high throughput multiplexing capability, 4-day turnaround time and simultaneous assessment of routinely tested and emerging markers make MiSeq highly applicable for clinical molecular testing in acute myeloid leukemia.

Next generation analysis of breast cancer genomes for precision medicine

For many years breast cancer classification has been based on histology and immune-histochemistry. New techniques, more strictly related to cancer biology, partially succeeded in fractionating patients, correlated to survival and better predicted the patient response to therapy. Nowadays, great expectations arise from massive parallel or high throughput next generation sequencing. Cancer genomics has already revolutionized our knowledge of breast cancer molecular pathology, paving the way to the development of new and more effective clinical protocols. This review is focused on the most recent advances in the field of cancer genomics and epigenomics, including DNA alterations and driver gene mutations, gene fusions, DNA methylation and miRNA expression.

Desktop transcriptome sequencing from archival tissue to identify clinically relevant translocations

Somatic mutations, often translocations or single nucleotide variations, are pathognomonic for certain types of cancers and are increasingly of clinical importance for diagnosis and prediction of response to therapy. Conventional clinical assays only evaluate 1 mutation at a time, and targeted tests are often constrained to identify only the most common mutations. Genome-wide or transcriptome-wide high-throughput sequencing (HTS) of clinical samples offers an opportunity to evaluate for all clinically significant mutations with a single test. Recently a "desktop version" of HTS has become available, but most of the experience to date is based on data obtained from high-quality DNA from frozen specimens. In this study, we demonstrate, as a proof of principle, that translocations in sarcomas can be diagnosed from formalin-fixed paraffin-embedded (FFPE) tissue with desktop HTS. Using the first generation MiSeq platform, full transcriptome sequencing was performed on FFPE material from archival blocks of 3 synovial sarcomas, 3 myxoid liposarcomas, 2 Ewing sarcomas, and 1 clear cell sarcoma. Mapping the reads to the "sarcomatome" (all known 83 genes involved in translocations and mutations in sarcoma) and using a novel algorithm for ranking fusion candidates, the pathognomonic fusions and the exact breakpoints were identified in all cases of synovial sarcoma, myxoid liposarcoma, and clear cell sarcoma. The Ewing sarcoma fusion gene was detectable in FFPE material only with a sequencing platform that generates greater sequencing depth. The results show that a single transcriptome HTS assay, from FFPE, has the potential to replace conventional molecular diagnostic techniques for the evaluation of clinically relevant mutations in cancer.

Towards clinical molecular diagnosis of inherited cardiac conditions: a comparison of bench-top genome DNA sequencers

BACKGROUND

Molecular genetic testing is recommended for diagnosis of inherited cardiac disease, to guide prognosis and treatment, but access is often limited by cost and availability. Recently introduced high-throughput bench-top DNA sequencing platforms have the potential to overcome these limitations.

METHODOLOGY/PRINCIPAL FINDINGS

We evaluated two next-generation sequencing (NGS) platforms for molecular diagnostics. The protein-coding regions of six genes associated with inherited arrhythmia syndromes were amplified from 15 human samples using parallelised multiplex PCR (Access Array, Fluidigm), and sequenced on the MiSeq (Illumina) and Ion Torrent PGM (Life Technologies). Overall, 97.9% of the target was sequenced adequately for variant calling on the MiSeq, and 96.8% on the Ion Torrent PGM. Regions missed tended to be of high GC-content, and most were problematic for both platforms. Variant calling was assessed using 107 variants detected using Sanger sequencing: within adequately sequenced regions, variant calling on both platforms was highly accurate (Sensitivity: MiSeq 100%, PGM 99.1%. Positive predictive value: MiSeq 95.9%, PGM 95.5%). At the time of the study the Ion Torrent PGM had a lower capital cost and individual runs were cheaper and faster. The MiSeq had a higher capacity (requiring fewer runs), with reduced hands-on time and simpler laboratory workflows. Both provide significant cost and time savings over conventional methods, even allowing for adjunct Sanger sequencing to validate findings and sequence exons missed by NGS.

CONCLUSIONS/SIGNIFICANCE

MiSeq and Ion Torrent PGM both provide accurate variant detection as part of a PCR-based molecular diagnostic workflow, and provide alternative platforms for molecular diagnosis of inherited cardiac conditions. Though there were performance differences at this throughput, platforms differed primarily in terms of cost, scalability, protocol stability and ease of use. Compared with current molecular genetic diagnostic tests for inherited cardiac arrhythmias, these NGS approaches are faster, less expensive, and yet more comprehensive.

Mutascope: sensitive detection of somatic mutations from deep amplicon sequencing

Summary: We present Mutascope, a sequencing analysis pipeline specifically developed for the identification of somatic variants present at low-allelic fraction from high-throughput sequencing of amplicons from matched tumor-normal specimen. Using datasets reproducing tumor genetic heterogeneity, we demonstrate that Mutascope has a higher sensitivity and generates fewer false-positive calls than tools designed for shotgun sequencing or diploid genomes.

Availability: Freely available on the web at http://sourceforge.net/projects/mutascope/.

Contact:oharismendy@ucsd.edu

Supplementary information:Supplementary data are available at Bioinformatics online.

Reducing sequence artifacts in amplicon-based massively parallel sequencing of formalin-fixed paraffin-embedded DNA by enzymatic depletion of uracil-containing templates

BACKGROUND

Formalin-fixed, paraffin-embedded (FFPE) tissues are routinely used for detecting mutational biomarkers in patients with cancer. A previous intractable challenge with FFPE DNA in genetic testing has been the high number of artifactual single-nucleotide changes (SNCs), particularly for the detection of low-level mutations. Pretreatment of FFPE DNA with uracil-DNA glycosylase (UDG) can markedly reduce these C:G>T:A SNCs with a small panel of amplicons. This procedure has implications for massively parallel sequencing approaches to mutation detection from DNA. We investigated whether sequence artifacts were problematic in amplicon-based massively parallel sequencing and what effect UDG pretreatment had on reducing these artifacts.

METHODS

We amplified selected amplicons from lung cancer FFPE DNAs using the TruSeq Cancer Panel. SNCs occurring at a frequency <10% were considered most likely to represent sequence artifacts and were enumerated for both UDG-treated and -untreated DNAs.

RESULTS

Massively parallel sequencing of FFPE DNA samples showed multiple SNCs, predominantly C:G>T:A changes, with a significant proportion occurring above the background sequencing error (defined as 1%). UDG pretreatment markedly reduced C:G>T:A SNCs without affecting the detection of true somatic mutations. However, C:G>T:A changes within CpG dinucleotides were often resistant to the UDG treatment as a consequence of 5-methyl cytosine being deaminated to thymine rather than uracil.

CONCLUSIONS

UDG pretreatment greatly facilitates the accurate discrimination of mutations in FFPE samples by use of amplicon-based approaches. This is particularly important when working with samples with low tumor purity or for the assessment of mutational heterogeneity in tumors.

Heterogeneity of colorectal cancer (CRC) in reference to KRAS proto-oncogene utilizing WAVE technology

BACKGROUND

New drugs targeting specific genes required for unregulated growth and metastases have improved survival rates for patients with metastatic colorectal cancer. Resistance to monoclonal antibodies specific for the epidermal growth factor receptor (EGFR) has been attributed to the presence of activating point mutations in the proto-oncogene KRAS. The use of EGFR inhibitor monotherapy in patients that have KRAS wild type has produced response rates of only 10-20%. The molecular basis for clinical resistance remains poorly understood. We propose two possible explanations to explain these low response rates; 1) levels of resistant CRC cells carrying mutated KRAS are below the sensitivity of standard direct sequencing modalities (<5%) or 2) the standard practice of analyzing a single area within a heterogeneous tumor is a practice that can overlook areas with mutated KRAS.

METHODS

In a collaborative effort with the surgical and molecular pathology departments, 3 formalin fixed paraffin embedded tissue blocks of human CRC were obtained from the human tissue bank maintained by the Lifespan Pathology Department and/or the human tissue bank maintained by the Molecular Pathology Core of the COBRE for Cancer Research Development. The three specimens previously demonstrated KRAS mutations detected by the Applied Biosystems Kit. The Wave system 4500 (high performance ion-pairing liquid chromatography (IP-HPLC)) was utilized to evaluate tissue for the presence of KRAS proto-oncogene mutations at codons 12 and 13.

RESULTS

Initially, the sensitivity of WAVE technology was compared with direct sequencing by evaluating a dilutional series. WAVE detected mutant alleles at levels of 2.5% compared to 20% performed with standard direct sequencing. Samples from three patients were evaluated by WAVE technology. Eight samples from patient 1 were analyzed. In two of eight samples, no mutations were detected at concentrations as low as 5%. In one sample a mutation was noted by WAVE and not by direct sequencing. All four samples from patient 2 tested positive for Exon 12/13 mutations. Of the seven samples from patient 3, five were positive for Exon 12/13 mutations and two were negative for Exon 12/13 mutations.

CONCLUSION

In these studies the analysis of three patients' colorectal cancer tissues were analyzed utilizing the WAVE technology. Results demonstrated a greater degree of sensitivity in mutation detection when compared to standard sequencing. These studies also demonstrated heterogeneity of expression of KRAS mutations between areas of the tissue samples at a genomic level. The low clinical response rates to EGFR inhibition might be explained by the variation in mutation presence, which was dependent upon the region examined. The heterogeneity demonstrated in these studies provides another phenotypic variant that will impact clinical care.

Single molecule molecular inversion probes for targeted, high-accuracy detection of low-frequency variation

The detection and quantification of genetic heterogeneity in populations of cells is fundamentally important to diverse fields, ranging from microbial evolution to human cancer genetics. However, despite the cost and throughput advances associated with massively parallel sequencing, it remains challenging to reliably detect mutations that are present at a low relative abundance in a given DNA sample. Here we describe smMIP, an assay that combines single molecule tagging with multiplex targeted capture to enable practical and highly sensitive detection of low-frequency or subclonal variation. To demonstrate the potential of the method, we simultaneously resequenced 33 clinically informative cancer genes in eight cell line and 45 clinical cancer samples. Single molecule tagging facilitated extremely accurate consensus calling, with an estimated per-base error rate of 8.4 × 10(-6) in cell lines and 2.6 × 10(-5) in clinical specimens. False-positive mutations in the single molecule consensus base-calls exhibited patterns predominantly consistent with DNA damage, including 8-oxo-guanine and spontaneous deamination of cytosine. Based on mixing experiments with cell line samples, sensitivity for mutations above 1% frequency was 83% with no false positives. At clinically informative sites, we identified seven low-frequency point mutations (0.2%-4.7%), including BRAF p.V600E (melanoma, 0.2% alternate allele frequency), KRAS p.G12V (lung, 0.6%), JAK2 p.V617F (melanoma, colon, two lung, 0.3%-1.4%), and NRAS p.Q61R (colon, 4.7%). We anticipate that smMIP will be broadly adoptable as a practical and effective method for accurately detecting low-frequency mutations in both research and clinical settings.

Variable clonal repopulation dynamics influence chemotherapy response in colorectal cancer

Intratumoral heterogeneity arises through the evolution of genetically diverse subclones during tumor progression. However, it remains unknown whether cells within single genetic clones are functionally equivalent. By combining DNA copy number alteration (CNA) profiling, sequencing, and lentiviral lineage tracking, we followed the repopulation dynamics of 150 single lentivirus-marked lineages from 10 human colorectal cancers through serial xenograft passages in mice. CNA and mutational analysis distinguished individual clones and showed that clones remained stable upon serial transplantation. Despite this stability, the proliferation, persistence, and chemotherapy tolerance of lentivirally marked lineages were variable within each clone. Chemotherapy promoted the dominance of previously minor or dormant lineages. Thus, apart from genetic diversity, tumor cells display inherent functional variability in tumor propagation potential, which contributes to both cancer growth and therapy tolerance.

Targeted, high-depth, next-generation sequencing of cancer genes in formalin-fixed, paraffin-embedded and fine-needle aspiration tumor specimens

Implementation of highly sophisticated technologies, such as next-generation sequencing (NGS), into routine clinical practice requires compatibility with common tumor biopsy types, such as formalin-fixed, paraffin-embedded (FFPE) and fine-needle aspiration specimens, and validation metrics for platforms, controls, and data analysis pipelines. In this study, a two-step PCR enrichment workflow was used to assess 540 known cancer-relevant variants in 16 oncogenes for high-depth sequencing in tumor samples on either mature (Illumina GAIIx) or emerging (Ion Torrent PGM) NGS platforms. The results revealed that the background noise of variant detection was elevated approximately twofold in FFPE compared with cell line DNA. Bioinformatic algorithms were optimized to accommodate this background. Variant calls from 38 residual clinical colorectal cancer FFPE specimens and 10 thyroid fine-needle aspiration specimens were compared across multiple cancer genes, resulting in an accuracy of 96.1% (95% CI, 96.1% to 99.3%) compared with Sanger sequencing, and 99.6% (95% CI, 97.9% to 99.9%) compared with an alternative method with an analytical sensitivity of 1% mutation detection. A total of 45 of 48 samples were concordant between NGS platforms across all matched regions, with the three discordant calls each represented at <10% of reads. Consequently, NGS of targeted oncogenes in real-life tumor specimens using distinct platforms addresses unmet needs for unbiased and highly sensitive mutation detection and can accelerate both basic and clinical cancer research.

Next-generation sequencing - feasibility and practicality in haematology

Next-generation sequencing platforms have evolved to provide an accurate and comprehensive means for the detection of molecular mutations in heterogeneous tumour specimens. Here, we review the feasibility and practicality of this novel laboratory technology. In particular, we focus on the utility of next-generation sequencing technology in characterizing haematological neoplasms and the landmark findings in key haematological malignancies. We also discuss deep-sequencing strategies to analyse the constantly increasing number of molecular markers applied for disease classification, patient stratification and individualized monitoring of minimal residual disease. Although many facets of this assay need to be taken into account, amplicon deep-sequencing has already demonstrated a promising technical performance and is being continuously developed towards routine application in diagnostic laboratories so that an impact on clinical practice can be achieved.

High-throughput sequencing of the melanoma genome

Next-generation sequencing technologies are now common for whole-genome, whole-exome and whole-transcriptome sequencing (RNA-seq) of tumors to identify point mutations, structural or copy number alterations and changes in gene expression. A substantial number of studies have already been performed for melanoma. One study analysed eight melanoma cell lines with RNA-Seq technology and identified 11 novel melanoma gene fusions. Whole-exome sequencing of seven melanoma cell lines identified overlapping gain of function mutations in MAP2K1 (MEK1) and MAP2K2 (MEK2) genes. Integrative sequencing of cutaneous melanoma metastases using different sequencing platforms revealed a new somatic point mutation in HRAS and a structural rearrangement affecting CDKN2C (a CDK4 inhibitor). These latter sequencing-based discoveries may be used to motivate the inclusion of the affected patients into clinical trials with specific signalling pathway inhibitors. Taken together, we are at the beginning of an era with new sequencing technologies providing a more comprehensive view of cancer mutational landscapes and hereby a better understanding of their pathogenesis. This will also open interesting perspectives for new treatment approaches and clinical trial designs.

LoFreq: a sequence-quality aware, ultra-sensitive variant caller for uncovering cell-population heterogeneity from high-throughput sequencing datasets

The study of cell-population heterogeneity in a range of biological systems, from viruses to bacterial isolates to tumor samples, has been transformed by recent advances in sequencing throughput. While the high-coverage afforded can be used, in principle, to identify very rare variants in a population, existing ad hoc approaches frequently fail to distinguish true variants from sequencing errors. We report a method (LoFreq) that models sequencing run-specific error rates to accurately call variants occurring in <0.05% of a population. Using simulated and real datasets (viral, bacterial and human), we show that LoFreq has near-perfect specificity, with significantly improved sensitivity compared with existing methods and can efficiently analyze deep Illumina sequencing datasets without resorting to approximations or heuristics. We also present experimental validation for LoFreq on two different platforms (Fluidigm and Sequenom) and its application to call rare somatic variants from exome sequencing datasets for gastric cancer. Source code and executables for LoFreq are freely available at http://sourceforge.net/projects/lofreq/.

Mutation detection by clonal sequencing of PCR amplicons and grouped read typing is applicable to clinical diagnostics

We describe a sensitive technique for mutation detection using clonal sequencing. We analyzed DNA extracted from 13 cancer cell lines and 35 tumor samples and applied a novel approach to identify disease-associated somatic mutations. By matching reads against an index of known variants, noise can be dramatically reduced, enabling the detection and quantification of those variants, even when they are present at less than 1% of the total sequenced population; this is comparable to, or better than, current diagnostic methods. Following the identification or exclusion of known variants, unmatched reads are grouped for BLAST searching to identify novel variants or contaminants. Known variants, novel variants, and contaminants were readily identified in tumor tissue using this approach. Our approach also enables an estimation of the per-base sequencing error rate, providing a confidence threshold for interpretation of the results in the clinic. This novel approach has immediate applicability to clinical testing for disease-associated genetic variants.

Improving indel detection specificity of the Ion Torrent PGM benchtop sequencer

The emergence of benchtop sequencers has made clinical genetic testing using next-generation sequencing more feasible. Ion Torrent's PGM^TM is one such benchtop sequencer that shows clinical promise in detecting single nucleotide variations (SNVs) and microindel variations (indels). However, the large number of false positive indels caused by the high frequency of homopolymer sequencing errors has impeded PGM^TM's usage for clinical genetic testing. An extensive analysis of PGM^TM data from the sequencing reads of the well-characterized genome of the Escherichia coli DH10B strain and sequences of the BRCA1 and BRCA2 genes from six germline samples was done. Three commonly used variant detection tools, SAMtools, Dindel, and GATK's Unified Genotyper, all had substantial false positive rates for indels. By incorporating filters on two major measures we could dramatically improve false positive rates without sacrificing sensitivity. The two measures were: B-Allele Frequency (BAF) and VARiation of the Width of gaps and inserts (VARW) per indel position. A BAF threshold applied to indels detected by UnifiedGenotyper removed ∼99% of the indel errors detected in both the DH10B and BRCA sequences. The optimum BAF threshold for BRCA sequences was determined by requiring 100% detection sensitivity and minimum false discovery rate, using variants detected from Sanger sequencing as reference. This resulted in 15 indel errors remaining, of which 7 indel errors were removed by selecting a VARW threshold of zero. VARW specific errors increased in frequency with higher read depth in the BRCA datasets, suggesting that homopolymer-associated indel errors cannot be reduced by increasing the depth of coverage. Thus, using a VARW threshold is likely to be important in reducing indel errors from data with higher coverage. In conclusion, BAF and VARW thresholds provide simple and effective filtering criteria that can improve the specificity of indel detection in PGM^TM data without compromising sensitivity.

Noninvasive identification and monitoring of cancer mutations by targeted deep sequencing of plasma DNA

Plasma of cancer patients contains cell-free tumor DNA that carries information on tumor mutations and tumor burden. Individual mutations have been probed using allele-specific assays, but sequencing of entire genes to detect cancer mutations in circulating DNA has not been demonstrated. We developed a method for tagged-amplicon deep sequencing (TAm-Seq) and screened 5995 genomic bases for low-frequency mutations. Using this method, we identified cancer mutations present in circulating DNA at allele frequencies as low as 2%, with sensitivity and specificity of >97%. We identified mutations throughout the tumor suppressor gene TP53 in circulating DNA from 46 plasma samples of advanced ovarian cancer patients. We demonstrated use of TAm-Seq to noninvasively identify the origin of metastatic relapse in a patient with multiple primary tumors. In another case, we identified in plasma an EGFR mutation not found in an initial ovarian biopsy. We further used TAm-Seq to monitor tumor dynamics, and tracked 10 concomitant mutations in plasma of a metastatic breast cancer patient over 16 months. This low-cost, high-throughput method could facilitate analysis of circulating DNA as a noninvasive "liquid biopsy" for personalized cancer genomics.