Deciphering the role of FUS::DDIT3 expression and tumor microenvironment in myxoid liposarcoma development.

BACKGROUND: Myxoid liposarcoma (MLS) displays a distinctive tumor microenvironment and is characterized by the FUS::DDIT3 fusion oncogene, however, the precise functional contributions of these two elements remain enigmatic in tumor development. METHODS: To study the cell-free microenvironment in MLS, we developed an experimental model system based on decellularized patient-derived xenograft tumors. We characterized the cell-free scaffold using mass spectrometry. Subsequently, scaffolds were repopulated using sarcoma cells with or without FUS::DDIT3 expression that were analyzed with histology and RNA sequencing. RESULTS: Characterization of cell-free MLS scaffolds revealed intact structure and a large variation of protein types remaining after decellularization. We demonstrated an optimal culture time of 3 weeks and showed that FUS::DDIT3 expression decreased cell proliferation and scaffold invasiveness. The cell-free MLS microenvironment and FUS::DDIT3 expression both induced biological processes related to cell-to-cell and cell-to-extracellular matrix interactions, as well as chromatin remodeling, immune response, and metabolism. Data indicated that FUS::DDIT3 expression more than the microenvironment determined the pre-adipocytic phenotype that is typical for MLS. CONCLUSIONS: Our experimental approach opens new means to study the tumor microenvironment in detail and our findings suggest that FUS::DDIT3-expressing tumor cells can create their own extracellular niche.

Tissue specificity of oncogenic BRAF targeted to lung and thyroid through a shared lineage factor.

Cells of origin in cancer determine tumor phenotypes, but whether lineage-defining transcription factors might influence tissue specificity of tumorigenesis among organs with similar developmental traits are unknown. We demonstrate here that tumor development and progression markedly differ in lung and thyroid targeted by Braf mutation in Nkx2.1CreERT2 mice heterozygous for Nkx2-1. In absence of tamoxifen, non-induced Nkx2.1CreERT2;BrafCA/+ mutants developed multiple full-blown lung adenocarcinomas with a latency of 1-3 months whereas thyroid tumors were rare and constrained, although minute BrafCA activation documented by variant allele sequencing was similar in both tissues. Induced oncogene activation accelerated neoplastic growth only in the lungs. By contrast, NKX2-1+ progenitor cells were equally responsive to constitutive expression of mutant Braf during lung and thyroid development. Both lung and thyroid cells transiently downregulated NKX2-1 in early tumor stages. These results indicate that BRAFV600E-induced tumorigenesis obey organ-specific traits that might be differentially modified by a shared lineage factor.

Extracellular vesicle DNA from human melanoma tissues contains cancer-specific mutations.

Liquid biopsies are promising tools for early diagnosis and residual disease monitoring in patients with cancer, and circulating tumor DNA isolated from plasma has been extensively studied as it has been shown to contain tumor-specific mutations. Extracellular vesicles (EVs) present in tumor tissues carry tumor-derived molecules such as proteins and nucleic acids, and thus EVs can potentially represent a source of cancer-specific DNA. Here we identified the presence of tumor-specific DNA mutations in EVs isolated from six human melanoma metastatic tissues and compared the results with tumor tissue DNA and plasma DNA. Tumor tissue EVs were isolated using enzymatic treatment followed by ultracentrifugation and iodixanol density cushion isolation. A panel of 34 melanoma-related genes was investigated using ultra-sensitive sequencing (SiMSen-seq). We detected mutations in six genes in the EVs (BRAF, NRAS, CDKN2A, STK19, PPP6C, and RAC), and at least one mutation was detected in all melanoma EV samples. Interestingly, the mutant allele frequency was higher in DNA isolated from tumor-derived EVs compared to total DNA extracted directly from plasma DNA, supporting the potential role of tumor EVs as future biomarkers in melanoma.

UMIErrorCorrect and UMIAnalyzer: Software for Consensus Read Generation, Error Correction, and Visualization Using Unique Molecular Identifiers.

BACKGROUND: Targeted sequencing using unique molecular identifiers (UMIs) enables detection of rare variant alleles in challenging applications, such as cell-free DNA analysis from liquid biopsies. Standard bioinformatics pipelines for data processing and variant calling are not adapted for deep-sequencing data containing UMIs, are inflexible, and require multistep workflows or dedicated computing resources. METHODS: We developed a bioinformatics pipeline using Python and an R package for data analysis and visualization. To validate our pipeline, we analyzed cell-free DNA reference material with known mutant allele frequencies (0%, 0.125%, 0.25%, and 1%) and public data sets. RESULTS: We developed UMIErrorCorrect, a bioinformatics pipeline for analyzing sequencing data containing UMIs. UMIErrorCorrect only requires fastq files as inputs and performs alignment, UMI clustering, error correction, and variant calling. We also provide UMIAnalyzer, a graphical user interface, for data mining, visualization, variant interpretation, and report generation. UMIAnalyzer allows the user to adjust analysis parameters and study their effect on variant calling. We demonstrated the flexibility of UMIErrorCorrect by analyzing data from 4 different targeted sequencing protocols. We also show its ability to detect different mutant allele frequencies in standardized cell-free DNA reference material. UMIErrorCorrect outperformed existing pipelines for targeted UMI sequencing data in terms of variant detection sensitivity. CONCLUSIONS: UMIErrorCorrect and UMIAnalyzer are comprehensive and customizable bioinformatics tools that can be applied to any type of library preparation protocol and enrichment chemistry using UMIs. Access to simple, generic, and open-source bioinformatics tools will facilitate the implementation of UMI-based sequencing approaches in basic research and clinical applications.

Digital Quantification of Chemical Oligonucleotide Synthesis Errors.

BACKGROUND: Chemically synthesized oligonucleotides are vital to most nucleic acids-based technologies and several applications are sensitive to oligonucleotide sequence errors. However, it is challenging to identify and quantify the types and amount of errors in synthetic oligonucleotides. METHODS: We applied a digital sequencing approach using unique molecular identifiers to quantify errors in chemically synthesized oligonucleotides from multiple manufacturers with different synthesis strategies, purity grades, batches, and sequence context. RESULTS: We detected both deletions and substitutions in chemical oligonucleotide synthesis, but deletions were 7 times more common. We found that 97.2% of all analyzed oligonucleotide molecules were intact across all manufacturers and purity grades, although the number of oligonucleotide molecules with deletions ranged between 0.2% and 11.7% for different types. Different batches of otherwise identical oligonucleotide types also varied significantly, and batch effect can impact oligonucleotide quality more than purification. We observed a bias of increased deletion rates in chemically synthesized oligonucleotides toward the 5'-end for 1 out of 2 sequence configurations. We also demonstrated that the performance of sequencing assays depends on oligonucleotide quality. CONCLUSIONS: Our data demonstrate that manufacturer, synthesis strategy, purity, batch, and sequence context all contribute to errors in chemically synthesized oligonucleotides and need to be considered when choosing and evaluating oligonucleotides. High-performance oligonucleotides are essential in numerous molecular applications, including clinical diagnostics.

The PEMDAC phase 2 study of pembrolizumab and entinostat in patients with metastatic uveal melanoma.

Preclinical studies have suggested that epigenetic therapy could enhance immunogenicity of cancer cells. We report the results of the PEMDAC phase 2 clinical trial (n = 29; NCT02697630) where the HDAC inhibitor entinostat was combined with the PD-1 inhibitor pembrolizumab in patients with metastatic uveal melanoma (UM). The primary endpoint was objective response rate (ORR), and was met with an ORR of 14%. The clinical benefit rate at 18 weeks was 28%, median progression free survival was 2.1 months and the median overall survival was 13.4 months. Toxicities were manageable, and there were no treatment-related deaths. Objective responses and/or prolonged survival were seen in patients with BAP1 wildtype tumors, and in one patient with an iris melanoma that exhibited a UV signature. Longer survival also correlated with low baseline ctDNA levels or LDH. In conclusion, HDAC inhibition and anti-PD1 immunotherapy results in durable responses in a subset of patients with metastatic UM.Trial registration ClinicalTrials.gov registration number: NCT02697630 (registered 3 March 2016). EudraCT registration number: 2016-002114-50.

Response to BRAF/MEK Inhibition in A598_T599insV BRAF Mutated Melanoma.

Approximately 50% of patients with metastatic melanoma harbor an activating BRAF mutation. Tumors with activating mutation BRAF gene proliferate excessively and can be treated with targeted BRAF-inhibitors in combination with MEK inhibitors. The most common BRAF mutation occurs at amino acid position 600. Other BRAF mutations are rare and their predictive value for treatment response to BRAF/MEK inhibition is low. Here, we report on a patient with a BRAF A598_T599insV mutated melanoma, a mutation that has only been described in one previous melanoma patient in which the treatment response to BRAF/MEK inhibition was transient. Our patient had a large ulcerated metastasis that showed a durable complete response implying that BRAF/MEK inhibition should be considered a treatment option for this mutation. We analyzed circulating cell-free tumor DNA (ctDNA) carrying the BRAF A598_T599insV mutation throughout treatment. The allele frequency of BRAF A598_T599insV decreased during regression of the tumors, indicating that this method has potential to monitor treatment response. Our case demonstrates durable response to BRAF/MEK inhibition in a melanoma patient carrying a BRAF A598_T599insV mutation. In addition, we show that allele frequency analysis of A598_T599insV mutation in blood using ultrasensitive sequencing can be used to monitor treatment response.

Corrigendum to "Considerations and quality controls when analyzing cell-free tumor DNA" [Biomol. Detect. Quantif. 17 (2019) 100078].

[This corrects the article DOI: 10.1016/j.bdq.2018.12.003.].

Detection of Circulating Tumor DNA in Plasma: A Potential Biomarker for Esophageal Adenocarcinoma.

BACKGROUND: Recent literature has demonstrated the potential of "liquid biopsy" and detection of circulating tumor (ct)DNA as a cancer biomarker. However, to date there is a lack of data specific to esophageal adenocarcinoma (EAC). This study was conducted to determine how detection and quantification of ctDNA changes with disease burden in patients with EAC and evaluate its potential as a biomarker in this population. METHODS: Blood samples were obtained from patients with stage I to IV EAC. Longitudinal blood samples were collected from a subset of patients. Imaging studies and pathology reports were reviewed to determine disease course. Tumor samples were sequenced to identify mutations. Mutations in plasma DNA were detected using custom, barcoded, patient-specific sequencing libraries. Mutations in plasma were quantified, and associations with disease stage and response to therapy were explored. RESULTS: Plasma samples from a final cohort of 38 patients were evaluated. Baseline plasma samples were ctDNA positive for 18 patients (47%) overall, with tumor allele frequencies ranging from 0.05% to 5.30%. Detection frequency of ctDNA and quantity of ctDNA increased with stage. Data from longitudinal samples indicate that ctDNA levels correlate with and precede evidence of response to therapy or recurrence. CONCLUSIONS: ctDNA can be detected in plasma of EAC patients and correlates with disease burden. Detection of ctDNA in early-stage EAC is challenging and may limit diagnostic applications. However, our data demonstrate the potential of ctDNA as a dynamic biomarker to monitor treatment response and disease recurrence in patients with EAC.

Multilaboratory Assessment of a New Reference Material for Quality Assurance of Cell-Free Tumor DNA Measurements.

We conducted a multilaboratory assessment to determine the suitability of a new commercially available reference material with 40 cancer variants in a background of wild-type DNA at four different variant allele frequencies (VAFs): 2%, 0.50%, 0.125%, and 0%. The variants include single nucleotides, insertions, deletions, and two structural variations selected for their clinical importance and to challenge the performance of next-generation sequencing (NGS) methods. Fragmented DNA was formulated to simulate the size distribution of circulating wild-type and tumor DNA in a synthetic plasma matrix. DNA was extracted from these samples and characterized with different methods and multiple laboratories. The various extraction methods had differences in yield, perhaps because of differences in chemistry. Digital PCR assays were used to measure VAFs to compare results from different NGS methods. Comparable VAFs were observed across the different NGS methods. This multilaboratory assessment demonstrates that the new reference material is an appropriate tool to determine the analytical parameters of different measurement methods and to ensure their quality assurance.

Impact of Polymerase Fidelity on Background Error Rates in Next-Generation Sequencing with Unique Molecular Identifiers/Barcodes.

Liquid biopsy and detection of tumor-associated mutations in cell-free circulating DNA often requires the ability to identify single nucleotide variants at allele frequencies below 0.1%. Standard sequencing protocols cannot achieve this level of sensitivity due to background noise from DNA damage and polymerase induced errors. Addition of unique molecular identifiers allows identification and removal of errors responsible for this background noise. Theoretically, high fidelity enzymes will also reduce error rates in barcoded NGS but this has not been thoroughly explored. We evaluated the impact of polymerase fidelity on the magnitude of error reduction at different steps of barcoded NGS library construction. We find that barcoding itself displays largest impact on error reduction, even with low fidelity polymerases. Use of high fidelity polymerases in the barcoding step of library construction further suppresses error in barcoded NGS, and allows detection of variant alleles below 0.1% allele frequency. However, the improvement in error correction is modest and is not directly proportional to polymerase fidelity. Depending on the specific application, other polymerase characteristics such as multiplexing capacity, PCR efficiency, buffer requirements and ability to amplify targets with high GC content may outweigh the relatively small additional decrease in error afforded by ultra-high fidelity polymerases.

Considerations and quality controls when analyzing cell-free tumor DNA.

Circulating cell-free tumor DNA (ctDNA) is a promising biomarker in cancer. Ultrasensitive technologies enable detection of low (< 0.1%) mutant allele frequencies, a pre-requisite to fully utilize the potential of ctDNA in cancer diagnostics. In addition, the entire liquid biopsy workflow needs to be carefully optimized to enable reliable ctDNA analysis. Here, we discuss important considerations for ctDNA detection in plasma. We show how each experimental step can easily be evaluated using simple quantitative PCR assays, including detection of cellular DNA contamination and PCR inhibition. Furthermore, ctDNA assay performance is also demonstrated to be affected by both DNA fragmentation and target sequence. Finally, we show that quantitative PCR is useful to estimate the required sequencing depth and to monitor DNA losses throughout the workflow. The use of quality control assays enables the development of robust and standardized workflows that facilitate the implementation of ctDNA analysis into clinical routine.

Plasma circulating tumor DNA as a potential tool for disease monitoring in head and neck cancer.

BACKGROUND: Recommendations for perioperative therapy in head and neck cancer are not explicit and recurrence occurs frequently. Circulating tumor DNA is an emerging cancer biomarker, but has not been extensively explored for detection of recurrence in head and neck cancer. METHODS: Patients diagnosed with head and neck squamous cell carcinoma were recruited into the study protocol. Tumors were sequenced to identify patient-specific mutations. Mutations were then identified in plasma circulating tumor DNA from pre-treatment blood samples and longitudinally during standard follow-up. Circulating tumor DNA status during follow-up was correlated to disease recurrence. RESULTS: Samples were taken from eight patients. Tumor mutations were verified in seven patients. Baseline circulating tumor DNA was positive in six patients. Recurrence occurred in four patients, two of whom had detectable circulating tumor DNA prior to recurrence. CONCLUSION: Circulating tumor DNA is a potential tool for disease and recurrence monitoring following curative therapy in head and neck cancer, allowing for better prognostication, and/or modification of treatment strategies.

Elevated pyrimidine dimer formation at distinct genomic bases underlies promoter mutation hotspots in UV-exposed cancers.

Sequencing of whole cancer genomes has revealed an abundance of recurrent mutations in gene-regulatory promoter regions, in particular in melanoma where strong mutation hotspots are observed adjacent to ETS-family transcription factor (TF) binding sites. While sometimes interpreted as functional driver events, these mutations are commonly believed to be due to locally inhibited DNA repair. Here, we first show that low-dose UV light induces mutations preferably at a known ETS promoter hotspot in cultured cells even in the absence of global or transcription-coupled nucleotide excision repair (NER). Further, by genome-wide mapping of cyclobutane pyrimidine dimers (CPDs) shortly after UV exposure and thus before DNA repair, we find that ETS-related mutation hotspots exhibit strong increases in CPD formation efficacy in a manner consistent with tumor mutation data at the single-base level. Analysis of a large whole genome cohort illustrates the widespread contribution of this effect to recurrent mutations in melanoma. While inhibited NER underlies a general increase in somatic mutation burden in regulatory elements including ETS sites, our data supports that elevated DNA damage formation at specific genomic bases is at the core of the prominent promoter mutation hotspots seen in skin cancers, thus explaining a key phenomenon in whole-genome cancer analyses.

Recurrent promoter mutations in melanoma are defined by an extended context-specific mutational signature.

Sequencing of whole tumor genomes holds the promise of revealing functional somatic regulatory mutations, such as those described in the TERT promoter. Recurrent promoter mutations have been identified in many additional genes and appear to be particularly common in melanoma, but convincing functional data such as influence on gene expression has been more elusive. Here, we show that frequently recurring promoter mutations in melanoma occur almost exclusively at cytosines flanked by a distinct sequence signature, TTCCG, with TERT as a notable exception. In active, but not inactive, promoters, mutation frequencies for cytosines at the 5' end of this ETS-like motif were considerably higher than expected based on a UV trinucleotide mutational signature. Additional analyses solidify this pattern as an extended context-specific mutational signature that mediates an exceptional position-specific vulnerability to UV mutagenesis, arguing against positive selection. We further use ultra-sensitive amplicon sequencing to demonstrate that cell cultures exposed to UV light quickly develop subclonal mutations specifically in affected positions. Our findings have implications for the interpretation of somatic mutations in regulatory regions, and underscore the importance of genomic context and extended sequence patterns to accurately describe mutational signatures in cancer.

Simple multiplexed PCR-based barcoding of DNA for ultrasensitive mutation detection by next-generation sequencing.

Detection of extremely rare variant alleles within a complex mixture of DNA molecules is becoming increasingly relevant in many areas of clinical and basic research, such as the detection of circulating tumor DNA in the plasma of cancer patients. Barcoding of DNA template molecules early in next-generation sequencing (NGS) library construction provides a way to identify and bioinformatically remove polymerase errors that otherwise make detection of these rare variants very difficult. Several barcoding strategies have been reported, but all require long and complex library preparation protocols. Simple, multiplexed, PCR-based barcoding of DNA for sensitive mutation detection using sequencing (SiMSen-seq) was developed to generate targeted barcoded libraries with minimal DNA input, flexible target selection and a very simple, short (∼4 h) library construction protocol. The protocol comprises a three-cycle barcoding PCR step followed directly by adaptor PCR to generate the library and then bead purification before sequencing. Thus, SiMSen-seq allows detection of variant alleles at <0.1% frequency with easy customization of library content (from 1 to 40+ PCR amplicons) and a protocol that can be implemented in any molecular biology laboratory. Here, we provide a detailed protocol for assay development and describe software to process the barcoded sequence reads.

Modelling and analysis of a gene-regulatory feed-forward loop with basal expression of the second regulator.

Efficient adaptation strategies to changing environmental conditions are essential for bacteria to survive and grow. Fundamental restructuring of their metabolism is usually mediated by corresponding gene regulation. Here, often several different environmental stimuli have to be integrated into a reasonable, energy-efficient response. Fast fluctuations and overshooting have to be filtered out. The gene regulatory network for the anaerobic adaptation of the pathogenic bacterium Pseudomonas aeruginosa is organized as a feed-forward loop (FFL), which is a three-gene network motif composed of two transcription factors (Anr for oxygen, NarxL for nitrate) and one target (Nar for nitrate reductase). The upstream transcription factor (Anr) induces the downstream transcription factor (NarXL). Both regulators act together positively by inducing the target (Nar) via a direct and indirect regulation path (coherent type-1 FFL). Since full promoter activity is only achieved when both transcription factors are present the target operon is expressed with a delay. Thus, in response to environmental stimuli (oxygen, nitrate), signals are mediated and processed in a way that short pulses are filtered out. In this study we analyze a special kind of FFL called FFLk by means of a family of ordinary differential equation models. The secondary FFL regulator (NarXL) is expressed constitutively but further induced in the presence of the upstream stimuli. This FFL modification has substantial influence on the response time and cost-benefit ratio mediated by environmental fluctuations. In order to find conditions where this regulatory network motif might be beneficial, we analyzed various models and environments. We describe the observed evolutional advantage of FFLk and its role in environmental adaptation and pathogenicity.