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.

PRADA: pipeline for RNA sequencing data analysis

PRAPI: Post-Transcriptional Regulation Analysis Pipeline for Iso-Seq

MPRAnalyze: Statistical Framework for Massively Parallel Reporter Assays

MAP-RSeq: Mayo Analysis Pipeline for RNA sequencing

UTAP: User-friendly Transcriptome Analysis Pipeline

RADPRE: A Computational Program for Identification of Differential Mrna Processing Including Alternative Polyadenylation

PRAM: a novel pooling approach for discovering intergenic transcripts from large-scale RNA sequencing experiments

TRAPLINE: a standardized and automated pipeline for RNA sequencing data analysis, evaluation and annotation

Metapipeline-DNA: A Comprehensive Germline & Somatic Genomics Nextflow Pipeline

QuickRNASeq Lifts Large-Scale RNA-seq Data Analyses to the Next Level of Automation and Interactive Visualization

RADIA: RNA and DNA Integrated Analysis for Somatic Mutation Detection

An interactive web-based application for Comprehensive Analysis of RNAi-screen Data

RNASEQR—a Streamlined and Accurate RNA-seq Sequence Analysis Program

SCALA: A complete solution for multimodal analysis of single-cell Next Generation Sequencing data

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

Accelerating Single-Cell Sequencing Data Analysis with SciDAP: A User-Friendly Approach

Pipeline for RNA sequencing data analysis by combination of Nextflow and R

Web-based platform for analysis of RNA folding from high throughput chemical probing data

L-RAPiT: A Cloud-Based Computing Pipeline for the Analysis of Long-Read RNA Sequencing Data

DNAscan: a fast, computationally and memory efficient bioinformatics pipeline for the analysis of DNA next-generation-sequencing data

SeekFusion - A Clinically Validated Fusion Transcript Detection Pipeline for PCR-Based Next-Generation Sequencing of RNA