Targeted-capture massively-parallel sequencing enables robust detection of clinically informative mutations from formalin-fixed tumours

Large-scale analysis of whole genome sequencing data from formalin-fixed paraffin-embedded cancer specimens demonstrates preservation of clinical utility

Whole genome sequencing (WGS) provides comprehensive, individualised cancer genomic information. However, routine tumour biopsies are formalin-fixed and paraffin-embedded (FFPE), damaging DNA, historically limiting their use in WGS. Here we analyse FFPE cancer WGS datasets from England's 100,000 Genomes Project, comparing 578 FFPE samples with 11,014 fresh frozen (FF) samples across multiple tumour types. We use an approach that characterises rather than discards artefacts. We identify three artefactual signatures, including one known (SBS57) and two previously uncharacterised (SBS FFPE, ID FFPE), and develop an "FFPEImpact" score that quantifies sample artefacts. Despite inferior sequencing quality, FFPE-derived data identifies clinically-actionable variants, mutational signatures and permits algorithmic stratification. Matched FF/FFPE validation cohorts shows good concordance while acknowledging SBS, ID and copy-number artefacts. While FF-derived WGS data remains the gold standard, FFPE-samples can be used for WGS if required, using analytical advancements developed here, potentially democratising whole cancer genomics to many.

Histologic and Genomic Analysis of Conjunctival SCC in African and American Cohorts Reveal UV Light and HPV Signatures and High Tumor Mutation Burden

Purpose

Conjunctival squamous cell carcinoma (conjSCC) is more prevalent and aggressive in sub-Saharan African countries compared with the rest of the world. This study aims to compare the genomic, immunophenotypic, and histologic features between patients from the United States and Ethiopia, to identify etiopathogenic mechanisms and unveil potential treatment strategies.

Methods

We compared histologic features and mutational profiles using whole exome sequencing, high-risk human papillomavirus (HPV) status, PD-L1 expression, and tumor-infiltrating lymphocytes in conjSCC tumors of patients from Ethiopia (ETH; n = 25) and the United States (from MD Anderson [the MDA cohort]; n = 29). Genomic alterations were compared with SCCs from other anatomic sites using data from The Cancer Genome Atlas.

Results

Solar elastosis was seen in 78% of ETH and 10% of MDA samples. Thicker tumors had higher density of CD8+ and CD3+ cells. HPV status was similar between the cohorts (ETH = 21% and MDA = 28%). The mean tumor mutation burden (TMB) was significantly higher in conjSCC (3.01/Mb, log10) and cutaneous SCC compared other SCC subtypes. ETH samples had higher TMB compared to the MDA cohort (3.34 vs. 2.73). Mutations in genes associated with ultraviolet light (UV) signature were most frequently encountered (SBS7b = 74% and SBS7a = 72%), with higher prevalence in the ETH cohort, whereas SBS2 and SBS13 signatures were more common among MDA HPV+ conjSCCs.

Conclusions

Our findings suggest that UV exposure may play a major role in conjSCC, with a higher prevalence in the ETH cohort compared with the MDA cohort, where HPV also contributes.

Characterization and mitigation of artifacts derived from NGS library preparation due to structure-specific sequences in the human genome

BACKGROUND

Hybridization capture-based targeted next generation sequencing (NGS) is gaining importance in routine cancer clinical practice. DNA library preparation is a fundamental step to produce high-quality sequencing data. Numerous unexpected, low variant allele frequency calls were observed in libraries using sonication fragmentation and enzymatic fragmentation. In this study, we investigated the characteristics of the artifact reads induced by sonication and enzymatic fragmentation. We also developed a bioinformatic algorithm to filter these sequencing errors.

RESULTS

We used pairwise comparisons of somatic single nucleotide variants (SNVs) and insertions and deletions (indels) of the same tumor DNA samples prepared using both ultrasonic and enzymatic fragmentation protocols. Our analysis revealed that the number of artifact variants was significantly greater in the samples generated using enzymatic fragmentation than using sonication. Most of the artifacts derived from the sonication-treated libraries were chimeric artifact reads containing both cis- and trans-inverted repeat sequences of the genomic DNA. In contrast, chimeric artifact reads of endonuclease-treated libraries contained palindromic sequences with mismatched bases. Based on these distinctive features, we proposed a mechanistic hypothesis model, PDSM (pairing of partial single strands derived from a similar molecule), by which these sequencing errors derive from ultrasonication and enzymatic fragmentation library preparation. We developed a bioinformatic algorithm to generate a custom mutation "blacklist" in the BED region to reduce errors in downstream analyses.

CONCLUSIONS

We first proposed a mechanistic hypothesis model (PDSM) of sequencing errors caused by specific structures of inverted repeat sequences and palindromic sequences in the natural genome. This new hypothesis predicts the existence of chimeric reads that could not be explained by previous models, and provides a new direction for further improving NGS analysis accuracy. A bioinformatic algorithm, ArtifactsFinder, was developed and used to reduce the sequencing errors in libraries produced using sonication and enzymatic fragmentation.

MET Amplification as a Resistance Driver to TKI Therapies in Lung Cancer: Clinical Challenges and Opportunities

Targeted therapy has emerged as an important pillar for the standard of care in oncogene-driven non-small cell lung cancer (NSCLC), which significantly improved outcomes of patients whose tumors harbor oncogenic driver mutations. However, tumors eventually develop resistance to targeted drugs, and mechanisms of resistance can be diverse. MET amplification has been proven to be a driver of resistance to tyrosine kinase inhibitor (TKI)-treated advanced NSCLC with its activation of EGFR, ALK, RET, and ROS-1 alterations. The combined therapy of MET-TKIs and EGFR-TKIs has shown outstanding clinical efficacy in EGFR-mutated NSCLC with secondary MET amplification-mediated resistance in a series of clinical trials. In this review, we aimed to clarify the underlying mechanisms of MET amplification-mediated resistance to tyrosine kinase inhibitors, discuss the ways and challenges in the detection and diagnosis of MET amplifications in patients with metastatic NSCLC, and summarize the recently published clinical data as well as ongoing trials of new combination strategies to overcome MET amplification-mediated TKI resistance.

One-Step Synthesis of Nanoliposomal Copper Diethyldithiocarbamate and Its Assessment for Cancer Therapy

The metal complex copper diethyldithiocarbamate (CuET) induces cancer cell death by inhibiting protein degradation and induces proteotoxic stress, making CuET a promising cancer therapeutic. However, no clinical formulation of CuET exists to date as the drug is insoluble in water and exhibits poor bioavailability. To develop a scalable formulation, nanoliposomal (LP) CuET was synthesized using ethanol injection as a facile one-step method that is suitable for large-scale manufacturing. The nanoparticles are monodispersed, colloidally stable, and approximately 100 nm in diameter with an encapsulation efficiency of over 80%. LP-CuET demonstrates excellent stability in plasma, minimal size change, and little drug release after six-month storage at various temperatures. Additionally, melanoma cell lines exhibit significant sensitivity to LP-CuET and cellular uptake occurs predominantly through endocytosis in YUMM 1.7 cancer cells. Intracellular drug delivery is mediated by vesicle acidification with more nanoparticles being internalized by melanoma cells compared with RAW 264.7 macrophages. Additionally, the nanoparticles preferentially accumulate in YUMM 1.7 tumors where they induce cancer cell death in vivo. The development and characterization of a stable and scalable CuET formulation illustrated in this study fulfils the requirements needed for a potent clinical grade formulation.

Characterization of DNA lesions associated with cell-free DNA by targeted deep sequencing

BACKGROUND

Recently, a next-generation sequencing (NGS)-based method has been used for the successful detection of circulating tumor DNA (ctDNA) in various cancer types. Thus, the use of NGS on liquid biopsies will improve cancer diagnosis and prognosis. However, the low-allelic fraction of ctDNA poses a challenge for the sensitive and specific detection of tumor variants in cell-free DNA (cfDNA). To distinguish true variants from false positives, the characteristics of errors that occur during sample preparation and sequencing need to be elucidated.

METHODS

We generated capture-based targeted deep sequencing data from plasma cfDNA and peripheral blood leucocyte (PBL) gDNA to profile background errors. To reveal cfDNA-associated DNA lesions, background error profiles from two sample types were compared in each nucleotide substitution class.

RESULTS

In this study, we determined the prevalence of single nucleotide substitutions in cfDNA sequencing data to identify DNA damage preferentially associated with cfDNA. On comparing sequencing errors between cfDNA and cellular genomic DNA (gDNA), we observed that the total substitution error rates in cfDNA were significantly higher than those in gDNA. When the substitution errors were divided into 12 substitution error classes, C:G>T:A substitution errors constituted the largest difference between cfDNA and gDNA samples. When the substitution error rates were estimated based on the location of DNA-fragment substitutions, the differences in error rates of most substitution classes between cfDNA and gDNA samples were observed only at the ends of the DNA fragments. In contrast, C:G>T:A substitution errors in the cfDNA samples were not particularly associated with DNA-fragment ends. All observations were verified in an independent dataset.

CONCLUSIONS

Our data suggested that cytosine deamination increased in cfDNA compared to that in cellular gDNA. Such an observation might be due to the attenuation of DNA damage repair before the release of cfDNA and/or the accumulation of cytosine deamination after it. These findings can contribute to a better understanding of cfDNA-associated DNA damage, which will enable the accurate analysis of somatic variants present in cfDNA at an extremely low frequency.

Evolution of late-stage metastatic melanoma is dominated by aneuploidy and whole genome doubling

Although melanoma is initiated by acquisition of point mutations and limited focal copy number alterations in melanocytes-of-origin, the nature of genetic changes that characterise lethal metastatic disease is poorly understood. Here, we analyze the evolution of human melanoma progressing from early to late disease in 13 patients by sampling their tumours at multiple sites and times. Whole exome and genome sequencing data from 88 tumour samples reveals only limited gain of point mutations generally, with net mutational loss in some metastases. In contrast, melanoma evolution is dominated by whole genome doubling and large-scale aneuploidy, in which widespread loss of heterozygosity sculpts the burden of point mutations, neoantigens and structural variants even in treatment-naïve and primary cutaneous melanomas in some patients. These results imply that dysregulation of genomic integrity is a key driver of selective clonal advantage during melanoma progression.

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:

A Pilot Study for the Feasibility of Exome-Sequencing in Circulating Tumor Cells Versus Single Metastatic Biopsies in Breast Cancer

The comparison of the landscape of somatic alterations in circulating tumor cells (CTCs) versus metastases is challenging. Here, we comprehensively characterized the somatic landscape in bulk (amplified and non-amplified), spike-in breast cancer cells, CTCs, and metastases from breast cancer patients using whole-exome sequencing (WES). We determined the level of genomic concordance for somatic nucleotide variants (SNVs), copy number alterations (CNAs), and structural variants (SVs). The variant allele fractions (VAFs) of somatic variants were remarkably similar between amplified and non-amplified cell line samples as technical replicates. In clinical samples, a significant fraction of somatic variants had low VAFs in CTCs compared to metastases. The most frequently recurrent gene mutations in clinical samples were associated with an elevated C > T mutational signature. We found complex rearrangement patterns including intra- and inter-chromosomal rearrangements, singleton, and recurrent gene fusions, and tandem duplications. We observed high molecular discordance for somatic alterations between paired samples consistent with marked heterogeneity of the somatic landscape. The most prevalent copy number calls were focal deletion events in CTCs and metastases. Our results demonstrate the feasibility of an integrated workflow for the identification of a complete repertoire of somatic alterations and highlight the intrapatient genomic differences that occur between CTCs and metastases.

Using postmortem formalin fixed paraffin-embedded tissues for molecular testing of sudden cardiac death: A cautionary tale of utility and limitations

For archived cases of previously young healthy individuals where cause of sudden death remains undetermined, formalin fixed paraffin-embedded tissues (FFPE) samples are often the only biological resource available for molecular testing. We aim to ascertain the validity of postmortem molecular analysis of 95 cardiac genes using the FFPE samples routinely processed in the offices of medical examiners - typical fixation time in formalin ranges from days to months. The study was conducted in the College of American Pathologists accredited Molecular Genetics Laboratory within the City of New York Office of Chief Medical Examiner. Twelve cases, with FFPE samples and corresponding non-formalin fixed samples (RNAlater-preserved tissues or bloodstain card), were chosen for testing results comparison. The methods of extracting DNA from FFPE samples using Covaris, Qiagen, and Promega products showed comparable results. The quality of the extracted DNA, the target-enriched DNA libraries of 95 cardiac genes using HaloPlex Target Enrichment system by Agilent Technologies, and sequencing results using Illumina Miseq instrument were evaluated. Compared to the sequencing results of the nonfixed samples, the FFPE samples were categorized into three groups: 1) Group 1 samples fixed in formalin 2-6 days, had greater than 55 % sequencing regions ≥30x and 94%-100% variant concordance. 2) Group 2 samples fixed in formalin for 8 days, showed intra-sample sequencing variations: the surface tissues showed 25%-27% extra variants (false positive) and 8.1%-9.7% missing variants (false negative), whereas the repeated core tissues showed reduced extra variants to 1.6 % and the false negative error was unchanged. 3) Group 3 samples fixed in formalin 29-136 days, had 2-55 % sequencing regions ≥30x, up to 52.2 % missed variants and up to 6.3 % extra variants. All reportable variants (pathogenic, likely pathogenic or variant of uncertain significance) identified in the nonfixed samples were also identified in FFPE, albeit three variants had low confidence variant calling. In summary, our study showed that postmortem molecular diagnostic testing using FFPE samples routinely processed by the medical examiners should be cautioned, as they are replete with false positive and negative results, particularly when sample fixation time is longer than 8 days. Saving non-formalin fixed samples for high fidelity molecular analysis is strongly encouraged.

Using ddPCR to assess the DNA yield of FFPE samples

Objectives

Detection of genomic alterations in diseases can be achieved with current molecular technologies. However, the molecules extracted from formalin-fixed, paraffin-embedded (FFPE) bio-samples are often limited possibly due to DNA fragmentation and crosslinking caused by the sample fixation and processing. The study objective was to design a droplet digital PCR (ddPCR) assay to assess the quality and quantity of DNA derived from various DNA extraction conditions on FFPE samples.

Methods

We used 10 μm-thick sections from 5 FFPE oral tumoral blocks, each consisting of 10–15 sections. The protocol variables tested included: 1) tissue staining; 2) duration and 3) temperature of post-digestion heat treatment; and 4) DNA extraction method. DNA quantity was assessed using the NanoDrop 2000 (Thermo Fisher Scientific, USA), the Qubit fluorometer (Thermo Fisher Scientific, USA), and a ddPCR-based assay. DNA quality was assessed using a ddPCR assay for the degree of fragmentation and the effectiveness of removing crosslinks with varying guanine-cytosine (GC)-content.

Results

Deparaffinization with xylene helped to increase the DNA yield. Tissue staining (methyl green staining, pH 6) prior to microdissection, comparing to no staining, caused additional DNA fragmentation. Compared to column-based method, DNA extracted with phenol chloroform and ethanol precipitation increased the degree of fragmentation and lowered the yield of amplifiable DNA. The cross-linking derived from GC-contents may not be the only factor impacting on the DNA quality.

Conclusions

Samples undergoing different pre-treatment conditions prior to extraction can impact the yield of amplifiable DNA. Our ddPCR assay can be used to assess for both DNA quantity and quality.

Assessment of concordance between fresh-frozen and formalin-fixed paraffin embedded tumor DNA methylation using a targeted sequencing approach

DNA methylation is altered in many types of disease, including metastatic colorectal cancer. However, the methylome has not yet been fully described in archival formalin-fixed paraffin embedded (FFPE) samples in the context of matched fresh-frozen (FF) tumor material at base-pair resolution using a targeted approach. Using next-generation sequencing, we investigated three pairs of matched FFPE and FF samples to determine the extent of their similarity. We identified a ‘bowing’ pattern specific to FFPE samples categorized by a lower CG proportion at the start of sequence reads. We have found no evidence that this affected methylation calling, nor concordance of results. We also found no significant increase in deamination, measured by C>T transitions, previously considered a result of crosslinking DNA by formalin fixation and a barrier to the use of FFPE in methylation studies. The methods used in this study have shown sensitivity of between 60-70% based on positions also methylated in colorectal cancer cell lines. We demonstrate that FFPE material is a useful source of tumor material for methylation studies using targeted sequencing.

Simple and Rapid Method to Obtain High-quality Tumor DNA from Clinical-pathological Specimens Using Touch Imprint Cytology

It is critical to determine the mutational status in cancer before administration and treatment of specific molecular targeted drugs for cancer patients. In the clinical setting, formalin-fixed paraffin-embedded (FFPE) tissues are widely used for genetic testing. However, FFPE DNA is generally damaged and fragmented during the fixation process with formalin. Therefore, FFPE DNA is sometimes not adequate for genetic testing because of low quality and quantity of DNA. Here we present a method of touch imprint cytology (TIC) to obtain genomic DNA from cancer cells, which can be observed under a microscope. Cell morphology and cancer cell numbers can be evaluated using TIC specimens. Furthermore, the extraction of genomic DNA from TIC samples can be completed within two days. The total amount and quality of TIC DNA obtained using this method was higher than that of FFPE DNA. This rapid and simple method allows researchers to obtain high-quality DNA for genetic testing (e.g., next generation sequencing analysis, digital PCR, and quantitative real time PCR) and to shorten the turnaround time for reporting results.

The reliable assurance of detecting somatic mutations in cancer-related genes by next-generation sequencing: the results of external quality assessment in China

To evaluate the proficiencies of laboratories utilizing next-generation sequencing (NGS) to detect somatic mutations in cancer-related genes, an external quality assessment (EQA) was implemented by the National Center for Clinical Laboratories of China in 2015. We prepared a panel of samples that comprised eight samples made by mixing synthetic mutated DNA fragments with normal human genomic DNA and one reference sample containing only genomic DNA. We validated our sample panel, and then distributed it to laboratories across China. We received complete results from 64 laboratories. The performances of 51.6 % (33/64) respondent labs were acceptable and 26.6 % (17/64) of the labs returned perfect results. In total, 449 mistakes were reported, including 201 false-negatives (201/449, 44.8 %) and 222 false-positives (222/449, 49.4 %) and 26 slightly discordant results (26/449, 5.8 %). We believe these unsatisfactory results and varied performances are mainly due to the enrichment methods used, the diverse sequencing chemistries of the different NGS platforms, and other errors within the sequencing process. The results indicate that our sample panel is suitable for use in EQA studies, and that further laboratory training in targeted NGS testing is urgently required. To address this, we propose a targeted NGS workflow with details on quality assurance procedures according to the current guidelines.

Assessment of Capture and Amplicon-Based Approaches for the Development of a Targeted Next-Generation Sequencing Pipeline to Personalize Lymphoma Management

Targeted next-generation sequencing panels are increasingly used to assess the value of gene mutations for clinical diagnostic purposes. For assay development, amplicon-based methods have been preferentially used on the basis of short preparation time and small DNA input amounts. However, capture sequencing has emerged as an alternative approach because of high testing accuracy. We compared capture hybridization and amplicon sequencing approaches using fresh-frozen and formalin-fixed, paraffin-embedded tumor samples from eight lymphoma patients. Next, we developed a targeted sequencing pipeline using a 32-gene panel for accurate detection of actionable mutations in formalin-fixed, paraffin-embedded tumor samples of the most common lymphocytic malignancies: chronic lymphocytic leukemia, diffuse large B-cell lymphoma, and follicular lymphoma. We show that hybrid capture is superior to amplicon sequencing by providing deep more uniform coverage and yielding higher sensitivity for variant calling. Sanger sequencing of 588 variants identified specificity limits of thresholds for mutation calling, and orthogonal validation on 66 cases indicated 93% concordance with whole-genome sequencing. The developed pipeline and assay identified at least one actionable mutation in 91% of tumors from 219 lymphoma patients and revealed subtype-specific mutation patterns and frequencies consistent with the literature. This pipeline is an accurate and sensitive method for identifying actionable gene mutations in routinely acquired biopsy materials, suggesting further assessment of capture-based assays in the context of personalized lymphoma management.

Application of Targeted Next-generation Sequencing, TruSeq Custom Amplicon Assay for Molecular Pathology Diagnostics on Formalin-fixed and Paraffin-embedded Samples

The implementation of targeted therapies revolutionized oncology. As the number of new oncogenic driver mutations, which provide molecular targets for prediction of effective and selective therapies, is increasing, the implementation of fast and reliable methods by molecular pathology labs is very important. Here we report our results with TruSeq Custom Amplicon assay performed on formalin-fixed and paraffin-embedded material. The oligo capture probes targeted the hotspot regions of 10 well-known oncogenes linked to clinical diagnosis and treatment of lung and colorectal adenocarcinomas, melanomas, and gastrointestinal stromal tumors. Fifteen previously genotyped formalin-fixed and paraffin-embedded DNA samples from different tumor types were selected for massively parallel sequencing. A bioinformatics pipeline was developed to identify high-quality variants and remove sequence artifacts. With the exception of 1 sample, which was of lower quality than the others, relevant mutations corresponding to tumor types could be reliable detected by the developed bioinformatical pipeline. This study indicates that the application of TruSeq Custom Amplicon assay is a promising tool in molecular pathology diagnostics, but it is important to standardize sample processing (including fixation, isolation procedure, sample selection based on quality assessment, and rigorous variant calling) to achieve the highest success rate and avoid false results.

Characterization of background noise in capture-based targeted sequencing data

Background

Targeted deep sequencing is increasingly used to detect low-allelic fraction variants; it is therefore essential that errors that constitute baseline noise and impose a practical limit on detection are characterized. In the present study, we systematically evaluate the extent to which errors are incurred during specific steps of the capture-based targeted sequencing process.

Results

We removed most sequencing artifacts by filtering out low-quality bases and then analyze the remaining background noise. By recognizing that plasma DNA is naturally fragmented to be of a size comparable to that of mono-nucleosomal DNA, we were able to identify and characterize errors that are specifically associated with acoustic shearing. Two-thirds of C:G > A:T errors and one quarter of C:G > G:C errors were attributed to the oxidation of guanine during acoustic shearing, and this was further validated by comparative experiments conducted under different shearing conditions. The acoustic shearing step also causes A > G and A > T substitutions localized to the end bases of sheared DNA fragments, indicating a probable association of these errors with DNA breakage. Finally, the hybrid selection step contributes to one-third of the remaining C:G > A:T and one-fifth of the C > T errors.

Conclusions

The results of this study provide a comprehensive summary of various errors incurred during targeted deep sequencing, and their underlying causes. This information will be invaluable to drive technical improvements in this sequencing method, and may increase the future usage of targeted deep sequencing methods for low-allelic fraction variant detection.

Electronic supplementary material

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

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.

The Utilization of Formalin Fixed-Paraffin-Embedded Specimens in High Throughput Genomic Studies

High throughput genomic assays empower us to study the entire human genome in short time with reasonable cost. Formalin fixed-paraffin-embedded (FFPE) tissue processing remains the most economical approach for longitudinal tissue specimen storage. Therefore, the ability to apply high throughput genomic applications to FFPE specimens can expand clinical assays and discovery. Many studies have measured the accuracy and repeatability of data generated from FFPE specimens using high throughput genomic assays. Together, these studies demonstrate feasibility and provide crucial guidance for future studies using FFPE specimens. Here, we summarize the findings of these studies and discuss the limitations of high throughput data generated from FFPE specimens across several platforms that include microarray, high throughput sequencing, and NanoString.

Performance comparison of two commercial human whole-exome capture systems on formalin-fixed paraffin-embedded lung adenocarcinoma samples

Background

Next Generation Sequencing (NGS) has become a valuable tool for molecular landscape characterization of cancer genomes, leading to a better understanding of tumor onset and progression, and opening new avenues in translational oncology. Formalin-fixed paraffin-embedded (FFPE) tissue is the method of choice for storage of clinical samples, however low quality of FFPE genomic DNA (gDNA) can limit its use for downstream applications.

Methods

To investigate the FFPE specimen suitability for NGS analysis and to establish the performance of two solution-based exome capture technologies, we compared the whole-exome sequencing (WES) data of gDNA extracted from 5 fresh frozen (FF) and 5 matched FFPE lung adenocarcinoma tissues using: SeqCap EZ Human Exome v.3.0 (Roche NimbleGen) and SureSelect XT Human All Exon v.5 (Agilent Technologies).

Results

Sequencing metrics on Illumina HiSeq were optimal for both exome systems and comparable among FFPE and FF samples, with a slight increase of PCR duplicates in FFPE, mainly in Roche NimbleGen libraries. Comparison of single nucleotide variants (SNVs) between FFPE-FF pairs reached overlapping values >90 % in both systems. Both WES showed high concordance with target re-sequencing data by Ion PGM™ in 22 lung-cancer genes, regardless the source of samples. Exon coverage of 623 cancer-related genes revealed high coverage efficiency of both kits, proposing WES as a valid alternative to target re-sequencing.

Conclusions

High-quality and reliable data can be successfully obtained from WES of FFPE samples starting from a relatively low amount of input gDNA, suggesting the inclusion of NGS-based tests into clinical contest. In conclusion, our analysis suggests that the WES approach could be extended to a translational research context as well as to the clinic (e.g. to study rare malignancies), where the simultaneous analysis of the whole coding region of the genome may help in the detection of cancer-linked variants.

Electronic supplementary material

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

Depletion of Mouse Cells from Human Tumor Xenografts Significantly Improves Downstream Analysis of Target Cells

The use of in vitro cell line models for cancer research has been a useful tool. However, it has been shown that these models fail to reliably mimic patient tumors in different assays(1). Human tumor xenografts represent the gold standard with respect to tumor biology, drug discovery, and metastasis research (2-4). Tumor xenografts can be derived from different types of material like tumor cell lines, tumor tissue from primary patient tumors(4) or serially transplanted tumors. When propagated in vivo, xenografted tissue is infiltrated and vascularized by cells of mouse origin. Multiple factors such as the tumor entity, the origin of xenografted material, growth rate and region of transplantation influence the composition and the amount of mouse cells present in tumor xenografts. However, even when these factors are kept constant, the degree of mouse cell contamination is highly variable. Contaminating mouse cells significantly impair downstream analyses of human tumor xenografts. As mouse fibroblasts show high plating efficacies and proliferation rates, they tend to overgrow cultures of human tumor cells, especially slowly proliferating subpopulations. Mouse cell derived DNA, mRNA, and protein components can bias downstream gene expression analysis, next-generation sequencing, as well as proteome analysis (5). To overcome these limitations, we have developed a fast and easy method to isolate untouched human tumor cells from xenografted tumor tissue. This procedure is based on the comprehensive depletion of cells of mouse origin by combining automated tissue dissociation with the benchtop tissue dissociator and magnetic cell sorting. Here, we demonstrate that human target cells can be can be obtained with purities higher than 96% within less than 20 min independent of the tumor type.

Genomic insights into head and neck cancer

Head and neck squamous cell carcinoma (HNSCC) is the sixth most common cancer worldwide and is frequently impervious to curative treatment efforts. Similar to other cancers associated with prolonged exposure to carcinogens, HNSCCs often have a high burden of mutations, contributing to substantial inter- and intra-tumor heterogeneity. The heterogeneity of this malignancy is further increased by the rising rate of human papillomavirus (HPV)-associated (HPV+) HNSCC, which defines an etiological subtype significantly different from the more common tobacco and alcohol associated HPV-negative (HPV-) HNSCC. Since 2011, application of large scale genome sequencing projects by The Cancer Genome Atlas (TCGA) network and other groups have established extensive datasets to characterize HPV- and HPV+ HNSCC, providing a foundation for advanced molecular diagnoses, identification of potential biomarkers, and therapeutic insights. Some genomic lesions are now appreciated as widely dispersed. For example, HPV- HNSCC characteristically inactivates the cell cycle suppressors TP53 (p53) and CDKN2A (p16), and often amplifies CCND1 (cyclin D), which phosphorylates RB1 to promote cell cycle progression from G1 to S. By contrast, HPV+ HNSCC expresses viral oncogenes E6 and E7, which inhibit TP53 and RB1, and activates the cell cycle regulator E2F1. Frequent activating mutations in PIK3CA and inactivating mutations in NOTCH1 are seen in both subtypes of HNSCC, emphasizing the importance of these pathways. Studies of large patient cohorts have also begun to identify less common genetic alterations, predominantly found in HPV- tumors, which suggest new mechanisms relevant to disease pathogenesis. Targets of these alterations including AJUBA and FAT1, both involved in the regulation of NOTCH/CTNNB1 signaling. Genes involved in oxidative stress, particularly CUL3, KEAP1 and NFE2L2, strongly associated with smoking, have also been identified, and are less well understood mechanistically. Application of sophisticated data-mining approaches, integrating genomic information with profiles of tumor methylation and gene expression, have helped to further yield insights, and in some cases suggest additional approaches to stratify patients for clinical treatment. We here discuss some recent insights built on TCGA and other genomic foundations.

Comparison of whole-exome sequencing of matched fresh and formalin fixed paraffin embedded melanoma tumours: implications for clinical decision making

The identification of recurrent driver mutations by whole-exome sequencing (WES) of fresh-frozen human cancers and the subsequent development of novel targeted therapies have recently transformed the treatment of many cancers including melanoma. In routine clinical practice, fresh-frozen tissue is rarely available and mutation testing usually needs to be carried out on archival formalin fixed, paraffin embedded (FFPE) tissue, from which DNA is typically fragmented, cross-linked and of lower quality. In this study we aimed to determine whether WES data generated from genomic DNA (gDNA) extracted from FFPE tissues can be produced reliably and of clinically-actionable standard. In this study of ten melanoma patients, we compared WES data produced from analysis of gDNA isolated from FFPE tumour tissue with that isolated from fresh-frozen tumour tissue from the same specimen. FFPE samples were sequenced using both Illumina's Nextera and NimbleGen SeqCap exome capture kits. To examine mutations between the two tissue sources and platforms, somatic mutations in the FFPE exomes were called using the matched fresh tissue sequence as a reference. Of the 10 FFPE DNA samples, seven Nextera and four SeqCap samples passed library preparation. On average, there were 5341 and 2246 variants lost in FFPE compared to matched fresh tissue utilising Nextera and SeqCap kits, respectively. In order to explore the feasibility of future clinical implementation of WES, FFPE variants in 27 genes of important clinical relevance in melanoma were assessed. The average concordance rate was 43.2% over a total of 1299 calls for the chosen genes in the FFPE DNA. For the current clinically most important melanoma mutations, 0/3 BRAF and 6/8 (75%) NRAS FFPE calls were concordant with the fresh tissue result, which was confirmed using a Sequenom OncoCarta Panel. The poor performance of FFPE WES indicates that specialised library construction to account for low quality DNA and further refinements will be necessary before this approach could be used for routine clinical decision making over currently preferred techniques.

The sialyl-glycolipid stage-specific embryonic antigen 4 marks a subpopulation of chemotherapy-resistant breast cancer cells with mesenchymal features

Introduction

Chemotherapy resistance resulting in incomplete pathologic response is associated with high risk of metastasis and early relapse in breast cancer. The aim of this study was to identify and evaluate biomarkers of treatment-resistant tumor cells.

Methods

We performed a cell surface marker screen in triple-negative breast cancer patient-derived xenograft models treated with standard care genotoxic chemotherapy. Global expression profiling was used to further characterize the identified treatment-resistant subpopulations.

Results

High expression of sialyl-glycolipid stage-specific embryonic antigen 4 (SSEA4) was found in residual tumor cells surviving chemotherapy and in samples from metastatic patients who relapsed after neoadjuvant chemotherapy. Gene and microRNA (miRNA) expression profiling linked SSEA4 positivity with a mesenchymal phenotype and a deregulation of drug resistance pathways. Functional assays demonstrated a direct link between epithelial–mesenchymal transition (EMT) and SSEA4 expression. Interestingly, SSEA4 expression, EMT, and drug resistance seemed to be regulated posttranscriptionally. Finally, high expression of CMP-N-acetylneuraminate-β-galactosamide-α-2,3-sialyltransferase 2 (ST3GAL2), the rate-limiting enzyme of SSEA4 synthesis, was found to be associated with poor clinical outcome in breast and ovarian cancer patients treated with chemotherapy.

Conclusions

In this study, we identified SSEA4 as highly expressed in a subpopulation of tumor cells resistant to multiple commonly used chemotherapy drugs, as well as ST3GAL2, the rate-limiting enzyme of SSEA4 synthesis, as a predictive marker of poor outcome for breast and ovarian cancer patients undergoing chemotherapy. Both biomarkers and additionally identified regulatory miRNAs may be used to further understand chemoresistance, to stratify patient groups in order to avoid ineffective and painful therapies, and to develop alternative treatment regimens for breast cancer patients.

Electronic supplementary material

The online version of this article (doi:10.1186/s13058-015-0652-6) contains supplementary material, which is available to authorized users.

Whole exome sequencing (WES) on formalin-fixed, paraffin-embedded (FFPE) tumor tissue in gastrointestinal stromal tumors (GIST)

BACKGROUND

Next generation sequencing (NGS) technology has been rapidly introduced into basic and translational research in oncology, but the reduced availability of fresh frozen (FF) tumor tissues and the poor quality of DNA extracted from formalin-fixed, paraffin-embedded (FFPE) has significantly impaired this process in the field of solid tumors. To evaluate if data generated from FFPE material can be reliably produced and potentially used in routine clinical settings, we performed whole exome sequencing (WES) from tumor samples of Gastrointestinal stromal tumors (GIST), either extracted FF or FFPE, and from matched normal DNA.

METHODS

We performed whole exome enrichment and sequencing at 100bp in paired end on four GIST samples, either from FFPE or fresh-frozen tissue, and from matched normal DNA.

RESULTS

The integrity of DNA extracted from FFPE was evaluated by a modified RAPD PCR method, thus identifying high quality (HQ) and low quality (LQ) FFPE. DNA library production and exome capture was feasible for both classes of FFPE, despite the smaller yield and insert size of LQ-FFPE. WES produced data of equal quality from FF and FFPE, while only HQ-FFPE yielded an amount of data comparable to FF samples. Bioinformatic analysis showed that the percentage of variants called both in FF and FFPE samples was very high in HQ-FFPE, reaching 94-96 % of the total number of called variants. Classification of somatic variants by nucleotide substitution type showed that HQ-FFPE and FF had similar mutational profiles, while LQ-FFPE samples carried a much higher number of mutations than the FF counterpart, with a significant enrichment of C > T/G > A substitutions. Focusing on potential disease-related variants allowed the discovery of additional somatic variants in GIST samples, apart from the known oncogenic driver mutation, both from sequencing of FF and FFPE material. False positive and false negative calls were present almost exclusively in the analysis of FFPE of low quality. On the whole this study showed that WES is feasible also on FFPE specimens and that it is possible to easily select FFPE samples of high quality that yield sequencing results comparable to the FF counterpart.

CONCLUSIONS

WES on FFPE material may represent an important and innovative source for GIST research and for other solid tumors, amenable of possible application in clinical practice.

HaloPlex Targeted Resequencing for Mutation Detection in Clinical Formalin-Fixed, Paraffin-Embedded Tumor Samples

In recent years, the advent of massively parallel next-generation sequencing technologies has enabled substantial advances in the study of human diseases. Combined with targeted DNA enrichment methods, high sequence coverage can be obtained for different genes simultaneously at a reduced cost per sample, creating unique opportunities for clinical cancer diagnostics. However, the formalin-fixed, paraffin-embedded (FFPE) process of tissue samples, routinely used in pathology departments, results in DNA fragmentation and nucleotide modifications that introduce a number of technical challenges for downstream biomolecular analyses. We evaluated the HaloPlex target enrichment system for somatic mutation detection in 80 tissue fractions derived from 20 clinical cancer cases with paired tumor and normal tissue available in both FFPE and fresh-frozen format. Several modifications to the standard method were introduced, including a reduced target fragment length and two strand capturing. We found that FFPE material can be used for HaloPlex-based target enrichment and next-generation sequencing, even when starting from small amounts of DNA. By specifically capturing both strands for each target fragment, we were able to reduce the number of false-positive errors caused by FFPE-induced artifacts and lower the detection limit for somatic mutations. We believe that the HaloPlex method presented here will be broadly applicable as a tool for somatic mutation detection in clinical cancer settings.

Development of a Targeted Multi-Disorder High-Throughput Sequencing Assay for the Effective Identification of Disease-Causing Variants

BACKGROUND

While next generation sequencing (NGS) is a useful tool for the identification of genetic variants to aid diagnosis and support therapy decision, high sequencing costs have limited its application within routine clinical care, especially in economically depressed areas. To investigate the utility of a multi-disease NGS based genetic test, we designed a custom sequencing assay targeting over thirty disease-associated areas including cardiac disorders, intellectual disabilities, hearing loss, collagenopathies, muscular dystrophy, Ashkenazi Jewish genetic disorders, and complex Mendelian disorders. We focused on these specific areas based on the interest of our collaborative clinical team, suggesting these diseases being the ones in need for the development of a sequencing-screening assay.

RESULTS

We targeted all coding, untranslated regions (UTR) and flanking intronic regions of 650 known disease-associated genes using the Roche-NimbleGen EZ SeqCapV3 capture system and sequenced on the Illumina HiSeq 2500 Rapid Run platform. Eight controls with known variants and one HapMap sample were first sequenced to assess the performance of the panel. Subsequently, as a proof of principle and to explore the possible utility of our test, we analyzed test disease subjects (n = 16). Eight had known Mendelian disorders and eight had complex pediatric diseases. In addition to assess whether copy number variation may be of utility as a companion assay relative to these specific disease areas, we used the Affymetrix Genome-Wide SNP Array 6.0 to analyze the same samples.

CONCLUSION

We identified potentially disease-associated variants: 22 missense, 4 nonsense, 1 frameshift, and 1 splice variants (16 previously identified, 12 novel among dbSNP and 15 novel among NHLBI Exome Variant Server). We found multi-disease targeted high-throughput sequencing to be a cost efficient approach in detecting disease-associated variants to aid diagnosis.

Appraisal of the technologies and review of the genomic landscape of ductal carcinoma in situ of the breast

Ductal carcinoma in situ is a biologically diverse entity. Whereas some lesions are cured by local surgical excision, others recur as in situ disease or progress to invasive carcinoma with subsequent potential for metastatic spread. Reliable prognostic biomarkers are therefore desirable for appropriate clinical management but remain elusive. In common with invasive breast cancer, ductal carcinoma in situ exhibits many genomic changes, predominantly copy number alterations. Although studies have revealed the genomic heterogeneity within individual ductal carcinoma in situ lesions and the association of certain copy number alterations with nuclear grade, none of the genomic changes defined so far is consistently associated with invasive transformation or recurrence risk in pure ductal carcinoma in situ. This article will review the current landscape of genomic alterations in ductal carcinoma in situ and their potential as prognostic biomarkers together with the technologies used to define these.

The efficacy of uracil DNA glycosylase pretreatment in amplicon-based massively parallel sequencing with DNA extracted from archived formalin-fixed paraffin-embedded esophageal cancer tissues

Advances in mutation testing for molecular-targeted cancer therapies have led to the increased use of archived formalin-fixed paraffin-embedded (FFPE) tumors. However, DNA extracted from FFPE tumors (FFPE DNA) is problematic for mutation testing, especially for amplicon-based massively parallel sequencing (MPS), owing to DNA fragmentation and artificial C:G > T:A single nucleotide variants (SNVs) caused by deamination of cytosine to uracil. Therefore, to reduce artificial C:G > T:A SNVs in amplicon-based MPS using FFPE DNA, we evaluated the efficacy of uracil DNA glycosylase (UDG) pretreatment, which can eliminate uracil-containing DNA molecules, with 126 archived FFPE esophageal cancer specimens. We also examined the association between the frequency of C:G > T:A SNVs and DNA quality, as assessed by a quantitative PCR (qPCR)-based assay. UDG pretreatment significantly lowered the frequency of C:G > T:A SNVs in highly fragmented DNA (by approximately 60%). This effect was not observed for good- to moderate-quality DNA, suggesting that a predictive assay (i.e., DNA quality assessment) needs to be performed prior to UDG pretreatment. These results suggest that UDG pretreatment is efficacious for mutation testing by amplicon-based MPS with fragmented DNA from FFPE samples.

Molecular methods for somatic mutation testing in lung adenocarcinoma: EGFR and beyond

Somatic mutational profiling in cancer has revolutionized the practice of clinical oncology. The discovery of driver mutations in non-small cell lung cancer (NSCLC) is an example of this. Molecular testing of lung adenocarcinoma is now considered standard of care and part of the diagnostic algorithm. This article provides an overview of the workflow of molecular testing in a clinical diagnostic laboratory discussing in particular novel assays that are currently in use for somatic mutation detection in NSCLC focussing on epidermal growth factor receptor (EGFR) mutations and anaplastic lymphoma kinase (ALK), ROS1 and RET rearrangements.

ClinSeK: a targeted variant characterization framework for clinical sequencing

Applying genomics to patient care demands sensitive, unambiguous and rapid characterization of a known set of clinically relevant variants in patients' samples, an objective substantially different from the standard discovery process, in which every base in every sequenced read must be examined. Further, the approach must be sufficiently robust as to be able to detect multiple and potentially rare variants from heterogeneous samples. To meet this critical objective, we developed a novel variant characterization framework, ClinSeK, which performs targeted analysis of relevant reads from high-throughput sequencing data. ClinSeK is designed for efficient targeted short read alignment and is capable of characterizing a wide spectrum of genetic variants from single nucleotide variation to large-scale genomic rearrangement breakpoints. Applying ClinSeK to over a thousand cancer patients demonstrated substantively better performance, in terms of accuracy, runtime and disk storage, for clinical applications than existing variant discovery tools. ClinSeK is freely available for academic use at http://bioinformatics.mdanderson.org/main/clinsek.

Sequence Artifacts in DNA from Formalin-Fixed Tissues: Causes and Strategies for Minimization

BACKGROUND

Precision medicine is dependent on identifying actionable mutations in tumors. Accurate detection of mutations is often problematic in formalin-fixed paraffin-embedded (FFPE) tissues. DNA extracted from formalin-fixed tissues is fragmented and also contains DNA lesions that are the sources of sequence artifacts. Sequence artifacts can be difficult to distinguish from true mutations, especially in the context of tumor heterogeneity, and are an increasing interpretive problem in this era of massively parallel sequencing. Understanding of the sources of sequence artifacts in FFPE tissues and implementation of preventative strategies are critical to improve the accurate detection of actionable mutations.

CONTENT

This mini-review focuses on DNA template lesions in FFPE tissues as the source of sequence artifacts in molecular analysis. In particular, fragmentation, base modification (including uracil and thymine deriving from cytosine deamination), and abasic sites are discussed as indirect or direct sources of sequence artifacts. We discuss strategies that can be implemented to minimize sequence artifacts and to distinguish true mutations from sequence artifacts. These strategies are applicable for the detection of actionable mutations in both single amplicon and massively parallel amplicon sequencing approaches.

SUMMARY

Because FFPE tissues are usually the only available material for DNA analysis, it is important to maximize the accurate informational content from FFPE DNA. Careful consideration of each step in the work flow is needed to minimize sequence artifacts. In addition, validation of actionable mutations either by appropriate experimental design or by orthogonal methods should be considered.

Sequencing approaches in cancer treatment

Use of sequencing approaches is an important aspect in the field of cancer genomics, where next-generation sequencing has already been utilized for targeting oncogenes or tumour-suppressor genes, that can be sequenced in a short time period. Alterations such as point mutations, insertions/deletions, copy number alterations, chromosomal rearrangements and epigenetic changes are encountered in cancer cell genomes, and application of various NGS technologies in cancer research will encounter such modifications. Rapid advancement in technology has led to exponential growth in the field of genomic analysis. The $1000 Genome Project (in which the goal is to sequence an entire human genome for $1000), and deep sequencing techniques (which have greater accuracy and provide a more complete analysis of the genome), are examples of rapid advancements in the field of cancer genomics. In this mini review, we explore sequencing techniques, correlating their importance in cancer therapy and treatment.

Sequence artefacts in a prospective series of formalin-fixed tumours tested for mutations in hotspot regions by massively parallel sequencing

Background

Clinical specimens undergoing diagnostic molecular pathology testing are fixed in formalin due to the necessity for detailed morphological assessment. However, formalin fixation can cause major issues with molecular testing, as it causes DNA damage such as fragmentation and non-reproducible sequencing artefacts after PCR amplification. In the context of massively parallel sequencing (MPS), distinguishing true low frequency variants from sequencing artefacts remains challenging. The prevalence of formalin-induced DNA damage and its impact on molecular testing and clinical genomics remains poorly understood.

Methods

The Cancer 2015 study is a population-based cancer cohort used to assess the feasibility of mutational screening using MPS in cancer patients from Victoria, Australia. While blocks were formalin-fixed and paraffin-embedded in different anatomical pathology laboratories, they were centrally extracted for DNA utilising the same protocol, and run through the same MPS platform (Illumina TruSeq Amplicon Cancer Panel). The sequencing artefacts in the 1-10% and the 10-25% allele frequency ranges were assessed in 488 formalin-fixed tumours from the pilot phase of the Cancer 2015 cohort. All blocks were less than 2.5 years of age (mean 93 days).

Results

Consistent with the signature of DNA damage due to formalin fixation, many formalin-fixed samples displayed disproportionate levels of C>T/G>A changes in the 1-10% allele frequency range. Artefacts were less apparent in the 10-25% allele frequency range. Significantly, changes were inversely correlated with coverage indicating high levels of sequencing artefacts were associated with samples with low amounts of available amplifiable template due to fragmentation. The degree of fragmentation and sequencing artefacts differed between blocks sourced from different anatomical pathology laboratories. In a limited validation of potentially actionable low frequency mutations, a NRAS G12D mutation in a melanoma was shown to be a false positive.

Conclusions

These findings indicate that DNA damage following formalin fixation remains a major challenge in laboratories working with MPS. Methodologies that assess, minimise or remove formalin-induced DNA damaged templates as part of MPS protocols will aid in the interpretation of genomic results leading to better patient outcomes.

Bioinformatics pipelines for targeted resequencing and whole-exome sequencing of human and mouse genomes: a virtual appliance approach for instant deployment

Targeted resequencing by massively parallel sequencing has become an effective and affordable way to survey small to large portions of the genome for genetic variation. Despite the rapid development in open source software for analysis of such data, the practical implementation of these tools through construction of sequencing analysis pipelines still remains a challenging and laborious activity, and a major hurdle for many small research and clinical laboratories. We developed TREVA (Targeted REsequencing Virtual Appliance), making pre-built pipelines immediately available as a virtual appliance. Based on virtual machine technologies, TREVA is a solution for rapid and efficient deployment of complex bioinformatics pipelines to laboratories of all sizes, enabling reproducible results. The analyses that are supported in TREVA include: somatic and germline single-nucleotide and insertion/deletion variant calling, copy number analysis, and cohort-based analyses such as pathway and significantly mutated genes analyses. TREVA is flexible and easy to use, and can be customised by Linux-based extensions if required. TREVA can also be deployed on the cloud (cloud computing), enabling instant access without investment overheads for additional hardware. TREVA is available at http://bioinformatics.petermac.org/treva/.