ScanITD: Detecting internal tandem duplication with robust variant allele frequency estimation

DTDHM: detection of tandem duplications based on hybrid methods using next-generation sequencing data

Background

Tandem duplication (TD) is a common and important type of structural variation in the human genome. TDs have been shown to play an essential role in many diseases, including cancer. However, it is difficult to accurately detect TDs due to the uneven distribution of reads and the inherent complexity of next-generation sequencing (NGS) data.

Methods

This article proposes a method called DTDHM (detection of tandem duplications based on hybrid methods), which utilizes NGS data to detect TDs in a single sample. DTDHM builds a pipeline that integrates read depth (RD), split read (SR), and paired-end mapping (PEM) signals. To solve the problem of uneven distribution of normal and abnormal samples, DTDHM uses the K-nearest neighbor (KNN) algorithm for multi-feature classification prediction. Then, the qualified split reads and discordant reads are extracted and analyzed to achieve accurate localization of variation sites. This article compares DTDHM with three other methods on 450 simulated datasets and five real datasets.

Results

In 450 simulated data samples, DTDHM consistently maintained the highest F1-score. The average F1-score of DTDHM, SVIM, TARDIS, and TIDDIT were 80.0%, 56.2%, 43.4%, and 67.1%, respectively. The F1-score of DTDHM had a small variation range and its detection effect was the most stable and 1.2 times that of the suboptimal method. Most of the boundary biases of DTDHM fluctuated around 20 bp, and its boundary deviation detection ability was better than TARDIS and TIDDIT. In real data experiments, five real sequencing samples (NA19238, NA19239, NA19240, HG00266, and NA12891) were used to test DTDHM. The results showed that DTDHM had the highest overlap density score (ODS) and F1-score of the four methods.

Conclusions

Compared with the other three methods, DTDHM achieved excellent results in terms of sensitivity, precision, F1-score, and boundary bias. These results indicate that DTDHM can be used as a reliable tool for detecting TDs from NGS data, especially in the case of low coverage depth and tumor purity samples.

Evaluation of hypereosinophilia in a case of FLT3-mutant acute myeloid leukemia treated with gilteritinib

Acute myeloid leukemias (AMLs) frequently harbor activating mutations in Fms-like tyrosine kinase 3 (FLT3). The use of FLT3 inhibitors (FLT3i) is the standard of care for treatment of newly diagnosed and relapsed patients with AML. Differentiation responses including clinical differentiation syndrome have been previously reported with FLT3i when used as single agents in relapsed disease. We present a case of hypereosinophilia in a patient on FLT3i therapy with persistent FLT3 polymerase chain reaction (PCR) positivity in peripheral blood. We sorted mature leukocytes by lineage to determine if the eosinophils were leukemia-derived. FLT3 PCR and next-generation sequencing analysis demonstrated monocytic differentiation of the FLT3-ITD leukemic clone with reactive hypereosinophilia that was derived from a preleukemic SF3B1, FLT3 wild-type clone. Our case is the first to definitively demonstrate the emergence of clonal FLT3-ITD monocytes with FLT3i and the first to demonstrate a differentiation response following decitabine, venetoclax, and gilteritinib triplet therapy.

FLT3-ITD Measurable Residual Disease Monitoring in Acute Myeloid Leukemia Using Next-Generation Sequencing

The in-frame internal tandem duplication (ITD) of the FMS-like tyrosine kinase 3 (FLT3) gene is an important negative prognostic marker in acute myeloid leukemia (AML). FLT3-ITD monitoring is essential for patients at relapse or those receiving FLT3-targeted therapies. Fragment analysis (FA) is commonly used to detect and quantify FLT3-ITDs; however, detecting low-burden FLT3-ITDs after a treatment is challenging. We, therefore, developed a customized, next-generation sequencing (NGS)-based FLT3-ITD assay that includes a new ITD-tracing algorithm, "SEED", optimized for measurable residual disease (MRD) monitoring. NGS-SEED showed an enhanced sensitivity (0.001%) and has a superior performance over conventional fragment analysis. We further investigated the prognostic impact of MRD analyzed by NGS-SEED in AML patients who underwent allogeneic hematopoietic stem cell transplantation (HSCT). Our assay showed that the MRD assessed before and after HSCT were significantly associated with a risk of relapse and a poor overall survival, respectively, in a time-dependent analysis. Thus, this report highlighted the prognostic value of serial MRD monitoring using a sensitive method in a clinical setting of AML patients with FLT3-ITD.

Frugal alignment-free identification of FLT3-internal tandem duplications with FiLT3r

BACKGROUND

Internal tandem duplications in the FLT3 gene, termed FLT3-ITDs, are useful molecular markers in acute myeloid leukemia (AML) for patient risk stratification and follow-up. FLT3-ITDs are increasingly screened through high-throughput sequencing (HTS) raising the need for robust and efficient algorithms. We developed a new algorithm, which performs no alignment and uses little resources, to identify and quantify FLT3-ITDs in HTS data.

RESULTS

Our algorithm (FiLT3r) focuses on the k-mers from reads covering FLT3 exons 14 and 15. We show that those k-mers bring enough information to accurately detect, determine the length and quantify FLT3-ITD duplications. We compare the performances of FiLT3r to state-of-the-art alternatives and to fragment analysis, the gold standard method, on a cohort of 185 AML patients sequenced with capture-based HTS. On this dataset FiLT3r is more precise (no false positive nor false negative) than the other software evaluated. We also assess the software on public RNA-Seq data, which confirms the previous results and shows that FiLT3r requires little resources compared to other software.

CONCLUSION

FiLT3r is a free software available at https://gitlab.univ-lille.fr/filt3r/filt3r . The repository also contains a Snakefile to reproduce our experiments. We show that FiLT3r detects FLT3-ITDs better than other software while using less memory and time.

Paediatric BCOR ‐associated sarcomas with a novel long spliced internal tandem duplication of BCOR exon 15

Clear cell sarcoma of the kidney (CCSK) and primitive myxoid mesenchymal tumour of infancy (PMMTI) are paediatric sarcomas that most commonly harbour internal tandem duplications (ITDs) of exon 15 of the BCOR gene, in the range of 87–114 base pairs (bp). Some cases, instead, have BCOR‐CCNB3 or YWHAE‐NUTM2 gene fusions. About 10% of cases lack any of these genetic alterations when tested by standard methods. Two cases of CCSK and one PMMTI lacking the aforementioned mutations were analysed using Archer FusionPlex technology. Two related BCOR exon 15 RNA transcripts with ITDs of lengths 388 and 96 bp were detected in each case; only the 388 bp transcript was identified when genomic DNA was sequenced. In silico analysis of this transcript revealed acceptor and donor splice sites indicating that, at the RNA level, the 388‐bp transcript was likely spliced to form the 96‐bp transcript. The results were confirmed by Sanger sequencing using primers targeting the ITD breakpoint. This novel and unusually long ITD segment is difficult to identify by DNA sequencing using typical primer design strategies flanking entire duplicated segments because it exceeds the typical read lengths of most sequencing platforms as well as the usual fragment lengths obtained from formalin‐fixed paraffin‐embedded material. As diagnosis of CCSK and PMMTI may be challenging by morphology and immunohistochemistry alone, it is important to identify mutations in these cases. Knowledge of this novel BCOR ITD is important in relation to primer design for detection by sequencing, and using RNA versus DNA for sequencing.

Accurate Detection and Quantification of FLT3 Internal Tandem Duplications in Clinical Hybrid Capture Next-Generation Sequencing Data

FLT3 internal tandem duplications (ITDs) are found in approximately one-third of patients with acute myeloid leukemia and have important prognostic and therapeutic implications that have supported their assessment in routine clinical practice. Conventional methods for assessing FLT3-ITD status and allele burden have been primarily limited to PCR fragment size analysis because of the inherent difficulty in detecting large ITD variants by next-generation sequencing (NGS). In this study, we assess the performance of publicly available bioinformatic tools for the detection and quantification of FLT3-ITDs in clinical hybridization-capture NGS data. We found that FLT3_ITD_ext had the highest overall accuracy for detecting FLT3-ITDs and was able to accurately quantify allele burden. Although all other tools evaluated were able to detect FLT3-ITDs reasonably well, allele burden was consistently underestimated. We were able to significantly improve quantification of FLT3-ITD allelic burden independent of the detection method by utilizing soft-clipped reads and/or ITD junctional sequences. In addition, we show that identifying mutant reads by previously identified junctional sequences further improves the sensitivity of detecting FLT3-ITDs in post-treatment samples. Our results demonstrate that FLT3-ITDs can be reliably detected in clinical NGS data using available bioinformatic tools. We further describe how accurate quantification of FLT3-ITD allele burden can be added on to existing clinical NGS pipelines for routine assessment of FLT3-ITD status in patients with acute myeloid leukemia.

Comprehensive review and evaluation of computational methods for identifying FLT3-internal tandem duplication in acute myeloid leukaemia

Internal tandem duplication (ITD) of FMS-like tyrosine kinase 3 (FLT3-ITD) constitutes an independent indicator of poor prognosis in acute myeloid leukaemia (AML). AML with FLT3-ITD usually presents with poor treatment outcomes, high recurrence rate and short overall survival. Currently, polymerase chain reaction and capillary electrophoresis are widely adopted for the clinical detection of FLT3-ITD, whereas the length and mutation frequency of ITD are evaluated using fragment analysis. With the development of sequencing technology and the high incidence of FLT3-ITD mutations, a multitude of bioinformatics tools and pipelines have been developed to detect FLT3-ITD using next-generation sequencing data. However, systematic comparison and evaluation of the methods or software have not been performed. In this study, we provided a comprehensive review of the principles, functionality and limitations of the existing methods for detecting FLT3-ITD. We further compared the qualitative and quantitative detection capabilities of six representative tools using simulated and biological data. Our results will provide practical guidance for researchers and clinicians to select the appropriate FLT3-ITD detection tools and highlight the direction of future developments in this field. Availability: A Docker image with several programs pre-installed is available at https://github.com/niu-lab/docker-flt3-itd to facilitate the application of FLT3-ITD detection tools.

ScanITD: Detecting internal tandem duplication with robust variant allele frequency estimation

BACKGROUND

Internal tandem duplications (ITDs) are tandem duplications within coding exons and are important prognostic markers and drug targets for acute myeloid leukemia (AML). Next-generation sequencing has enabled the discovery of ITD at single-nucleotide resolution. ITD allele frequency is used in the risk stratification of patients with AML; higher ITD allele frequency is associated with poorer clinical outcomes. However, the ITD allele frequency data are often unavailable to treating physicians and the detection of ITDs with accurate variant allele frequency (VAF) estimation remains challenging for short-read sequencing.

RESULTS

Here we present the ScanITD approach, which performs a stepwise seed-and-realignment procedure for ITD detection with accurate VAF prediction. The evaluations on simulated and real data demonstrate that ScanITD outperforms 3 state-of-the-art ITD detectors, especially for VAF estimation. Importantly, ScanITD yields better accuracy than general-purpose structural variation callers for predicting ITD size range duplications.

CONCLUSIONS

ScanITD enables the accurate identification of ITDs with robust VAF estimation. ScanITD is written in Python and is open-source software that is freely accessible at https://github.com/ylab-hi/ScanITD.