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

Development and Analytical Validation of a DNA Dual-Strand Approach for the US Food and Drug Administration-Approved Next-Generation Sequencing-Based Praxis Extended RAS Panel for Metastatic Colorectal Cancer Samples

A next-generation sequencing method was developed that can distinguish single-stranded modifications from low-frequency somatic mutations present on both strands of DNA in formalin-fixed paraffin-embedded colorectal cancer samples. We applied this method for analytical validation of the Praxis Extended RAS Panel, a US Food and Drug Administration-approved companion diagnostic for panitumumab, on the Illumina MiSeqDx platform. With the use of the TruSeq amplicon workflow, both strands of DNA from the starting material were interrogated independently. Mutations were reported only if found on both strands; artifacts usually present on only one strand would not be reported. A total of 56 mutations were targeted within the KRAS and NRAS genes. A minimum read depth of 1800× per amplicon is required per sample but averaged >30,000× at maximum multiplexing levels. Analytical validation studies were performed to determine the simultaneous detection of mutations on both strands, reproducibility, assay detection level, precision of the assay across various factors, and the impact of interfering substances. In conclusion, this assay can clearly distinguish single-stranded artifacts from low-frequency mutations. Furthermore, the assay is accurate, precise, and reproducible, can achieve consistent detection of a mutation at 5% mutation frequency, exhibits minimal impact from tested interfering substances, and can simultaneously detect 56 mutations in a single run using 10 samples plus controls.

Reliability of Whole-Exome Sequencing for Assessing Intratumor Genetic Heterogeneity

Multi-region sequencing is used to detect intratumor genetic heterogeneity (ITGH) in tumors. To assess whether genuine ITGH can be distinguished from sequencing artifacts, we performed whole-exome sequencing (WES) on three anatomically distinct regions of the same tumor with technical replicates to estimate technical noise. Somatic variants were detected with three different WES pipelines and subsequently validated by high-depth amplicon sequencing. The cancer-only pipeline was unreliable, with about 69% of the identified somatic variants being false positive. Even with matched normal DNA for which 82% of the somatic variants were detected reliably, only 36%-78% were found consistently in technical replicate pairs. Overall, 34%-80% of the discordant somatic variants, which could be interpreted as ITGH, were found to constitute technical noise. Excluding mutations affecting low-mappability regions or occurring in certain mutational contexts was found to reduce artifacts, yet detection of subclonal mutations by WES in the absence of orthogonal validation remains unreliable.

Comparative Genomics Reveals Shared Mutational Landscape in Canine Hemangiosarcoma and Human Angiosarcoma

Angiosarcoma is a highly aggressive cancer of blood vessel-forming cells with few effective treatment options and high patient mortality. It is both rare and heterogenous, making large, well-powered genomic studies nearly impossible. Dogs commonly suffer from a similar cancer, called hemangiosarcoma, with breeds like the golden retriever carrying heritable genetic factors that put them at high risk. If the clinical similarity of canine hemangiosarcoma and human angiosarcoma reflects shared genomic etiology, dogs could be a critically needed model for advancing angiosarcoma research. We assessed the genomic landscape of canine hemangiosarcoma via whole-exome sequencing (47 golden retriever hemangiosarcomas) and RNA sequencing (74 hemangiosarcomas from multiple breeds). Somatic coding mutations occurred most frequently in the tumor suppressor TP53 (59.6% of cases) as well as two genes in the PI3K pathway: the oncogene PIK3CA (29.8%) and its regulatory subunit PIK3R1 (8.5%). The predominant mutational signature was the age-associated deamination of cytosine to thymine. As reported in human angiosarcoma, CDKN2A/B was recurrently deleted and VEGFA, KDR, and KIT recurrently gained. We compared the canine data to human data recently released by The Angiosarcoma Project, and found many of the same genes and pathways significantly enriched for somatic mutations, particularly in breast and visceral angiosarcomas. Canine hemangiosarcoma closely models the genomic landscape of human angiosarcoma of the breast and viscera, and is a powerful tool for investigating the pathogenesis of this devastating disease. IMPLICATIONS: We characterize the genomic landscape of canine hemangiosarcoma and demonstrate its similarity to human angiosarcoma.

Systematic evaluation of PAXgene® tissue fixation for the histopathological and molecular study of lung cancer

Whilst adequate for most existing pathological tests, formalin is generally considered a poor DNA preservative and use of alternative fixatives may prove advantageous for molecular testing of tumour material; an increasingly common approach to identify targetable driver mutations in lung cancer patients. We collected paired PAXgene® tissue-fixed and formalin-fixed samples of block-sized tumour and lung parenchyma, Temno-needle core tumour biopsies and fine needle tumour aspirates (FNAs) from non-small cell lung cancer resection specimens. Traditionally processed formalin fixed paraffin wax embedded (FFPE) samples were compared to paired PAXgene® tissue fixed paraffin-embedded (PFPE) samples. We evaluated suitability for common laboratory tests (H&E staining and immunohistochemistry) and performance for downstream molecular investigations relevant to lung cancer, including RT-PCR and next generation DNA sequencing (NGS). Adequate and comparable H&E staining was seen in all sample types and nuclear staining was preferable in PAXgene® fixed Temno tumour biopsies and tumour FNA samples. Immunohistochemical staining was broadly comparable. PFPE samples enabled greater yields of less-fragmented DNA than FFPE comparators. PFPE samples were also superior for PCR and NGS performance, both in terms of quality control metrics and for variant calling. Critically we identified a greater number of genetic variants in the epidermal growth factor receptor gene when using PFPE samples and the Ingenuity® Variant Analysis pipeline. In summary, PFPE samples are adequate for histopathological diagnosis and suitable for the majority of existing laboratory tests. PAXgene® fixation is superior for DNA and RNA integrity, particularly in low-yield samples and facilitates improved NGS performance, including the detection of actionable lung cancer mutations for precision medicine in lung cancer samples.

Mutation of the PTCH1 gene predicts recurrence of breast cancer

Breast cancer is the most common cancer in women, and some patients develop recurrence after standard therapy. Effective predictors are urgently needed to detect recurrence earlier. The activation of Hedgehog signaling in breast cancer is correlated with poor prognosis. PTCH1 is an essential membrane receptor of Hedgehog. However, there are few reports about mutations in Hedgehog genes in breast cancer. We conducted a comprehensive study via an experimental and bioinformatics approach to detect mutated genes in breast cancer. Twenty-two breast cancer patients who developed recurrence within 24 months postoperatively were enrolled with 22 control cancer patients. Targeted deep sequencing was performed to assess the mutations among individuals with breast cancer using a panel of 143 cancer-associated genes. Bioinformatics and public databases were used to predict the protein functions of the mutated genes. Mutations were identified in 44 breast cancer specimens, and the most frequently mutated genes were BRCA2, APC, ATM, BRCA1, NF1, TET2, TSC1, TSC2, NOTCH1, MSH2, PTCH1, TP53, PIK3CA, FBXW7, and RB1. Mutation of these genes was correlated with protein phosphorylation and autophosphorylation, such as peptidyl-tyrosine and protein kinase C phosphorylation. Among these highly mutated genes, mutations of PTCH1 were associated with poor prognosis and increased recurrence of breast cancer, especially mutations in exons 22 and 23. The public sequencing data from the COSMIC database were exploited to predict the functions of the mutations. Our findings suggest that mutation of PTCH1 is correlated with early recurrence of breast cancer patients and will become a powerful predictor for recurrence of breast cancer.

Current and emerging biomarkers in metastatic colorectal cancer

Background

The incorporation of novel biomarkers into therapy selection for patients with metastatic colorectal cancer (mcrc) has significantly improved outcomes. Optimal treatment planning now takes into account diverse characteristics of patients and their tumours to create personalized therapeutic plans.

Discussion

This review is split into two sections. In the first section, we review the prognostic and predictive significance of expanded RAS mutation testing, BRAF mutations, ERBB2 (her2) amplification, microsatellite instability (msi) and deficient mismatch repair (dmmr) protein, NTRK fusions, PIK3CA mutations, and met amplifications. The therapeutic implication of each of those biomarkers for personalizing therapies for each patient with mcrc is discussed. In the second section, we touch on testing methods and considerations of relevance to clinicians when they interpret companion diagnostics meant to guide therapy selection. The advantages and pitfalls of various methods are evaluated, and we also look at the potential of liquid biopsies and circulating tumour dna (ctdna) to change the landscape of therapeutic choice and biologic understanding of the disease.

Summary

Routine testing for extended RAS, BRAF, dmmr or high msi, and NTRK fusions is necessary to determine the best sequencing of chemotherapy and biologic agents for patients with mcrc. Although next-generation sequencing and ctdna are increasingly being adopted, other techniques such as immunohistochemistry retain their relevance in detection of her2 amplification, NTRK fusions, and dmmr.

Tumor mutational burden assessed by targeted NGS predicts clinical benefit from immune checkpoint inhibitors in non-small cell lung cancer

In non‐small cell lung cancer (NSCLC), immune checkpoint inhibitors (ICIs) significantly improve overall survival (OS). Tumor mutational burden (TMB) has emerged as a predictive biomarker for patients treated with ICIs. Here, we evaluated the predictive power of TMB measured by the Oncomine™ Tumor Mutational Load targeted sequencing assay in 76 NSCLC patients treated with ICIs. TMB was assessed retrospectively in 76 NSCLC patients receiving ICI therapy. Clinical data (RECIST 1.1) were collected and patients were classified as having either durable clinical benefit (DCB) or no durable benefit (NDB). Additionally, genetic alterations and PD‐L1 expression were assessed and compared with TMB and response rate. TMB was significantly higher in patients with DCB than in patients with NDB (median TMB = 8.5 versus 6.0 mutations/Mb, Mann–Whitney p = 0.0244). 64% of patients with high TMB (cut‐off = third tertile, TMB ≥ 9) were responders (DCB) compared to 33% and 29% of patients with intermediate and low TMB, respectively (cut‐off = second and first tertile, TMB = 5–9 and TMB ≤ 4, respectively). TMB‐high patients showed significantly longer progression‐free survival (PFS) and OS (log‐rank test p = 0.0014 for PFS and 0.0197 for OS). While identifying different subgroups of patients, combining PD‐L1 expression and TMB increased the predictive power (from AUC 0.63 to AUC 0.65). Our results show that the TML panel is an effective tool to stratify patients for ICI treatment. A combination of biomarkers might maximize the predictive precision for patient stratification. Our study supports TMB evaluation through targeted NGS in NSCLC patient samples as a tool to predict response to ICI therapy. We offer recommendations for a reliable and cost‐effective assessment of TMB in a routine diagnostic setting. © 2019 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of Pathological Society of Great Britain and Ireland.

A novel genomic panel as an adjunctive diagnostic tool for the characterization and profiling of breast Fibroepithelial lesions

Background

Known collectively as breast fibroepithelial lesions (FELs), the common fibroadenomas (FAs) and the rarer phyllodes tumors (PTs) are a heterogenous group of biphasic neoplasms. Owing to limited tissue availability, inter-observer variability, overlapping histological features and heterogeneity of these lesions, diagnosing them accurately on core biopsies is challenging. As the choice management option depends on the histological diagnosis; a novel 16-gene panel assay was developed to improve the accuracy of preoperative diagnosis on core biopsy specimens.

Methods

Using this 16-gene panel, targeted amplicon-based sequencing was performed on 275 formalin-fixed, paraffin-embedded (FFPE) breast FEL specimens, archived at the Singapore General Hospital, from 2008 to 2012.

Results

In total, 167 FAs, 24 benign, 14 borderline and 6 malignant PTs, were profiled. Compared to FAs, PTs had significantly higher mutation rates in the TERT promoter (p <  0.001), RARA (p <  0.001), FLNA, RB1 and TP53 (p = 0.002, 0.020 and 0.018, respectively). In addition to a higher mutational count (p <  0.001), TERT promoter (p <  0.001), frameshift, nonsense and splice site (p = 0.001, < 0.001 and 0.043, respectively) mutations were also frequently observed in PTs.

A multivariate logistic regression model was built using these as variables and a predictive scoring system was developed. It classifies a FEL at low or high risk (score <  1 and ≥ 1, respectively) of being a PT. This scoring system has good discrimination (ROC area = 0.773, 95% CI: 0.70 to 0.85), calibration (p = 0.945) and is significant in predicting PTs (p <  0.001).

Conclusion

This novel study demonstrates the ability to extract DNA of sufficient quality and quantity for targeted sequencing from FFPE breast core biopsy specimens, along with their successful characterization and profiling using our customized 16-gene panel. Prospective work includes validating the utility of this promising 16-gene panel assay as an adjunctive diagnostic tool in clinical practice.

The molecular genetic make-up of male breast cancer

Male breast cancer (MBC) is extremely rare and accounts for less than 1% of all breast malignancies. Therefore, clinical management of MBC is currently guided by research on the disease in females. In this study, DNA obtained from 45 formalin-fixed paraffin-embedded (FFPE) MBCs with and 90 MBCs (52 FFPE and 38 fresh-frozen) without matched normal tissues was subjected to massively parallel sequencing targeting all exons of 1943 cancer-related genes. The landscape of mutations and copy number alterations was compared to that of publicly available estrogen receptor (ER)-positive female breast cancers (smFBCs) and correlated to prognosis. From the 135 MBCs, 90% showed ductal histology, 96% were ER-positive, 66% were progesterone receptor (PR)-positive, and 2% HER2-positive, resulting in 50, 46 and 4% luminal A-like, luminal B-like and basal-like cases, respectively. Five patients had Klinefelter syndrome (4%) and 11% of patients harbored pathogenic BRCA2 germline mutations. The genomic landscape of MBC to some extent recapitulated that of smFBC, with recurrent PIK3CA (36%) and GATA3 (15%) somatic mutations, and with 40% of the most frequently amplified genes overlapping between both sexes. TP53 (3%) somatic mutations were significantly less frequent in MBC compared to smFBC, whereas somatic mutations in genes regulating chromatin function and homologous recombination deficiency-related signatures were more prevalent. MDM2 amplifications were frequent (13%), correlated with protein overexpression (P = 0.001) and predicted poor outcome (P = 0.007). In conclusion, despite similarities in the genomic landscape between MBC and smFBC, MBC is a molecularly unique and heterogeneous disease requiring its own clinical trials and treatment guidelines.

Utility of Histologic and Histochemical Screening for 16S Ribosomal RNA Gene Sequencing of Formalin-Fixed, Paraffin-Embedded Tissue for Bacterial Endocarditis

OBJECTIVES

16S ribosomal RNA (rRNA) sequencing is a powerful but expensive tool for the identification of bacteria in culture-negative endocarditis. Histologic criteria to screen formalin-fixed, paraffin-embedded (FFPE) specimens for testing are evaluated.

METHODS

Sixty-eight cases of infective endocarditis and controls were histologically reviewed and analyzed by 16S rRNA gene sequencing.

RESULTS

Sequencing identified a specific pathogenic organism in 33 (49%) of 68 cases with acute inflammation and in 0 of 10 controls (P = .004). Visualization of organisms by Gram or Grocott methenamine silver stains had the strongest association with positive sequencing, while antibiotic treatment effect and acid decalcification decreased sensitivity. Molecular identifications were concordant with blood culture results in 90% of the cases, and a positive sequencing result was obtained in approximately half of the cases with negative valve cultures.

CONCLUSIONS

Histologic screening criteria are extremely helpful for identifying cases likely to be positive by molecular testing and can provide significant cost savings in filtering out low-yield specimens.

Bioinformatics and DNA-extraction strategies to reliably detect genetic variants from FFPE breast tissue samples

Background

Archived formalin fixed paraffin embedded (FFPE) samples are valuable clinical resources to examine clinically relevant morphology features and also to study genetic changes. However, DNA quality and quantity of FFPE samples are often sub-optimal, and resulting NGS-based genetics variant detections are prone to false positives. Evaluations of wet-lab and bioinformatics approaches are needed to optimize variant detection from FFPE samples.

Results

As a pilot study, we designed within-subject triplicate samples of DNA derived from paired FFPE and fresh frozen breast tissues to highlight FFPE-specific artifacts. For FFPE samples, we tested two FFPE DNA extraction methods to determine impact of wet-lab procedures on variant calling: QIAGEN QIAamp DNA Mini Kit (“QA”), and QIAGEN GeneRead DNA FFPE Kit (“QGR”). We also used negative-control (NA12891) and positive control samples (Horizon Discovery Reference Standard FFPE). All DNA sample libraries were prepared for NGS according to the QIAseq Human Breast Cancer Targeted DNA Panel protocol and sequenced on the HiSeq 4000. Variant calling and filtering were performed using QIAGEN Gene Globe Data Portal. Detailed variant concordance comparisons and mutational signature analysis were performed to investigate effects of FFPE samples compared to paired fresh frozen samples, along with different DNA extraction methods.

In this study, we found that five times or more variants were called with FFPE samples, compared to their paired fresh-frozen tissue samples even after applying molecular barcoding error-correction and default bioinformatics filtering recommended by the vendor. We also found that QGR as an optimized FFPE-DNA extraction approach leads to much fewer discordant variants between paired fresh frozen and FFPE samples. Approximately 92% of the uniquely called FFPE variants were of low allelic frequency range (< 5%), and collectively shared a “C > T|G > A” mutational signature known to be representative of FFPE artifacts resulting from cytosine deamination. Based on control samples and FFPE-frozen replicates, we derived an effective filtering strategy with associated empirical false-discovery estimates.

Conclusions

Through this study, we demonstrated feasibility of calling and filtering genetic variants from FFPE tissue samples using a combined strategy with molecular barcodes, optimized DNA extraction, and bioinformatics methods incorporating genomics context such as mutational signature and variant allelic frequency.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-6056-8) contains supplementary material, which is available to authorized users.

Sources of erroneous sequences and artifact chimeric reads in next generation sequencing of genomic DNA from formalin-fixed paraffin-embedded samples

Tissues used in pathology laboratories are typically stored in the form of formalin-fixed, paraffin-embedded (FFPE) samples. One important consideration in repurposing FFPE material for next generation sequencing (NGS) analysis is the sequencing artifacts that can arise from the significant damage to nucleic acids due to treatment with formalin, storage at room temperature and extraction. One such class of artifacts consists of chimeric reads that appear to be derived from non-contiguous portions of the genome. Here, we show that a major proportion of such chimeric reads align to both the ‘Watson’ and ‘Crick’ strands of the reference genome. We refer to these as strand-split artifact reads (SSARs). This study provides a conceptual framework for the mechanistic basis of the genesis of SSARs and other chimeric artifacts along with supporting experimental evidence, which have led to approaches to reduce the levels of such artifacts. We demonstrate that one of these approaches, involving S1 nuclease-mediated removal of single-stranded fragments and overhangs, also reduces sequence bias, base error rates, and false positive detection of copy number and single nucleotide variants. Finally, we describe an analytical approach for quantifying SSARs from NGS data.

EGFR T790M detection in formalin-fixed paraffin-embedded tissues of patients with lung cancer using RNA-based in situ hybridization: A preliminary feasibility study

Background

Following drug resistance in patients with lung cancer treated by EGFR TKIs, a biopsy is required to obtain sufficient cancer tissue for T790M detection in order to select potential beneficiaries suitable for third‐generation EGFR TKIs, such as osimertinib. The purpose of this study was to explore the feasibility of using a new in situ analysis technique based on RNA target sequences to detect EGFR T790M in lung cancer.

Methods

A total of 28 formalin‐fixed paraffin‐embedded (FFPE) samples from 24 lung adenocarcinoma patients archived in Peking Union Medical College Hospital from 2015 to 2016 were collected. The BaseScope T790M detection technique by in situ hybridization on FFPE slides was used to analyze the mutation of EGFR T790M, and then the results were compared with the data acquired by Scorpions ARMS assay, which is the so‐called gold standard for EGFR gene mutation testing. The sensitivity and specificity of the BaseScope T790M detection technique were preliminarily evaluated.

Results

Of the 28 FFPE specimens, the average proportion of T790M‐positive cells was 35.78% ± 17.68% in 18 samples with EGFR T790M, confirmed by Scorpions ARMS assay, Compared with real‐time PCR assay, the sensitivity and specificity of BaseScope T790M were all 100% in our cohort.

Conclusion

BaseScope T790M assay could be completed on only one FFPE slide and the visualized molecular result overplayed with histomorphological information perfectly, so it may be the alternative method for EGFR T790M evaluation. BaseScope assay has potential clinical utility, and it will be necessary to carry out validation with a large number of cases.

Pyrosequencing versus methylation-specific PCR for assessment of MGMT methylation in tumor and blood samples of glioblastoma patients

Circulating biomarkers in blood may provide an interesting alternative to risky tissue biopsies in the diagnosis and follow-up of glioblastoma patients. We have assessed MGMT methylation status in blood and tissue samples from unresected glioblastoma patients who had been included in the randomized GENOM-009 trial. Paired blood and tissue samples were assessed by methylation-specific PCR (MSP) and pyrosequencing (PYR). After establishing the minimum PYR cut-off that could yield a significant difference in overall survival, we assessed the sensitivity, specificity, positive predictive value and negative predictive value (NPV) of the analyses. Methylation could be detected in cfDNA by both MSP and PYR but with low concordance with results in tissue. Sensitivity was low for both methods (31% and 38%, respectively), while specificity was higher for MSP in blood than for PYR in plasma (96% vs 76%) and NPV was similar (56 vs 57%). Concordance of results in tissue by MSP and PYR was 84.3% (P < 0.001) and correlated with outcome. We conclude that detection of cfDNA in the blood of glioblastoma patients can be an alternative when tumor tissue is not available but methods for the detection of cfDNA in blood must improve before it can replace analysis in tumor tissue.

Mosaic abnormalities of the skin: review and guidelines from the European Reference Network for rare skin diseases

BACKGROUND

Cutaneous mosaicism is an area of dermatology in which there has been an explosion of knowledge within the current decade. This has led to fundamental changes in the understanding of the conditions in this field, and to an ongoing paradigm shift in the approach to management of mosaic skin disorders.

OBJECTIVES

To lay out the general principles of mosaicism as they are currently understood, summarize the known cutaneous mosaic abnormalities of the skin with associated phenotypic and genotypic information, review the latest trials on targeted therapies and propose guidelines for the general approach to a patient with suspected mosaicism.

METHODS

This was a consensus expert review as part of the European Reference Network project (ERN-Skin).

CONCLUSIONS

This study provides clinicians with a practical approach to the patient with suspected mosaicism, redefines mosaicism for the modern genetic era, and proposes a new classification system based on genetic mechanism. What's already known about this topic? Cutaneous mosaicism is a complex field of dermatology that encompasses most birthmarks, and many rare syndromes. Some cutaneous patterns are known to be seen in mosaicism. Very few treatment options are available for most mosaic abnormalities of the skin. Recent high-sensitivity genetic techniques have led to an explosion of knowledge about genotype and phenotype in the literature. What does this study add? Expert consensus from the European Reference Network project. Review of knowledge of confirmed mosaic abnormalities of the skin, including cutaneous phenotype, extracutaneous associated features and genotype. Proposed new classification of mosaic abnormalities of the skin by genetic mechanism and therefore inheritance potential. Practical tips on correct sample collection and genetic investigation. Review of trials of targeted therapies. Guidelines for a practical clinical approach to the patient with suspected mosaicism.

Is the Mirror Image Method Really Useful in Tumor Tissue Bank Quality Control?

Oncology research projects are highly dependent on the quality of tumor samples stored in the biobank. Microscopic control is important to ensure the quality of the frozen sample (Does the sample correspond to tumor tissue? Does the sample contain a sufficient number of tumor cells for molecular analysis?). The aim of this study was to evaluate the value of the mirror image method in quality control of colonic adenocarcinoma samples stored in a tumor bank. Microscopic concordance for the differentiation grade, malignant and normal cell percentages, necrosis, mucinous component, and ulceration was assessed on 82 colon adenocarcinoma banked samples and their paired, formalin-fixed, paraffin-embedded mirror controls. Molecular concordance for KRAS status was evaluated in 76 of these 82 cases. Morphological correspondence between frozen and mirror samples was good for the mucinous component (intraclass correlation coefficient [ICC] = 0.81), moderate for differentiation (Cohen's kappa coefficient [k] = 0.67), fair for malignant cells (ICC = 0.44), and poor for ulceration (k = 0.08), normal tissue (ICC = 0.36), and necrosis (ICC = 0.13) percentages. Molecular correspondence for KRAS status was almost perfect (95% correspondence, k = 0.88) between frozen and mirror samples. In conclusion, the mirror sample method is not a good alternative for microscopic and molecular control of frozen colonic adenocarcinoma samples.

Association Between Preanalytical Factors and Tumor Mutational Burden Estimated by Next-Generation Sequencing-Based Multiplex Gene Panel Assay

BACKGROUND

Tumor mutational burden (TMB) measured via next-generation sequencing (NGS)-based gene panel is a promising biomarker for response to immune checkpoint inhibitors (ICIs) in solid tumors. However, little is known about the preanalytical factors that can affect the TMB score.

MATERIALS AND METHODS

Data of 199 patients with solid tumors who underwent multiplex NGS gene panel (OncoPrime), which was commercially provided by a Clinical Laboratory Improvement Amendments-licensed laboratory and covered 0.78 megabase (Mb) of capture size relevant to the TMB calculation, were reviewed. Associations between the TMB score and preanalytical factors, including sample DNA quality, sample type, sampling site, and storage period, were analyzed. Clinical outcomes of patients with a high TMB score (≥10 mutations per megabase) who received anti-programmed cell death protein 1 antibodies (n = 22) were also analyzed.

RESULTS

Low DNA library concentration (<5 nM), formalin-fixed paraffin-embedded tissue (FFPE), and the prolonged sample storage period (range, 0.9-58.1 months) correlated with a higher TMB score. After excluding low DNA library samples from the analysis, FFPE samples, but not the sample storage period, exhibited a marked correlation with a high TMB score. Of 22 patients with a high TMB score, we observed the partial response in 2 patients (9.1%).

CONCLUSION

Our results indicate that the TMB score estimated via NGS-based gene panel could be affected by the DNA library concentration and sample type. These factors could potentially increase the false-positive and/or artifactual variant calls. As each gene panel has its own pipeline for variant calling, it is unknown whether these factors have a significant effect in other platforms.

IMPLICATIONS FOR PRACTICE

A high tumor mutational burden score, as estimated via next-generation sequencing-based gene panel testing, should be carefully interpreted as it could be affected by the DNA library concentration and sample type.

Whole-genome sequencing based on formalin-fixed paraffin-embedded endomyocardial biopsies for genetic studies on outcomes after heart transplantation

BACKGROUND

Whole-genome sequencing (WGS) of heart transplant recipient- and donor-derived cardiac biopsies may facilitate organ matching, graft failure prediction, and immunotolerance research. The objective of this study was to determine the feasibility of WGS based on formalin-fixed paraffin-embedded endomyocardial biopsies.

METHODS AND RESULTS

The study included serial donor- and recipient samples from patients who had undergone heart transplantation at Skane University Hospital, Lund, Sweden, between 1988 and 2009. DNA extraction and WGS were conducted. Additional WGS sequencing quality metrics and coverage were obtained with the Genome Analysis Toolkit (GATK). 455 endomyocardial samples from 37 heart transplant recipients were acquired from routine rejection monitoring and stored as formalin-fixed paraffin-embedded samples. They were analyzed after 3-26 years of storage. DNA was extracted from 114 samples and WGS was run on 85 samples. DNA extraction yielded 313 ng (IQR 96-601) for all samples. A coverage of 11.3x (IQR 9.0-15.9) was recorded for all WGS samples. Three samples stored for > 25 years yielded a coverage of > 25x. Data were generated for 1.7 billion reads per sample (IQR 1.4-2.7). A Transition/Transversion (TiTv) ratio of 2.09 ± 0.05 was calculated for all WGS samples. No associations were found among storage time, DNA yield, or sequencing quality metrics.

CONCLUSIONS

The present study demonstrated the feasibility of whole-genome sequencing based on endomyocardial biopsies. This process could enable large-scale retrospective genomic studies using stored histopathological samples.

Tumor mutational signatures in sebaceous skin lesions from individuals with Lynch syndrome

BACKGROUND

Muir-Torre syndrome is defined by the development of sebaceous skin lesions in individuals who carry a germline mismatch repair (MMR) gene mutation. Loss of expression of MMR proteins is frequently observed in sebaceous skin lesions, but MMR-deficiency alone is not diagnostic for carrying a germline MMR gene mutation.

METHODS

Whole exome sequencing was performed on three MMR-deficient sebaceous lesions from individuals with MSH2 gene mutations (Lynch syndrome) and three MMR-proficient sebaceous lesions from individuals without Lynch syndrome with the aim of characterizing the tumor mutational signatures, somatic mutation burden, and microsatellite instability status. Thirty predefined somatic mutational signatures were calculated for each lesion.

RESULTS

Signature 1 was ubiquitous across the six lesions tested. Signatures 6 and 15, associated with defective DNA MMR, were significantly more prevalent in the MMR-deficient lesions from the MSH2 carriers compared with the MMR-proficient non-Lynch sebaceous lesions (mean ± SD=41.0 ± 8.2% vs. 2.3 ± 4.0%, p = 0.0018). Tumor mutation burden was, on average, significantly higher in the MMR-deficient lesions compared with the MMR-proficient lesions (23.3 ± 11.4 vs. 1.8 ± 0.8 mutations/Mb, p = 0.03). All four sebaceous lesions observed in sun exposed areas of the body demonstrated signature 7 related to ultraviolet light exposure.

CONCLUSION

Tumor mutational signatures 6 and 15 and somatic mutation burden were effective in differentiating Lynch-related from non-Lynch sebaceous lesions.

Integration of transcriptional and mutational data simplifies the stratification of peripheral T-cell lymphoma

The histological diagnosis of peripheral T-cell lymphoma (PTCL) can represent a challenge, particularly in the case of closely related entities such as angioimmunoblastic T-lymphoma (AITL), PTCL-not otherwise specified (PTCL-NOS), and ALK-negative anaplastic large-cell lymphoma (ALCL). Although gene expression profiling and next generations sequencing have been proven to define specific features recurrently associated with distinct entities, genomic-based stratifications have not yet led to definitive diagnostic criteria and/or entered into the routine clinical practice. Herein, to improve the current molecular classification between AITL and PTCL-NOS, we analyzed the transcriptional profiles from 503 PTCLs stratified according to their molecular configuration and integrated them with genomic data of recurrently mutated genes (RHOA ^G17V , TET2, IDH2 ^R172 , and DNMT3A) in 53 cases (39 AITLs and 14 PTCL-NOSs) included in the series. Our analysis unraveled that the mutational status of RHOA ^G17V , TET2, and DNMT3A poorly correlated, individually, with peculiar transcriptional fingerprints. Conversely, in IDH2 ^R172 samples a strong transcriptional signature was identified that could act as a surrogate for mutational status. The integrated analysis of clinical, mutational, and molecular data led to a simplified 19-gene signature that retains high accuracy in differentiating the main nodal PTCL entities. The expression levels of those genes were confirmed in an independent cohort profiled by RNA-sequencing.

Evolutionary history of human colitis-associated colorectal cancer

OBJECTIVE

IBD confers an increased lifetime risk of developing colorectal cancer (CRC), and colitis-associated CRC (CA-CRC) is molecularly distinct from sporadic CRC (S-CRC). Here we have dissected the evolutionary history of CA-CRC using multiregion sequencing.

DESIGN

Exome sequencing was performed on fresh-frozen multiple regions of carcinoma, adjacent non-cancerous mucosa and blood from 12 patients with CA-CRC (n=55 exomes), and key variants were validated with orthogonal methods. Genome-wide copy number profiling was performed using single nucleotide polymorphism arrays and low-pass whole genome sequencing on archival non-dysplastic mucosa (n=9), low-grade dysplasia (LGD; n=30), high-grade dysplasia (HGD; n=13), mixed LGD/HGD (n=7) and CA-CRC (n=19). Phylogenetic trees were reconstructed, and evolutionary analysis used to reveal the temporal sequence of events leading to CA-CRC.

RESULTS

10/12 tumours were microsatellite stable with a median mutation burden of 3.0 single nucleotide alterations (SNA) per Mb, ~20% higher than S-CRC (2.5 SNAs/Mb), and consistent with elevated ageing-associated mutational processes. Non-dysplastic mucosa had considerable mutation burden (median 47 SNAs), including mutations shared with the neighbouring CA-CRC, indicating a precancer mutational field. CA-CRCs were often near triploid (40%) or near tetraploid (20%) and phylogenetic analysis revealed that copy number alterations (CNAs) began to accrue in non-dysplastic bowel, but the LGD/HGD transition often involved a punctuated 'catastrophic' CNA increase.

CONCLUSIONS

Evolutionary genomic analysis revealed precancer clones bearing extensive SNAs and CNAs, with progression to cancer involving a dramatic accrual of CNAs at HGD. Detection of the cancerised field is an encouraging prospect for surveillance, but punctuated evolution may limit the window for early detection.

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.

Targeted capture-based NGS is superior to multiplex PCR-based NGS for hereditary BRCA1 and BRCA2 gene analysis in FFPE tumor samples

BACKGROUND

With the introduction of Olaparib treatment for BRCA-deficient recurrent ovarian cancer, testing for somatic and/or germline mutations in BRCA1/2 genes in tumor tissues became essential for treatment decisions. In most cases only formalin-fixed paraffin-embedded (FFPE) samples, containing fragmented and chemically modified DNA of minor quality, are available. Thus, multiplex PCR-based sequencing is most commonly applied in routine molecular testing, which is predominantly focused on the identification of known hot spot mutations in oncogenes.

METHODS

We compared the overall performance of an adjusted targeted capture-based enrichment protocol and a multiplex PCR-based approach for calling of pathogenic SNVs and InDels using DNA extracted from 13 FFPE tissue samples. We further applied both strategies to seven blood samples and five matched FFPE tumor tissues of patients with known germline exon-spanning deletions and gene-wide duplications in BRCA1/2 to evaluate CNV detection based solely on panel NGS data. Finally, we analyzed DNA from FFPE tissues of 11 index patients from families suspected of having hereditary breast and ovarian cancer, of whom no blood samples were available for testing, in order to identify underlying pathogenic germline BRCA1/2 mutations.

RESULTS

The multiplex PCR-based protocol produced inhomogeneous coverage among targets of each sample and between samples as well as sporadic amplicon drop out, leading to insufficiently or non-covered nucleotides, which subsequently hindered variant detection. This protocol further led to detection of PCR-artifacts that could easily have been misinterpreted as pathogenic mutations. No such limitations were observed by application of an adjusted targeted capture-based protocol, which allowed for CNV calling with 86% sensitivity and 100% specificity. All pathogenic CNVs were confirmed in the five matched FFPE tumor samples from patients carrying known pathogenic germline mutations and we additionally identified somatic loss of the second allele in BRCA1/2. Furthermore we detected pathogenic BRCA1/2 variants in four the eleven FFPE samples from patients of whom no blood was available for analysis.

CONCLUSIONS

We demonstrate that an adjusted targeted capture-based enrichment protocol is superior to commonly applied multiplex PCR-based protocols for reliable BRCA1/2 variant detection, including CNV-detection, using FFPE tumor samples.

Use of FFPE-derived DNA in next generation sequencing: DNA extraction methods

Archival tissues represent a rich resource for clinical genomic studies, particularly when coupled with comprehensive medical records. Use of these in next generation sequencing (NGS) is a priority. Nine formalin-fixed paraffin-embedded (FFPE) DNA extraction methods were evaluated using twelve FFPE samples of varying tissue types. Quality assessment included total yield, percent dsDNA, fragment analysis and multiplex PCR. After assessment, three tissue types from four FFPE DNA methods were selected for NGS downstream evaluation, targeted and whole exome sequencing. In addition, two low input library protocols were evaluated for WES. Analysis revealed average coverage across the target regions for WES was ~20-30X for all four FFPE DNA extraction methods. For the targeted panels, the highest molecular tag coverage was obtained with the Kingfisher FFPE extraction method. The genotype concordance was 99% for the commonly called variant positions between all four extraction methods with the targeted PCR NGS panel and 96% with WES. Assessing quality of extracted DNA aids in selecting the optimal NGS approach, and the choice of both DNA extraction and library preparation approaches can impact the performance of archival tissue in NGS.

DNA extraction from FFPE tissue samples - a comparison of three procedures

Aim of the study

One of the critical steps in molecular oncology diagnostics is obtaining high quality genomic DNA. Therefore, it is important to evaluate and compare the techniques used to extract DNA from tissue samples. Since formalin-fixed, paraffin-embedded (FFPE) tissues are routinely used for both retrospective and prospective studies, we compared three commercially available methods of nucleic acid extraction in terms of quantity and quality of isolated DNA.

Material and methods

Slides prepared from 42 FFPE blocks were macro-dissected. Resulting material was divided and processed simultaneously using three extraction kits: QIAamp DNA FFPE Tissue Kit (QIAGEN), Cobas DNA Sample Preparation Kit (Roche Molecular Systems) and Maxwell 16 FFPE Plus LEV DNA Purification Kit (Promega). Subsequently, quantity and quality of obtained DNA samples were analysed spectrophotometrically (NanoDrop 2000, Thermo Scientific). Results of quantitative analysis were confirmed by a fluorometric procedure (Qubit 3.0 Fluorometer, Life Technologies).

Results

The results demonstrated that the yields of total DNA extracted using either Maxwell or Cobas methods were significantly higher compared to the QIAamp method (p < 0.001). The Maxwell Extraction Kit delivered DNA samples of the highest quality (p < 0.01). However, the highest total yield of extracted DNA was achieved with the Cobas technique, which may be due to a higher volume of eluate compared to the Maxwell method.

Conclusions

To our knowledge, this is the first paper which directly compares three extraction methods: Cobas, Maxwell and QIAamp. The data herein provide information required for the selection of a protocol that best suits the needs of the overall study design in terms of the quantity and quality of the extracted DNA.

Next-Generation Sequencing (NGS) Methods Show Superior or Equivalent Performance to Non-NGS Methods on BRAF, EGFR, and KRAS Proficiency Testing Samples

Context.—

There has been a rapid expansion of next-generation sequencing (NGS)–based assays for the detection of somatic variants in solid tumors. However, limited data are available regarding the comparative performance of NGS and non-NGS assays using standardized samples across a large number of laboratories.

Objective.—

To compare the performance of NGS and non-NGS assays using well-characterized proficiency testing samples provided by the College of American Pathologists (CAP) Molecular Oncology Committee. A secondary goal was to compare the use of preanalytic and postanalytic practices.

Design.—

A total of 17 343 responses were obtained from participants in the BRAF, EGFR, KRAS, and the Multigene Tumor Panel surveys across 84 different proficiency testing samples interrogating 16 variants and 3 wild-type sequences. Performance and preanalytic/postanalytic practices were analyzed by method.

Results.—

While both NGS and non-NGS achieved an acceptable response rate of greater than 95%, the overall performance of NGS methods was significantly better than that of non-NGS methods for the identification of variants in BRAF (overall 97.8% versus 95.6% acceptable responses, P = .001) and EGFR (overall 98.5% versus 97.3%, P = .01) and was similar for KRAS (overall 98.8% and 97.6%, P = .10). There were specific variant differences, but in all discrepant cases, NGS methods outperformed non-NGS methods. NGS laboratories also more consistently used preanalytic and postanalytic practices suggested by the CAP checklist requirements than non-NGS laboratories.

Conclusions.—

The overall analytic performance of both methods was excellent. For specific BRAF and EGFR variants, NGS outperformed non-NGS methods and NGS laboratories report superior adherence to suggested laboratory practices.

Mutational load in carotid body tumor

Background

Carotid body tumor (CBT) is a rare neoplasm arising from paraganglion located near the bifurcation of the carotid artery. There is great intra-tumor heterogeneity, and CBT development could be associated with both germline and somatic allelic variants. Studies on the molecular genetics of CBT are limited, and the molecular mechanisms of its pathogenesis are not fully understood. This work is focused on the estimation of mutational load (ML) in CBT.

Methods

Using the NextSeq 500 platform, we performed exome sequencing of tumors with matched lymph node tissues and peripheral blood obtained from six patients with CBT. To obtain reliable results in tumors with low ML, we developed and successfully applied a complex approach for the analysis of sequencing data. ML was evaluated as the number of somatic variants per megabase (Mb) of the target regions covered by the Illumina TruSeq Exome Library Prep Kit.

Results

The ML in CBT varied in the range of 0.09–0.28/Mb. Additionally, we identified several pathogenic/likely pathogenic somatic and germline allelic variants across six patients studied (including TP53 variants).

Conclusions

Using the developed approach, we estimated the ML in CBT, which is much lower than in common malignant tumors. Identified variants in known paraganglioma/pheochromocytoma-causative genes and novel genes could be associated with the pathogenesis of CBT. The obtained results expand our knowledge of the mutation process in CBT as well as the biology of tumor development.

The use of technical replication for detection of low-level somatic mutations in next-generation sequencing

Accurate genome-wide detection of somatic mutations with low variant allele frequency (VAF, <1%) has proven difficult, for which generalized, scalable methods are lacking. Herein, we describe a new computational method, called RePlow, that we developed to detect low-VAF somatic mutations based on simple, library-level replicates for next-generation sequencing on any platform. Through joint analysis of replicates, RePlow is able to remove prevailing background errors in next-generation sequencing analysis, facilitating remarkable improvement in the detection accuracy for low-VAF somatic mutations (up to ~99% reduction in false positives). The method is validated in independent cancer panel and brain tissue sequencing data. Our study suggests a new paradigm with which to exploit an overwhelming abundance of sequencing data for accurate variant detection.

Somatic mutations of low allele frequencies are often difficult to detect. Here, the authors develop RePlow, a computational method that leverages technical replication for detecting low-level somatic mutations using next-generation sequencing.

The prevalence of EBV and CMV DNA in epithelial ovarian cancer

Background

The underlying cause of epithelial ovarian cancer (EOC) is unknown. It has been theorized that infectious agents could contribute to ovarian tumorigenesis.

Objective

To investigate the potential role of oncogenic viral infection in EOC, we examined the prevalence of Epstein-Barr Virus (EBV) DNA and cytomegalovirus (CMV) DNA in EOC tissue samples.

Methods

Formalin-fixed, paraffin-imbedded (FFPE) tumor tissue samples from 198 patients included in the Danish Pelvic Mass Study were studied: 163 with serous adenocarcinomas, 15 with endometrioid adenocarcinomas, 11 with mucinous adenocarcinomas, and nine with clear-cell carcinomas. For controls in the EBV analysis, we used 176 tissue samples from patients diagnosed with benign mucinous cystadenomas. EBV and CMV genotyping was performed by real-time polymerase chain reaction with CMV and EBV CE-IVD approved kits. In-situ hybridization (ISH) was performed on the EBV positive samples.

Results

Sufficient DNA material was obtained in 191 and 174 tissue samples from cases and controls, respectively. Ten of 191 case samples (5.2%) and one of 174 control samples (0.5%) were positive for EBV DNA (P value = 0.011). CMV DNA was detected in only one case sample (0.5%). ISH confirmed that three of the samples were of stromal origin, while the remaining seven tested negative for EBV.

Conclusions

This study is the first to demonstrate a higher prevalence of EBV DNA in tissue samples from patients with EOC than in a benign control group. However, the cellular origin of seven of the samples could not be determined by ISH analysis. Our study did not support an association between CMV and EOC.

Cytological preparations for molecular analysis: A review of technical procedures, advantages and limitations for referring samples for testing

Minimally invasive procedures such as endobronchial ultrasound-guided transbronchial needle aspiration (EBUS-TBNA) must yield not only good quality and quantity of material for morphological assessment, but also an adequate sample for analysis of molecular markers to guide patients to appropriate targeted therapies. In this context, cytopathologists worldwide should be familiar with minimum requirements for refereeing cytological samples for testing. The present manuscript is a review with comprehensive description of the content of the workshop entitled Cytological preparations for molecular analysis: pre-analytical issues for EBUS TBNA, presented at the 40th European Congress of Cytopathology in Liverpool, UK. The present review emphasises the advantages and limitations of different types of cytology substrates used for molecular analysis such as archival smears, liquid-based preparations, archival cytospin preparations and FTA (Flinders Technology Associates) cards, as well as their technical requirements/features. These various types of cytological specimens can be successfully used for an extensive array of molecular studies, but the quality and quantity of extracted nucleic acids rely directly on adequate pre-analytical assessment of those samples. In this setting, cytopathologists must not only be familiar with the different types of specimens and associated technical procedures, but also correctly handle the material provided by minimally invasive procedures, ensuring that there is sufficient amount of material for a precise diagnosis and correct management of the patient through personalised care.

Next generation sequencing-based gene panel tests for the management of solid tumors

Next generation sequencing (NGS) has been an invaluable tool to put genomic sequencing into clinical practice. The incorporation of clinically relevant target sequences into NGS-based gene panel tests has generated practical diagnostic tools that enable individualized cancer-patient care. The clinical utility of gene panel testing includes investigation of the genetic basis for an individual's response to therapy, such as signaling pathways associated with a response to specific therapies, microsatellite instability and a hypermutated phenotype, and deficiency in the DNA double-strand break repair pathway. In this review, we describe the concept of precision cancer medicine using target sequences in gene panel tests as well as the importance of the control of sample quality in routine NGS-based genomic testing. We describe geographic and ethnic differences in cancer genomes, and discuss issues that need to be addressed in the future based on our experiences in Japan.

Current views on tumor mutational burden in patients with non-small cell lung cancer treated by immune checkpoint inhibitors

In the last few years, the treatment of patients with non-small cell lung cancer (NSCLC) has impressively benefitted from immunotherapy, in particular from the inhibition of immune checkpoints such as programmed cell death-1 (PD-1) and programmed cell death ligand-1 (PD-L1). However, despite the significant survival benefit for some patients with advanced NSCLC, the objective response rates (ORRs) remain relatively low no more than 20-30% with a large proportion of patients demonstrating primary resistance. Although the selection of NSCLC patients for the first-line treatment is currently guided by the expression of PD-L1 in tumor cells as detected by immunohistochemistry, this is not the case for the second-line setting. Moreover, the sensitivity and specificity of PD-L1 expression is modest which has prompted the search for additional predictive biomarkers. In this context, the assessment of the tumor mutational burden (TMB), defined as the total number of nonsynonymous mutations in the coding regions of genes, has recently emerged as an additional powerful biomarker to select patients for immunotherapy. The purpose of our review is to highlight the recent advances as well as the challenges and perspectives in the field of TMB and immunotherapy for patients with NSCLC.

The Potential Clinical Utility of Circulating Tumor DNA in Esophageal Adenocarcinoma: From Early Detection to Therapy

Esophageal adenocarcinoma (EAC) is a lethal cancer requiring improved screening strategies and treatment options due to poor detection methods, aggressive progression, and therapeutic resistance. Emerging circulating tumor DNA (ctDNA) technologies may offer a unique non-invasive strategy to better characterize the highly heterogeneous cancer and more clearly establish the genetic modulations leading to disease progression. The presented review describes the potential advantages of ctDNA methodologies as compared to current clinical strategies to improve clinical detection, enhance disease surveillance, evaluate prognosis, and personalize treatment. Specifically, we describe the ctDNA-targetable genetic markers of prognostic significance to stratify patients into risk of progression from benign to malignant disease and potentially offer cost-effective screening of established cancer. We also describe the application of ctDNA to more effectively characterize the heterogeneity and particular mutagenic resistance mechanisms in real-time to improve prognosis and therapeutic monitoring strategies. Lastly, we discuss the inconsistent clinical responses to currently approved therapies for EAC and the role of ctDNA to explore the dynamic regulation of novel targeted and immunotherapies to personalize therapy and improve patient outcomes. Although there are clear limitations of ctDNA technologies for immediate clinical deployment, this review presents the prospective role of such applications to potentially overcome many of the notable hurdles to treating EAC patients. A deeper understanding of complex EAC tumor biology may result in the progress toward improved clinical outcomes.

A novel approach to stabilize fetal cell-free DNA fraction in maternal blood samples for extended period of time

This study was undertaken to evaluate a novel method for stabilizing and preserving the original proportion of cell-free fetal DNA (cffDNA) in maternal blood for extended periods of time without using crosslinking agents, such as formaldehyde, which compromise DNA integrity and extraction efficiency. Blood was drawn from pregnant donors into K₃EDTA and Blood Exo DNA ProTeck® (ProTeck) tubes. Blood drawn into both tubes were aliquoted and stored at three different temperatures. At indicated times sample aliquots were processed for cell-free DNA (cfDNA) extraction. Plasma cfDNA and cffDNA quantified by droplet digital PCR (ddPCR) assay which amplify RASSF1A gene promoter region. ProTeck reagent is formaldehyde free and inhibits blood cell metabolism in blood samples during storage. Cell-free DNA concentration increased over time in blood plasma stored in K₃EDTA tubes at 4, 22 and 30°C. Blood stored in ProTeck tubes, cfDNA concentration was stable at 4, 22 and 30°C for 21, 28 and 7 days, respectively. In K₃EDTA tubes cffDNA proportion decreases steadily over time whereas in ProTeck tubes cffDNA proportion remained stable. This novel technology stabilizes cffDNA proportion in maternal blood samples at 4, 22 and 30°C for 21, 28 and 7 days, respectively.

The evolving role of receptors as predictive biomarkers for metastatic breast cancer

INTRODUCTION

In breast cancer, estrogen receptor (ER) and human epidermal growth factor receptor 2 (HER2) are essential biomarkers to predict response to endocrine and anti-HER2 therapies, respectively. In metastatic breast cancer, the use of these receptors and targeted therapies present additional challenges: temporal heterogeneity, together with limited sampling methodologies, hinders receptor status assessment, and the constant evolution of the disease invariably leads to resistance to treatment. Areas covered: This review summarizes the genomic abnormalities in ER and HER2, such as mutations, amplifications, translocations, and alternative splicing, emerging as novel biomarkers that provide an insight into underlying mechanisms of resistance and hold potential predictive value to inform treatment selection. We also describe how liquid biopsies for sampling of circulating markers and ultrasensitive detection technologies have emerged which complement ongoing efforts for biomarker discovery and analysis. Expert commentary: While evidence suggests that genomic aberrations in ER and HER2 could contribute to meeting the pressing need for better predictive biomarkers, efforts need to be made to standardize assessment methods and better understand the resistance mechanisms these markers denote. Taking advantage of emerging technologies, research in upcoming years should include prospective trials incorporating these predictors into the study design to validate their potential clinical value.

Methods of measurement for tumor mutational burden in tumor tissue

Immunotherapies based on immune checkpoint inhibitors are emerging as an innovative treatment for different types of advanced cancers. While the utility of immune checkpoint inhibitors has been clearly demonstrated, the response rate is highly variable across individuals. Due to the cost and toxicity of these immunotherapies, a critical challenge in this field is the identification of predictive biomarkers to discriminate which patients may respond to immunotherapy. Recently, a high tumor mutational burden (TMB) has been identified as a genetic signature that is associated with a favorable outcome for immune checkpoint inhibitor therapy. The TMB is defined as the total number of nonsynonymous mutations per coding area of a tumor genome. Initially, it was determined using whole exome sequencing, but due to the high costs and long turnaround time of this method, targeted panel sequencing is currently being explored to measure TMB. In the near future, TMB evaluation may play an important role in immuno-oncology, but its implementation in a routine setting involves robust analytical and clinical validation. Standardization is also needed in order to make informed decisions about patients. This review presents the methodologies employed for determining TMB and discusses the factors that may have an impact on its measurement.

The genomics of desmoplastic small round cell tumor reveals the deregulation of genes related to DNA damage response, epithelial-mesenchymal transition, and immune response

BACKGROUND

Desmoplastic small round cell tumor (DSRCT) is a rare, aggressive, and poorly investigated simple sarcoma with a low frequency of genetic deregulation other than an Ewing sarcoma RNA binding protein 1 (EWSR1)-Wilm's tumor suppressor (WT1) translocation. We used whole-exome sequencing to interrogate six consecutive pre-treated DSRCTs whose gene expression was previously investigated.

METHODS

DNA libraries were prepared from formalin-fixed, paraffin-embedded archival tissue specimens following the Agilent SureSelectXT2 target enrichment protocol and sequenced on Illumina NextSeq 500. Raw sequence data were aligned to the reference genome with Burrows-Wheeler Aligner algorithm. Somatic mutations and copy number alterations (CNAs) were identified using MuTect2 and EXCAVATOR2, respectively. Biological functions associated with altered genes were investigated through Ingenuity Pathway Analysis (IPA) software.

RESULTS

A total of 137 unique somatic mutations were identified: 133 mutated genes were case-specific, and 2 were mutated in two cases but in different positions. Among the 135 mutated genes, 27% were related to two biological categories: DNA damage-response (DDR) network that was also identified through IPA and mesenchymal-epithelial reverse transition (MErT)/epithelial-mesenchymal transition (EMT) already demonstrated to be relevant in DSRCT. The mutated genes in the DDR network were involved in various steps of transcription and particularly affected pre-mRNA. Half of these genes encoded RNA-binding proteins or DNA/RNA-binding proteins, which were recently recognized as a new class of DDR players. CNAs in genes/gene families, involved in MErT/EMT and DDR, were recurrent across patients and mostly segregated in the MErT/EMT category. In addition, recurrent gains of regions in chromosome 1 involving many MErT/EMT gene families and loss of one arm or the entire chromosome 6 affecting relevant immune-regulatory genes were recorded.

CONCLUSIONS

The emerging picture is an extreme inter-tumor heterogeneity, characterized by the concurrent deregulation of the DDR and MErT/EMT dynamic and plastic programs that could favour genomic instability and explain the refractory DSRCT profile.

Mutational Evolution in Relapsed Diffuse Large B-Cell Lymphoma

Current genomic models in diffuse large B-cell lymphoma (DLBCL) are based on single tumor biopsies, which might underestimate heterogeneity. Data on mutational evolution largely remains unknown. An exploratory study using whole exome sequencing on paired (primary and relapse) formalin fixed paraffin embedded DLBCL biopsies (n = 14) of 6 patients was performed to globally assess the mutational evolution and to identify gene mutations specific for relapse samples from patients treated with rituximab, cyclophosphamide, doxorubicin, vincristine, and prednisolone. A minority of the mutations detected in the primary sample (median 7.6%, range 4.8–66.2%) could not be detected in the matching relapse sample. Relapsed DLBCL samples showed a mild increase of mutations (median 12.5%, range 9.4–87.6%) as compared to primary tumor biopsies. We identified 264 genes possibly related to therapy resistance, including tyrosine kinases (n = 18), (transmembrane) glycoproteins (n = 73), and genes involved in the JAK-STAT pathway (n = 7). Among the potentially resistance related genes were PIM1, SOCS1, and MYC, which have been reported to convey a risk for treatment failure. In conclusion, we show modest temporal heterogeneity between paired tumor samples with the acquisition of new mutations and identification of genes possibly related to therapy resistance. The mutational evolution could have implications for treatment decisions and development of novel targeted drugs.

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.

Relationship between formalin reagent and success rate of targeted sequencing analysis using formalin fixed paraffin embedded tissues

BACKGROUND

Tumor genetic alterations are determined to aid in selecting therapy and predicting prognosis. In routine clinical practice, targeted sequencing analysis is performed using formalin-fixed paraffin embedded (FFPE) tissues. However, successful genetic analysis remains challenging because FFPE DNA is fragmented during the sample preparation process.

METHODS

Real-time PCR was performed to assess DNA quality and quantities. Targeted sequencing was performed using FFPE tissues fixed with different types of formalin.

RESULTS

DNA was less fragmented from samples fixed in low formalin concentration (10% vs. 20%) and neutral buffered conditions (neutral buffered vs. non-neutral). DNA fragmentation increased over the fixation time. In a preliminary test study, we compared fixation using 10% neutral buffered formalin (n = 180) and 20% formalin (n = 26). The success rate of targeted analysis was higher using 10% neutral formalin (98.3%; 177/180) compared with 20% formalin (34.6%; 9/26). In a validation study with additional formalin-fixed paraffin embedded tissues fixed with 10% neutral buffered formalin (n = 860), we reproduced these results and achieved a high success rate for targeted sequencing analysis (98.4%; 846/860).

CONCLUSION

Our data show that 10% neutral buffered formalin is recommended for fixation of formalin-fixed paraffin embedded samples to achieve high success rate of targeted sequencing analysis.

Association analysis using somatic mutations

Somatic mutations drive the growth of tumor cells and are pivotal biomarkers for many cancer treatments. Genetic association analysis using somatic mutations is an effective approach to study the functional impact of somatic mutations. However, standard regression methods are not appropriate for somatic mutation association studies because somatic mutation calls often have non-ignorable false positive rate and/or false negative rate. While large scale association analysis using somatic mutations becomes feasible recently—thanks for the improvement of sequencing techniques and the reduction of sequencing cost—there is an urgent need for a new statistical method designed for somatic mutation association analysis. We propose such a method with computationally efficient software implementation: Somatic mutation Association test with Measurement Errors (SAME). SAME accounts for somatic mutation calling uncertainty using a likelihood based approach. It can be used to assess the associations between continuous/dichotomous outcomes and individual mutations or gene-level mutations. Through simulation studies across a wide range of realistic scenarios, we show that SAME can significantly improve statistical power than the naive generalized linear model that ignores mutation calling uncertainty. Finally, using the data collected from The Cancer Genome Atlas (TCGA) project, we apply SAME to study the associations between somatic mutations and gene expression in 12 cancer types, as well as the associations between somatic mutations and colon cancer subtype defined by DNA methylation data. SAME recovered some interesting findings that were missed by the generalized linear model. In addition, we demonstrated that mutation-level and gene-level analyses are often more appropriate for oncogene and tumor-suppressor gene, respectively.

Cancer is a genetic disease that is driven by the accumulation of somatic mutations. Association studies using somatic mutations is a powerful approach to identify the potential impact of somatic mutations on molecular or clinical features. One challenge for such tasks is the non-ignorable somatic mutation calling errors. We have developed a statistical method to address this challenge and applied our method to study the gene expression traits associated with somatic mutations in 12 cancer types. Our results show that some somatic mutations affect gene expression in several cancer types. In particular, we show that the associations between gene expression traits and TP53 gene level mutation reveal some similarities across a few cancer types.

Solid-phase enzyme catalysis of DNA end repair and 3' A-tailing reduces GC-bias in next-generation sequencing of human genomic DNA

The use of next-generation sequencing (NGS) has been instrumental in advancing biological research and clinical diagnostics. To fully utilize the power of NGS, complete, uniform coverage of the entire genome is required. In this study, we identified the primary sources of bias observed in sequence coverage across AT-rich regions of the human genome with existing amplification-free DNA library preparation methods. We have found evidence that a major source of bias is the inefficient processing of AT-rich DNA in end repair and 3' A-tailing, causing under-representation of extremely AT-rich regions. We have employed immobilized DNA modifying enzymes to catalyze end repair and 3' A-tailing reactions, to notably reduce the GC bias observed with existing library construction methods.

Reducing DNA damage by formaldehyde in liquid-based cytology preservation solutions to enable the molecular testing of lung cancer specimens

BACKGROUND

Liquid-based cytology (LBC) is a useful cytopathological method, and LBC lung adenocarcinoma specimens may be used for genetic analysis in the near future. In the current study, the authors determined whether LBC specimens can be used for epidermal growth factor receptor (EGFR) mutation analysis in human lung adenocarcinoma cell lines.

METHODS

Genomic DNA was extracted from 3 lung adenocarcinoma cell lines that were fixed in LBC preservation solution using 2 protocols (one for cultured cells and one for tissues) of a DNA extraction kit. Different fixation times were tested for each protocol: 30 minutes, 1 hour, and 1 to 9 days. As controls, cells also were fixed in 10% formalin or 95% ethanol. The authors investigated the effect of fixation time on DNA fragmentation, polymerase chain reaction (PCR) amplification, and EGFR mutation detection.

RESULTS

The DNA yield of LBC specimens tended to decrease depending on fixation time. When using the DNA extraction protocol for tissues, PCR amplification was successful after 9 days of fixation, although extracted genomic DNA that was fixed for >1 hour demonstrated fragmentation. Mutation analyses using the Cycleave PCR method were successful after 7 days of fixation. The DNA extraction protocol for tissues was appropriate for lung adenocarcinoma cell lines that were stored for >1 day in a preservative solution. The results of the current study demonstrated that EGFR mutations can be detected on day 7 using lung adenocarcinoma cell lines fixed in CytoRich Red preservative.

CONCLUSIONS

When LBC specimens are used for targeted molecular genetic testing, the appropriate preservative solution and extraction protocol first should be determined.

Modified SureSelectQXT Target Enrichment Protocol for Illumina Multiplexed Sequencing of FFPE Samples

Background

Personalised medicine is nowadays a major objective in oncology. Molecular characterization of tumours through NGS offers the possibility to find possible therapeutic targets in a time- and cost-effective way. However, the low quality and complexity of FFPE DNA samples bring a series of disadvantages for massive parallel sequencing techniques compared to high-quality DNA samples (from blood cells, cell cultures, etc.).

Results

We performed several experiments to understand the behaviour of FFPE DNA samples during the construction of SureSelect^QXT libraries. First, we designed a quality checkpoint for FFPE DNA samples based on the quantification of their amplification capability (qcPCR). We observed that FFPE DNA samples can be classified according to DIN value and qcPCR concentration into unusable, or low-quality (LQ) and good-quality (GQ) DNA. For GQ samples, we increased the amount of input DNA to 150 ng and the digestion time to 30 min, whereas for LQ samples, we used 50 ng of DNA as input but we decreased the digestion time to 1 min. In all cases, we increased the cycles of the pre-hyb PCR to 10 but decreased the cycles of the post-hyb PCR to 8. In addition, we confirmed that using half of the volume of reagents can be beneficial. Finally, in order to obtain better results, we designed a decision flow-chart to achieve a seeding concentration of 12–14 pM for MiSeq Reagent Kit v2.

Conclusions

Our experiments allowed us to unveil the behaviour of low-quality FFPE DNA samples during the construction of SureSelect^QXT libraries. Sequencing results showed that, using our modified SureSelect^QXT protocol, the final percentage of usable reads for low-quality samples was increased more than three times allowing to reach median depth/million reads values of 76.35. This value is equivalent to ~ 0.9 and ~ 0.7 of the values obtained for good-quality FFPE and high-quality DNA respectively.

Electronic supplementary material

The online version of this article (10.1186/s12575-018-0084-7) contains supplementary material, which is available to authorized users.

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.