Amplicon sequencing of colorectal cancer: variant calling in frozen and formalin-fixed samples

Clinical Bioinformatician Body of Knowledge-Bioinformatics and Software Core: A Report of the Association for Molecular Pathology

With the evolution of next-generation sequencing-based testing in molecular diagnostics laboratories, the clinical role of bioinformaticians has also evolved. The Association for Molecular Pathology's Clinical Bioinformatician Body of Knowledge aims to define the various roles the clinical bioinformatician operates individually or within a clinical bioinformatics team, along with proficiencies and skill sets that may be required or desirable across these roles. One of the most common professional responsibilities of a clinical bioinformatician is to implement bioinformatics pipelines, either vendor supplied or custom built for the assays in the molecular diagnostics laboratory, along with analysis and quality control of clinical genomics data. This second article in the series describes the various stages in the life cycle of a clinical bioinformatics pipeline and the considerations, areas of expertise, and skill sets required in each stage. This information may help laboratory professionals to better work with clinical bioinformaticians and laboratory directors to hire the appropriate expertise based on the specific needs of the laboratory.

Comparison of next‑generation sequencing quality metrics and concordance in the detection of cancer‑specific molecular alterations between formalin‑fixed paraffin‑embedded and fresh‑frozen samples in comprehensive genomic profiling with the Illumina® TruSight Oncology 500 assay

Next-generation sequencing (NGS) technology is routinely employed to detect clinically significant variants in the field of precision medicine. Formalin-fixed paraffin-embedded (FFPE) tissues remain a widely used source of genetic material for diagnostic purposes due to their long-term storage stability, which preserves the architecture of tumour tissues. However, the degradation of nucleic acids (NAs) that occurs during the fixation process can lead to unreliable results or hinder analysis. Considering the challenges posed by the quality and quantity of nucleic acids extracted from FFPE samples, the present study aimed to compare paired fresh-frozen (FF) and FFPE samples. The present study also aimed to assess the concordance rate of quality control metrics of sequenced libraries and identify variants and biomarkers between FFPE and FF tissues. The Illumina® TruSight Oncology 500 assay (TSO 500; Illumina, Inc.) was used to conduct comprehensive genomic profiling of samples from patients with various cancer types. All identified alterations were annotated using the Pierian Dx Clinical Genomics Workspace version 6.20 (PierianDx CGW). A total of 138 DNA and 138 RNA analyses on 69 paired samples were performed. These results demonstrated the significant potential of FF tissue as a primary source of higher-quality genetic material to detect small variants, microsatellite instabilities and the tumour mutational burden, using the TSO 500 assay compared with FFPE samples. This approach could effectively reduce the issues associated with poor-quality NAs extracted from FFPE samples. Based on our findings of lower concordance in the detection of splice variants, fusions and copy number variants in paired samples, we recommend that future studies using the TSO 500 assay should focus directly on this issue.

Nadir paranazal sinüs kanserlerinde yeni tanımlanan reseptör tirozin kinaz mutasyonları ve potansiyel fonksiyonel etkileri

Amaç: Paranazal sinüs kanserleri oldukça nadir görülen heterojen bir hastalık grubudur. Maksiler sinüs skuamoz hücreli karsinomu, paranazal sinüs kanserlerinin anatomik ve histolojik olarak en yaygın alt tipidir. Bu kanserin genetik profiline dair bilginin sınırlı olması, hastaların hedefli tedavi seçeneklerinden yararlanamamasına neden olmaktadır. Çalışmamızda bu nadir kanserdeki reseptör tirozin kinaz mutasyonlarının tanımlanması ve mutasyonların olası fonksiyonel etkilerinin tahmin edilmesi amaçlanmıştır. Gereç ve Yöntem: Bu amaçla 30 olgunun tümörüne ait FFPE dokulardan DNA izolasyonu gerçekleştirildi, olguların mutasyon profili yeni nesil sekanslama yöntemi ve biyoinformatik değerlendirme ile belirlendi. Belirlenen patojenik/ olası patojenik varyantların fonksiyonel etkileri farklı in silico araçlar yardımıyla tahminlendi. Bulgular: Olgularının tamamında en az bir adet patojenik/olası patojenik KIT, PDFGRA ve RETmutasyonu belirlendi. KIT geninin katalitik bölgesindeki mutasyonların kinaz aktivitesini arttıracağı tahmin edildi. PDFGRA genindeki p.P567P ve p.D1074D mutasyonları, 30 olgunun tamamında ve SRA veritabanından elde edilen normal dokulara ait okumaların tümünde belirlendi. Sonuç: Reseptör tirozin kinaz mutasyonlarının paranazal sinüs kanserlerinde de önemli rol oynayabileceğinin belirlenmiş olması özellikle artmış kinaz aktivitesini hedefleyen tedavi yaklaşımlarını bu olguların erişimine sunma potansiyeli taşıması bakımından oldukça önemlidir.

Minimizing Sample Failure Rates for Challenging Clinical Tumor Samples

Identification of somatic variants in cancer by high-throughput sequencing has become common clinical practice, largely because many of these variants may be predictive biomarkers for targeted therapies. However, there can be high sample quality control (QC) failure rates for some tests that prevent the return of results. Stem-loop inhibition mediated amplification (SLIMamp) is a patented technology that has been incorporated into commercially available cancer next-generation sequencing testing kits. The claimed advantage is that these kits can interrogate challenging formalin-fixed, paraffin-embedded tissue samples with low tumor purity, poor-quality DNA, and/or low-input DNA, resulting in a high sample QC pass rate. The study aimed to substantiate that claim using Pillar Biosciences oncoReveal Solid Tumor Panel. Forty-eight samples that had failed one or more preanalytical QC sample parameters for whole-exome sequencing from the Australian Translational Genomics Center's ISO15189-accredited diagnostic genomics laboratory were acquired. XING Genomic Services performed an exploratory data analysis to characterize the samples and then tested the samples in their ISO15189-accredited laboratory. Clinical reports could be generated for 37 (77%) samples, of which 29 (60%) contained clinically actionable or significant variants that would not otherwise have been identified. Eleven samples were deemed unreportable, and the sequencing data were likely dominated by artifacts. A novel postsequencing QC metric was developed that can discriminate between clinically reportable and unreportable samples.

Performance of 16S Metagenomic Profiling in Formalin-Fixed Paraffin-Embedded versus Fresh-Frozen Colorectal Cancer Tissues

Formalin-fixed, paraffin-embedded (FFPE) tissues represent the most widely available clinical material to study colorectal cancer (CRC). However, the accuracy and clinical validity of FFPE microbiome profiling in CRC is uncertain. Here, we compared the microbial composition of 10 paired fresh-frozen (FF) and FFPE CRC tissues using 16S rRNA sequencing and RNA-ISH. Both sample types showed different microbial diversity and composition. FF samples were enriched in archaea and representative CRC-associated bacteria, such as Firmicutes, Bacteroidetes and Fusobacteria. Conversely, FFPE samples were mainly enriched in typical contaminants, such as Sphingomonadales and Rhodobacterales. RNA-ISH in FFPE tissues confirmed the presence of CRC-associated bacteria, such as Fusobacterium and Bacteroides, as well as Propionibacterium allowing discrimination between tumor-associated and contaminant taxa. An internal quality index showed that the degree of similarity within sample pairs inversely correlated with the dominance of contaminant taxa. Given the importance of FFPE specimens for larger studies in human cancer genomics, our findings may provide useful indications on potential confounding factors to consider for accurate and reproducible metagenomics analyses.

The Mutational Concordance of Fixed Formalin Paraffin Embedded and Fresh Frozen Gastro-Oesophageal Tumours Using Whole Exome Sequencing

1.

BACKGROUND

The application of massively parallel sequencing has led to the identification of aberrant druggable pathways and somatic mutations within therapeutically relevant genes in gastro-oesophageal cancer. Given the widespread use of formalin-fixed paraffin-embedded (FFPE) samples in the study of this disease, it would be beneficial, especially for the purposes of biomarker evaluation, to assess the concordance between comprehensive exome-wide sequencing data from archival FFPE samples originating from a prospective clinical study and those derived from fresh-frozen material. 2.

METHODS

We analysed whole-exome sequencing data to define the mutational concordance of 16 matched fresh-frozen and FFPE gastro-oesophageal tumours (N = 32) from a prospective clinical study. We assessed DNA integrity prior to sequencing and then identified coding mutations in genes that have previously been implicated in other cancers. In addition, we calculated the mutant-allele heterogeneity (MATH) for these samples. 3.

RESULTS

Although there was increased degradation of DNA in FFPE samples compared with frozen samples, sequencing data from only two FFPE samples failed to reach an adequate mapping quality threshold. Using a filtering threshold of mutant read counts of at least ten and a minimum of 5% variant allele frequency (VAF) we found that there was a high median mutational concordance of 97% (range 80.1-98.68%) between fresh-frozen and FFPE gastro-oesophageal tumour-derived exomes. However, the majority of FFPE tumours had higher mutant-allele heterogeneity (MATH) scores when compared with corresponding frozen tumours (p < 0.001), suggesting that FFPE-based exome sequencing is likely to over-represent tumour heterogeneity in FFPE samples compared to fresh-frozen samples. Furthermore, we identified coding mutations in 120 cancer-related genes, including those associated with chromatin remodelling and Wnt/β-catenin and Receptor Tyrosine Kinase signalling. 4.

CONCLUSIONS

These data suggest that comprehensive genomic data can be generated from exome sequencing of selected DNA samples extracted from archival FFPE gastro-oesophageal tumour tissues within the context of prospective clinical trials.

Why Formalin-fixed, Paraffin-embedded Biospecimens Must Be Used in Genomic Medicine: An Evidence-based Review and Conclusion

Fresh-frozen tissue is the “gold standard” biospecimen type for next-generation sequencing (NGS). However, collecting frozen tissue is usually not feasible because clinical workflows deliver formalin-fixed, paraffin-embedded (FFPE) tissue blocks. Some clinicians and researchers are reticent to embrace the use of FFPE tissue for NGS because FFPE tissue can yield low quantities of degraded DNA, containing formalin-induced mutations. We describe the process by which formalin-induced deamination can lead to artifactual cytosine (C) to thymine (T) and guanine (G) to adenine (A) (C:G > T:A) mutation calls and perform a literature review of 17 publications that compare NGS data from patient-matched fresh-frozen and FFPE tissue blocks. We conclude that although it is indeed true that sequencing data from FFPE tissue can be poorer than those from frozen tissue, any differences occur at an inconsequential magnitude, and FFPE biospecimens can be used in genomic medicine with confidence:

Pisces: an accurate and versatile variant caller for somatic and germline next-generation sequencing data

Motivation

Next-generation sequencing technology is transitioning quickly from research labs to clinical settings. The diagnosis and treatment selection for many acquired and autosomal conditions necessitate a method for accurately detecting somatic and germline variants.

Results

We have developed Pisces, a rapid, versatile and accurate small-variant calling suite designed for somatic and germline amplicon sequencing applications. Accuracy is achieved by four distinct modules, each incorporating a number of novel algorithmic strategies.

Availability and implementation

Pisces is distributed under an open source license and can be downloaded from https://github.com/Illumina/Pisces. Pisces is available on the BaseSpace™ SequenceHub. It is distributed on Illumina sequencing platforms such as the MiSeq™ and is included in the Praxis™ Extended RAS Panel test which was recently approved by the FDA.

Supplementary information

Supplementary data are available at Bioinformatics online.

Extended Molecular Profiling Improves Stratification and Prediction of Survival After Resection of Colorectal Liver Metastases

OBJECTIVE

The aim of this study was to assess the effect of cancer-related genes and their mutations analyzed by next-generation sequencing (NGS) on the oncological outcome after resection of colorectal liver metastases (CRLM).

BACKGROUND

Traditional prognostic scores include clinical and pathological parameters of primary tumor and metastases. The modified clinical risk score (m-CS), based on size of metastases, primary tumor nodal status, and RAS mutation status outperformed traditional scores. We hypothesized to further improve the scoring system based on the results of NGS.

METHODS

Cancer tissues of 139 patients with CRLM were used for NGS. The work-up included the analysis of recurrent somatic mutations and copy number changes of 720 genes. Clinical data were extracted from a prospectively collected institutional liver database.

RESULTS

Depending on significance, the following cancer-related genes and their alterations (%) were further investigated: APC (86%), TP53 (78%), KRAS (29%), SMAD4 (15%), PIK3CA (14%), BRAF (8%), ERBB2 (6%), SMAD3 (5%), SMAD2 (4%), and NRAS (4%). The most predictive parameters for poor oncological outcome were alterations in the SMAD family (P = 0.0186) and RAS-RAF pathway (P = 0.032). Refining the m-CS by replacing RAS with RAS-RAF pathway and adding SMAD family resulted in an extended clinical risk score which is highly predictive for oncological outcome (P < 0.0001).

CONCLUSION

In conclusion, mutations of the SMAD family revealed a strong prognostic effect after surgery for CRLM. Integration of alterations of the SMAD family as well as the RAS/RAF pathway resulted in a new, still simple but highly prognostic score.

Comparison of Fresh Frozen Tissue With Formalin-Fixed Paraffin-Embedded Tissue for Mutation Analysis Using a Multi-Gene Panel in Patients With Colorectal Cancer

Background: Next generation sequencing (NGS)-based multi-gene panel tests have been performed to predict the treatment response and prognosis in patients with colorectal cancer (CRC). Whether the multi-gene mutation results of formalin-fixed paraffin-embedded (FFPE) tissues are identical to those of fresh frozen tissues remains unknown. Methods: A 22-gene panel with 103 hotspots was used to detect mutations in paired fresh frozen tissue and FFPE tissue from 118 patients with CRC. Results: In our study, 117 patients (99.2%) had one or more variants, with 226 variants in FFPE tissue and 221 in fresh frozen tissue. Of the 129 variants identified in this study, 96 variants were present in both FFPE and fresh frozen tissues; 27 variants were found in FFPE tissues only; 6 variants were found only in fresh frozen tissues. The mutation results demonstrated >94.0% concordance in all variants, with Kappa coefficient >0.500 in 64.3% (83/129) of variants. At the gene level, concordance ranged from 73.8 to 100.0%, with Kappa coefficient >0.500 in 81.3% (13/16) of genes. Conclusions: The results of mutation analysis performed with a multi-gene panel and FFPE and fresh frozen tissue were highly concordant in patients with CRC, at both the variant and gene levels. There were, however, some important differences in mutation results between the two tissue types. Therefore, fresh frozen tissue should not routinely be replaced with FFPE tissue for mutation analysis with a multi-gene panel. Rather, FFPE tissue is a reasonable alternative for fresh frozen tissue when the latter is unavailable.

Systematic comparison of somatic variant calling performance among different sequencing depth and mutation frequency

In the past decade, treatments for tumors have made remarkable progress, such as the successful clinical application of targeted therapies. Nowadays, targeted therapies are based primarily on the detection of mutations, and next-generation sequencing (NGS) plays an important role in relevant clinical research. The mutation frequency is a major problem in tumor mutation detection and increasing sequencing depth is a widely used method to improve mutation calling performance. Therefore, it is necessary to evaluate the effect of different sequencing depth and mutation frequency as well as mutation calling tools. In this study, Strelka2 and Mutect2 tools were used in detecting the performance of 30 combinations of sequencing depth and mutation frequency. Results showed that the precision rate kept greater than 95% in most of the samples. Generally, for higher mutation frequency (≥20%), sequencing depth ≥200X is sufficient for calling 95% mutations; for lower mutation frequency (≤10%), we recommend improving experimental method rather than increasing sequencing depth. Besides, according to our results, although Strelka2 and Mutect2 performed similarly, the former performed slightly better than the latter one at higher mutation frequency (≥20%), while Mutect2 performed better when the mutation frequency was lower than 10%. Besides, Strelka2 was 17 to 22 times faster than Mutect2 on average. Our research will provide a useful and comprehensive guideline for clinical genomic researches on somatic mutation identification through systematic performance comparison among different sequencing depths and mutation frequency.

Whole mitochondrial genome sequencing highlights mitochondrial impact in gastric cancer

Mitochondria are organelles that perform major roles in cellular operation. Thus, alterations in mitochondrial genome (mtGenome) may lead to mitochondrial dysfunction and cellular deregulation, influencing carcinogenesis. Gastric cancer (GC) is one of the most incident and mortal types of cancer in Brazil, particularly in the Amazon region. Here, we sequenced and compared the whole mtGenome extracted from FFPE tissue samples of GC patients (tumor and internal control – IC) and cancer-free individuals (external control – EC) from this region. We found 3-fold more variants and up to 9-fold more heteroplasmic regions in tumor when compared to paired IC samples. Moreover, tumor presented more heteroplasmic variants when compared to EC, while IC and EC showed no significant difference when compared to each other. Tumor also presented substantially more variants in the following regions: MT-RNR1, MT-ND5, MT-ND4, MT-ND2, MT-DLOOP1 and MT-CO1. In addition, our haplogroup results indicate an association of Native American ancestry (particularly haplogroup C) to gastric cancer development. To the best of our knowledge, this is the first study to sequence the whole mtGenome from FFPE samples and to apply mtGenome analysis in association to GC in Brazil.

pTrimmer: An efficient tool to trim primers of multiplex deep sequencing data

Background

With the widespread use of multiple amplicon-sequencing (MAS) in genetic variation detection, an efficient tool is required to remove primer sequences from short reads to ensure the reliability of downstream analysis. Although some tools are currently available, their efficiency and accuracy require improvement in trimming large scale of primers in high throughput target genome sequencing. This issue is becoming more urgent considering the potential clinical implementation of MAS for processing patient samples. We here developed pTrimmer that could handle thousands of primers simultaneously with greatly improved accuracy and performance.

Result

pTrimmer combines the two algorithms of k-mers and Needleman-Wunsch algorithm, which ensures its accuracy even with the presence of sequencing errors. pTrimmer has an improvement of 28.59% sensitivity and 11.87% accuracy compared to the similar tools. The simulation showed pTrimmer has an ultra-high sensitivity rate of 99.96% and accuracy of 97.38% compared to cutPrimers (70.85% sensitivity rate and 58.73% accuracy). And the performance of pTrimmer is notably higher. It is about 370 times faster than cutPrimers and even 17,000 times faster than cutadapt per threads. Trimming 2158 pairs of primers from 11 million reads (Illumina PE 150 bp) takes only 37 s and no more than 100 MB of memory consumption.

Conclusions

pTrimmer is designed to trim primer sequence from multiplex amplicon sequencing and target sequencing. It is highly sensitive and specific compared to other three similar tools, which could help users to get more reliable mutational information for downstream analysis.

Electronic supplementary material

The online version of this article (10.1186/s12859-019-2854-x) contains supplementary material, which is available to authorized users.

Fixation Effects on Variant Calling in a Clinical Resequencing Panel

Formalin fixation is the standard method for the preservation of tissue for diagnostic purposes, including pathologic review and molecular assays. However, this method is known to cause artifacts that can affect the accuracy of molecular genetic test results. We assessed the applicability of alternative fixatives to determine whether these perform significantly better on next-generation sequencing assays, and whether adequate morphology is retained for primary diagnosis, in a prospective study using a clinical-grade, laboratory-developed targeted resequencing assay. Several parameters relating to sequencing quality and variant calling were examined and quantified in tumor and normal colon epithelial tissues. We identified an alternative fixative that suppresses many formalin-related artifacts while retaining adequate morphology for pathologic review.

Using plasma cell-free DNA to monitor the chemoradiotherapy course of cervical cancer

The liquid biopsy is being integrated into cancer diagnostics and surveillance. However, critical questions still remain, such as how to precisely evaluate cancer mutation burden and interpret the corresponding clinical implications. Herein, we evaluated the role of peripheral blood cell-free DNA (cfDNA) in characterizing the dynamic mutation alterations of 48 cancer driver genes from cervical cancer patients. We performed targeted deep sequencing on 93 plasma cfDNA from 57 cervical cancer patients and from this developed an algorithm, allele fraction deviation (AFD), to monitor in an unbiased manner the dynamic changes of genomic aberrations. Differing treatments, including chemotherapy (n = 22), radiotherapy (n = 14) and surgery (n = 15), led to a significant decrease in AFD values (Wilcoxon, p = 0.029). The decrease of cfDNA AFD values was accompanied by shrinkage in the size of the tumor in most patients. However, in a subgroup of patients where cfDNA AFD values did not reflect a reduction in tumor size, there was a detection of progressive disease (metastasis). Furthermore, a low AFD value at diagnosis followed a later increase of AFD value also successfully predicted relapse. These results show that plasma cfDNA, together with targeted deep sequencing, may help predict treatment response and disease development in cervical cancer.

Clinical whole-genome sequencing from routine formalin-fixed, paraffin-embedded specimens: pilot study for the 100,000 Genomes Project

PURPOSE

Fresh-frozen (FF) tissue is the optimal source of DNA for whole-genome sequencing (WGS) of cancer patients. However, it is not always available, limiting the widespread application of WGS in clinical practice. We explored the viability of using formalin-fixed, paraffin-embedded (FFPE) tissues, available routinely for cancer patients, as a source of DNA for clinical WGS.

METHODS

We conducted a prospective study using DNAs from matched FF, FFPE, and peripheral blood germ-line specimens collected from 52 cancer patients (156 samples) following routine diagnostic protocols. We compared somatic variants detected in FFPE and matching FF samples.

RESULTS

We found the single-nucleotide variant agreement reached 71% across the genome and somatic copy-number alterations (CNAs) detection from FFPE samples was suboptimal (0.44 median correlation with FF) due to nonuniform coverage. CNA detection was improved significantly with lower reverse crosslinking temperature in FFPE DNA extraction (80 °C or 65 °C depending on the methods). Our final data showed somatic variant detection from FFPE for clinical decision making is possible. We detected 98% of clinically actionable variants (including 30/31 CNAs).

CONCLUSION

We present the first prospective WGS study of cancer patients using FFPE specimens collected in a routine clinical environment proving WGS can be applied in the clinic.

Describing the Reportable Range Is Important for Reliable Treatment Decisions: A Multiple Laboratory Study for Molecular Tumor Profiling Using Next-Generation Sequencing

Because interpretation of next-generation sequencing (NGS) data remains challenging, optimization of the NGS process is needed to obtain correct sequencing results. Therefore, extensive validation and continuous monitoring of the quality is essential. NGS performance was compared with traditional detection methods and technical quality of nine NGS technologies was assessed. First, nine formalin-fixed, paraffin-embedded patient samples were analyzed by 114 laboratories by using different detection methods. No significant differences in performance were observed between analyses with NGS and traditional techniques. Second, two DNA control samples were analyzed for a selected number of variants by 26 participants with the use of nine different NGS technologies. Quality control metrics were analyzed from raw data files and a survey about routine procedures. Results showed large differences in coverages, but observed variant allele frequencies in raw data files were in line with predefined variant allele frequencies. Many false negative results were found because of low-quality regions, which were not reported as such. It is recommended to disclose the reportable range, the fraction of targeted genomic regions for which calls of acceptable quality can be generated, to avoid any errors in therapy decisions. NGS can be a reliable technique, only if essential quality control during analysis is applied and reported.

Exploring DNA quality of single cells for genome analysis with simultaneous whole-genome amplification

Single cell genome analysis methods are powerful tools to define features of single cells and to identify differences between them. Since the DNA amount of a single cell is very limited, cellular DNA usually needs to be amplified by whole-genome amplification before being subjected to further analysis. A single nucleus only contains two haploid genomes. Thus, any DNA damage that prevents amplification results in loss of damaged DNA sites and induces an amplification bias. Therefore, the assessment of single cell DNA quality is urgently required. As of today, there is no simple method to determine the quality of a single cell DNA in a manner that will still retain the entire cellular DNA for amplification and downstream analysis. Here, we describe a method for whole-genome amplification with simultaneous quality control of single cell DNA by using a competitive spike-in DNA template.

Biomarkers in colorectal liver metastases

BACKGROUND

Despite a 5-year overall survival rate of 58 per cent after liver resection for colorectal liver metastases (CLMs), more than half of patients develop recurrence, highlighting the need for accurate risk stratification and prognostication. Traditional prognostic factors have been superseded by newer outcome predictors, including those defined by the molecular origin of the primary tumour.

METHODS

This review synthesized findings in the literature using the PubMed database of articles in the English language published between 1998 and 2017 on prognostic and predictive biomarkers in patients undergoing resection of CLMs.

RESULTS

Responses to preoperative chemotherapy define prognosis in patients undergoing CLM resection. There are differences by embryological origin too. Somatic mutations in the proto-oncogenes KRAS and NRAS are associated with positive surgical margins and tumour regrowth after ablation. Other mutations (such as BRAF) and co-occurring mutations in RAS/TP53 and APC/PIK3CA have emerged as important biomarkers that determine an individual patient's tumour biology and may be used to predict outcome after CLM resection.

CONCLUSION

Knowledge of somatic mutations can guide the use of preoperative therapy, extent of surgical margin and selection for ablation alone.

cutPrimers: A New Tool for Accurate Cutting of Primers from Reads of Targeted Next Generation Sequencing

Cutting of primers from reads is an important step of processing targeted amplicon-based next generation sequencing data. Existing tools are adapted for cutting of one or several primer/adapter sequences from reads and removing all of their occurrences. Also most of the existing tools use kmers and may cut only part of primers or primers with studied sequence of gene. Because of this, use of such programs leads to incorrect trimming, reduction of coverage, and increase in the number of false-positive and/or false-negative results. We have developed a new tool named cutPrimers for accurate cutting of any number of primers from reads. Using sequencing reads that were obtained during study of BRCA1/2 genes, we compared it with cutadapt, AlienTrimmer, and BBDuk. All of them trimmed reads in such a way that coverage of at least two amplicons decreased to unacceptable level (<30 reads) comparing with reads trimmed with cutPrimers. At the same time, Trimmomatic and AlienTrimmer cut all occurrences of primer sequences, so the length of the remaining reads was less than prospective.

Non-reproducible sequence artifacts in FFPE tissue: an experience report

BACKGROUND

Recent advances in sequencing technologies supported the development of molecularly targeted therapy in cancer patients. Thus, genomic analyses are becoming a routine part in clinical practice and accurate detection of actionable mutations is essential to assist diagnosis and therapy choice. However, this is often challenging due to major problems associated with DNA from formalin-fixed paraffin-embedded tissue which is usually the primary source for genetic testing.

OBJECTIVES

Here we want to share our experience regarding major problems associated with FFPE DNA used for PCR-based sequencing as illustrated by the mutational analysis of ERBB4 in melanoma. We want to focus on two major problems including extensive DNA fragmentation and hydrolytic deamination as source of non-reproducible sequence artifacts. Further, we provide potential explanations and possible strategies to minimize these difficulties and improve the detection of targetable mutations.

METHODS

Genomic DNA from formalin-fixed paraffin-embedded tumor samples was isolated followed by PCR amplification, Sanger sequencing and statistical analysis.

RESULTS

Analysis of Sanger sequencing data revealed a total of 46 ERBB4 mutations in 27 of 96 samples including the identification of 11 mutations at three previously unknown mutational hotspots. Unfortunately, we were not able to confirm any assumed hotspot mutation within repeated sequencing of relevant amplicons suggesting the detection of sequence artifacts most likely caused by DNA lesions associated with FFPE tissues.

CONCLUSION

Since DNA from FFPE tissue is usually the primary source for mutational analyses, appropriate measures must be implemented in the workflow to assess DNA damage in formalin-fixed tissue to ensure accurate detection of actionable mutations and minimize the occurrence of sequence artifacts.

Use of the QIAGEN GeneReader NGS system for detection of KRAS mutations, validated by the QIAGEN Therascreen PCR kit and alternative NGS platform

BACKGROUND

The detection of somatic mutations in primary tumors is critical for the understanding of cancer evolution and targeting therapy. Multiple technologies have been developed to enable the detection of such mutations. Next generation sequencing (NGS) is a new platform that is gradually becoming the technology of choice for genotyping cancer samples, owing to its ability to simultaneously interrogate many genomic loci at massively high efficiency and increasingly lower cost. However, multiple barriers still exist for its broader adoption in clinical research practice, such as fragmented workflow and complex bioinformatics analysis and interpretation.

METHODS

We performed validation of the QIAGEN GeneReader NGS System using the QIAact Actionable Insights Tumor Panel, focusing on clinically meaningful mutations by using DNA extracted from formalin-fixed paraffin-embedded (FFPE) colorectal tissue with known KRAS mutations. The performance of the GeneReader was evaluated and compared to data generated from alternative technologies (PCR and pyrosequencing) as well as an alternative NGS platform. The results were further confirmed with Sanger sequencing.

RESULTS

The data generated from the GeneReader achieved 100% concordance with reference technologies. Furthermore, the GeneReader workflow provides a truly integrated workflow, eliminating artifacts resulting from routine sample preparation; and providing up-to-date interpretation of test results.

CONCLUSION

The GeneReader NGS system offers an effective and efficient method to identify somatic (KRAS) cancer mutations.

Assessment of the quality of DNA from various formalin-fixed paraffin-embedded (FFPE) tissues and the use of this DNA for next-generation sequencing (NGS) with no artifactual mutation

Formalin-fixed, paraffin-embedded (FFPE) tissues used for pathological diagnosis are valuable for studying cancer genomics. In particular, laser-capture microdissection of target cells determined by histopathology combined with FFPE tissue section immunohistochemistry (IHC) enables precise analysis by next-generation sequencing (NGS) of the genetic events occurring in cancer. The result is a new strategy for a pathological tool for cancer diagnosis: 'microgenomics'. To more conveniently and precisely perform microgenomics, we revealed by systematic analysis the following three details regarding FFPE DNA compared with paired frozen tissue DNA. 1) The best quality of FFPE DNA is obtained by tissue fixation with 10% neutral buffered formalin for 1 day and heat treatment of tissue lysates at 95°C for 30 minutes. 2) IHC staining of FFPE tissues decreases the quantity and quality of FFPE DNA to one-fourth, and antigen retrieval (at 120°C for 15 minutes, pH 6.0) is the major reason for this decrease. 3) FFPE DNA prepared as described herein is sufficient for NGS. For non-mutated tissue specimens, no artifactual mutation occurs during FFPE preparation, as shown by precise comparison of NGS of FFPE DNA and paired frozen tissue DNA followed by validation. These results demonstrate that even FFPE tissues used for routine clinical diagnosis can be utilized to obtain reliable NGS data if appropriate conditions of fixation and validation are applied.

Analysis of PIK3CA mutations and PI3K pathway proteins in advanced gastric cancer

BACKGROUND

Although surgery and chemotherapy have extended advanced gastric cancer patient survival, some patients still experience relapse and metastasis. We postulated that PI3K pathway proteins could be prognostic biomarkers for the advanced gastric cancer patients.

METHODS

A retrospective cohort of 160 advanced gastric cancer patients receiving potentially curative surgery with/without chemotherapy was investigated for PIK3CA mutation and PI3K pathway protein level in the context of overall survival and relapse-free survival.

RESULTS

Thirteen patients (13 of 111, 11.7%) had PIK3CA mutations in codon 545, whereas one patient (1 of 94, 1.1%) had a mutation in PIK3CA codon 1047. PI3K pathway protein immunohistochemistry demonstrated that phosphorylated AKT positive [p-AKT (+)] patients in the surgery-only group had a good prognosis in terms of overall survival and relapse-free survival. No significant association between PIK3CA mutations and PI3K pathway protein level was seen.

CONCLUSIONS

This study revealed that (1) PIK3CA hotspot mutations occurred with low frequency in gastric cancer; (2) PIK3CA hotspot mutations were not directly associated with PI3K pathway activation; and (3) p-AKT (+) may be a biomarker for better outcomes for gastric cancer patients undergoing gastrectomy regardless of the PIK3CA mutation status.

Spatio-temporal mutation profiles of case-matched colorectal carcinomas and their metastases reveal unique de novo mutations in metachronous lung metastases by targeted next generation sequencing

Background

Targeted next generation sequencing (tNGS) has become part of molecular pathology diagnostics for determining RAS mutation status in colorectal cancer (CRC) patients as predictive tool for decision on EGFR-targeted therapy. Here, we investigated mutation profiles of case-matched tissue specimens throughout the disease course of CRC, to further specify RAS-status dynamics and to identify de novo mutations associated with distant metastases.

Methods

Case-matched formalin-fixed and paraffin-embedded (FFPE) resection specimens (n = 70; primary tumours, synchronous and/or metachronous liver and/or lung metastases) of 14 CRC cases were subjected to microdissection of normal colonic epithelial, primary and metastatic tumour cells, their DNA extraction and an adapted library protocol for limited DNA using the 48 gene TruSeq Amplicon Cancer Panel^TM, MiSeq sequencing and data analyses (Illumina).

Results

By tNGS primary tumours were RAS wildtype in 5/14 and mutated in 9/14 (8/9 KRAS exon 2; 1/9 NRAS Exon 3) of cases. RAS mutation status was maintained in case-matched metastases throughout the disease course, albeit with altered allele frequencies. Case-matched analyses further identified a maximum of three sequence variants (mainly in APC, KRAS, NRAS, TP53) shared by all tumour specimens throughout the disease course per individual case. In addition, further case-matched de novo mutations were detected in synchronous and/or metachronous liver and/or lung metastases (e.g. in APC, ATM, FBXW7, FGFR3, GNAQ, KIT, PIK3CA, PTEN, SMAD4, SMO, STK11, TP53, VHL). Moreover, several de novo mutations were more frequent in synchronous (e.g. ATM, KIT, PIK3CA, SMAD4) or metachronous (e.g. FBXW7, SMO, STK11) lung metastases. Finally, some de novo mutations occurred only in metachronous lung metastases (CDKN2A, FGFR2, GNAS, JAK3, SRC).

Conclusion

Together, this study employs an adapted FFPE-based tNGS approach to confirm conservation of RAS mutation status in primary and metastatic tissue specimens of CRC patients. Moreover, it identifies genes preferentially mutated de novo in late disease stages of metachronous CRC lung metastases, several of which might be actionable by targeted therapies.

Electronic supplementary material

The online version of this article (doi:10.1186/s12943-016-0549-8) contains supplementary material, which is available to authorized users.

A rational two-step approach to KRAS mutation testing in colorectal cancer using high resolution melting analysis and pyrosequencing

Background

KRAS mutation testing is mandatory in the management of metastatic colorectal cancer prior to treatment with anti-EGFR antibodies as patients whose tumors express mutant KRAS do not benefit from these agents. Although the U.S. Food and Drug Administration has recently approved two in-vitro diagnostics kits for determination of KRAS status, there is generally no consensus on the preferred method and new tests are continuously being developed. Most of these techniques focus on the hotspot mutations at codons 12 and 13 of the KRAS gene.

Methods

We describe a two-step approach to KRAS codon 12/13 mutation testing involving high resolution melting analysis (HRM) followed by pyrosequencing using the Therascreen KRAS Pyro kit (Qiagen) of only those samples that are not clearly identified as KRAS wildtype or mutant by HRM. First, we determined KRAS status in a panel of 61 colorectal cancer samples using both methods to compare technical performance and concordance of results. Subsequently, we evaluated practicability and costs of our concept in an independent set of 120 colorectal cancer samples in a routine diagnostic setting.

Results

HRM and pyrosequencing appeared to be equally sensitive, allowing for clear detection of mutant alleles at a mutant allele frequency ≥12.5 %. Pyrosequencing yielded more exploitable results due to lower input requirements and a lower rate of analysis failures. KRAS codon 12/13 status was called concordantly for 98.2 % (56/57) of all samples that could be successfully analysed by both methods and 100 % (19/19) of samples that were identified mutant by HRM. Reviewing the actual effort and expenses for KRAS mutation testing in our laboratory revealed, that the selective use of pyrosequencing for only those samples that could not be analysed by HRM increased the fraction of valid results from 87.5 % for HRM alone to 99.2 % (119/120) while allowing for a net reduction of operational costs of >75 % compared to pyrosequencing alone.

Conclusions

Combination of HRM and pyrosequencing in a two-step diagnostic procedure constitutes a reliable and economic analysis platform for KRAS mutation testing in colorectal cancer in a clinical setting.

Electronic supplementary material

The online version of this article (doi:10.1186/s12885-016-2589-2) contains supplementary material, which is available to authorized users.

Assessment of clinically related outcomes and biomarker analysis for translational integration in colorectal cancer (ACROBATICC): study protocol for a population-based, consecutive cohort of surgically treated colorectal cancers and resected colorectal liver metastasis

Background

More accurate predictive and prognostic biomarkers for patients with colorectal cancer (CRC) primaries or colorectal liver metastasis (CLM) are needed. Outside clinical trials, the translational integration of emerging pathways and novel techniques should facilitate exploration of biomarkers for improved staging and prognosis.

Methods

An observational study exploring predictive and prognostic biomarkers in a population-based, consecutive cohort of surgically treated colorectal cancers and resected colorectal liver metastases. Long-term outcomes will be cancer-specific survival, recurrence-free survival and overall survival at 5 years from diagnosis. Beyond routine clinicopathological and anthropometric characteristics and laboratory and biochemistry results, the project allows for additional blood samples and fresh-frozen tumour and normal tissue for investigation of circulating tumour cells (CTCs) and novel biomarkers (e.g. immune cells, microRNAs etc.). Tumour specimens will be investigated by immunohistochemistry in full slides. Extracted DNA/RNA will be analysed for genomic markers using specific PCR techniques and next-generation sequencing (NGS) panels. Flow cytometry will be used to characterise biomarkers in blood. Collaboration is open and welcomed, with particular interest in mutual opportunities for validation studies.

Status and perspectives

The project is ongoing and recruiting at an expected rate of 120–150 patients per year, since January 2013. A project on circulating tumour cells (CTCs) has commenced, with analysis being prepared. Investigating molecular classes beyond the TNM staging is under way, including characteristics of microsatellite instability (MSI) and elevated microsatellite alterations in selected tetranucleotides (EMAST). Hot spot panels for known mutations in CRC are being investigated using NGS. Immune-cell characteristics are being performed by IHC and flow cytometry in tumour and peripheral blood samples. The project has ethical approval (REK Helse Vest, #2012/742), is financially supported with a Ph.D.-Grant (EMAST project; Folke Hermansen Cancer Fund) and a CTC-project (Norwegian Research Council; O. Nordgård). The ACROBATICC clinical and molecular biobank repository will serve as a long-term source for novel exploratory analysis and invite collaborators for mutual validation of promising biomarker results. The project aims to generate results that can help better discern prognostic groups in stage II/III cancers; explore prognostic and predictive biomarkers, and help detail the biology of colorectal liver metastasis for better patient selection and tailored treatment. The project is registered at http://www.ClinicalTrials.gov NCT01762813.

Advances in clinical next-generation sequencing: target enrichment and sequencing technologies

The huge parallel sequencing capabilities of next generation sequencing technologies have made them the tools of choice to characterize genomic aberrations for research and diagnostic purposes. For clinical applications, screening the whole genome or exome is challenging owing to the large genomic area to be sequenced, associated costs, complexity of data, and lack of known clinical significance of all genes. Consequently, routine screening involves limited markers with established clinical relevance. This process, referred to as targeted genome sequencing, requires selective enrichment of the genomic areas comprising these markers via one of several primer or probe-based enrichment strategies, followed by sequencing of the enriched genomic areas. Here, the authors review current target enrichment approaches and next generation sequencing platforms, focusing on the underlying principles, capabilities, and limitations of each technology along with validation and implementation for clinical testing.