Abstract:Background: The somatic mutations that drive acute myeloid leukaemia (AML) are highly heterogenous. Identifying the mutations that drive individual cases can help determine patient prognosis and therapy. For this reason, genetic testing for prognostically important mutations in leukaemic DNA is now routine in many diagnostic laboratories. Also, the analysis of gene-expression profiles from AML RNA can provide additional clinically useful information that cannot be inferred from DNA sequencing. Whilst it is both expensive and impractical to carry out both types of sequencing and analyses, we hypothesised that the two could be combined more cheaply and effectively if the presence of mutations in prognostically important gene mutations could be identified from RNA-seq data. However, computational methodologies for robust detection of the diverse types of somatic mutations found in AML such as substitutions, indels, tandem duplications and translocations, are not currently available. Aims: To develop an stand-alone, lightweight and use-friendly software for the identification of clinically relevant mutations from AML RNA-seq data. Methods: To ensure efficient mapping of RNA-seq reads, we hash-indexed the DNA sequences of target genes using 10-mer sliding windows and implemented the ”seed and extend“ algorithm for read alignment. Point mutations and small indels were detected from reads with imperfect alignments and tandem duplications were detected from reads spanning the duplication junction. Translocations were detected from reads whose ends belonged to preselected fusion partner genes. Results: To benchmark our software we used RNA-seq data from 151 whole-exome/genome-sequenced AML samples studied by The Cancer Genome Atlas Research Network. We show that our software reliably calls clinically important mutations affecting the genes NPM1 (4-nt insertion), FLT3 (substitutions and internal tandem duplications, ITD), MLL partial tandem duplications (PTD), as well as substitutions in CEBPA, IDH1/2, TP53 and RUNX1. Furthermore, we identified gene fusions including PML-RARA, MYH11-CBFB, RUNX1-RUNX1T1, BCR-ABL1 and NUP98-NSD1. Our software is fast and memory efficient and is able to identify the above mutations in less than 20 minutes starting with RNA-seq FastQ files of 100 million 50 bp paired-end reads, using a standard modern laptop computer. In addition, the software operates through a graphical user interface making it accessible to users without programming knowledge. Summary/Conclusion: We demonstrated that clinically important somatic mutations that drives AML can be reliably detected from RNA-seq data alone using our software. As our approach can be readily combined with conventional gene expression analyses of the same RNA-seq dataset, it can be used to generate data with enhanced clinical utility that can improve prognostication and guide patient treatment. As RNA sequencing is a straightforward procedure, our approach can readily enter clinical laboratories, where it can significantly reduce experimental costs and accelerate diagnostic work-ups.

LFMD: a new likelihood-based method to detect low-frequency mutations without molecular tags

Detecting and Quantitating Low Fraction DNA Variants with Low-Depth Sequencing

DELFMUT: Duplex Sequencing-Oriented Depth Estimation Model for Stable Detection of Low-Frequency Mutations.

High efficiency error suppression for accurate detection of low-frequency variants

Application of Next-Generation Sequencing in the Detection of Low-Abundance Mutations.

Selective multiplexed enrichment for the detection and quantitation of low-fraction DNA variants via low-depth sequencing

Extended Enrichment for Ultrasensitive Detection of Low‐Frequency Mutations by Long Blocker Displacement Amplification

Ultra-specific discrimination of single-nucleotide mutations using sequestration-assisted molecular beacons.

Detecting Ultralow Frequency Mutation in Circulating Cell‐Free DNA of Early‐Stage Nonsmall Cell Lung Cancer Patients with Unique Molecular Identifiers

PS1040 DETECTION OF MRD MAY PREDICT THE OUTCOME OF PATIENTS WITH ACUTE MYELOID LEUKEMIA WITH BIALLELIC CEBPA MUTATIONS

A Novel Approach to Detect Large Indels from Targeted Sequencing Data in Clinical Cancer Setting

Detection of low-frequency mutations in clinical samples by increasing mutation abundance via the excision of wild-type sequences

Enhanced Error Suppression for Accurate Detection of Low‐Frequency Variants

Ultrafast and Highly Specific Detection of One-Base Mutated Cell-Free DNA at a Very Low Abundance

Influence of Sequencing Depth on the Fidelity and Sensitivity of 1%-5% Low-Frequency Mutation Detection and Recommendation for Standardization of Sequencing Depth

Evaluating the performance of low-frequency variant calling tools for the detection of variants from short-read deep sequencing data

BayVarC: an ultra-sensitive ctDNA variant caller using Bayesian approach

Ultrasensitive single-nucleotide polymorphism detection using target-recycled ligation, strand displacement and enzymatic amplification.

Ultrasensitive and high-efficiency screen of de novo low-frequency mutations by o2n-seq

Calibration-free NGS quantitation of mutations below 0.01% VAF

Detection of ultra-rare mutations by next-generation sequencing