UMI-Gen: A UMI-based read simulator for variant calling evaluation in paired-end sequencing NGS libraries

Liquid Biopsy in B and T Cell Lymphomas: From Bench to Bedside

Liquid biopsy through the analysis of circulating tumor DNA (ctDNA) is emerging as a powerful and non-invasive tool complementing tissue biopsy in lymphoma management. Whilst tissue biopsy remains the diagnostic gold standard, it fails to detect the molecular heterogeneity of the tumor's multiple compartments and poses challenges for sequential disease monitoring. In diffuse large-B-cell lymphoma (DLBCL), ctDNA facilitates non-invasive genotyping by identifying hallmark mutations (e.g., MYD88, CD79B, EZH2), enabling molecular cluster classification. Dynamic changes in ctDNA levels during DLBCL treatment correlate strongly with progression-free survival and overall survival, underscoring its value as a predictive and prognostic biomarker. In Hodgkin's lymphoma, characterized by a scarcity of malignant cells in tissue biopsies, ctDNA provides reliable molecular insights into tumor biology, response to therapy, and relapse risk. In primary central nervous system lymphoma, the detection of MYD88 L265P in ctDNA offers a highly sensitive, specific, and minimally invasive diagnostic option. Likewise, in aggressive T-cell lymphomas, ctDNA supports molecular profiling, aligns with tumor burden, and shows high concordance with tissue-based results. Ongoing and future clinical trials will be critical for validating and standardizing ctDNA applications, ultimately integrating liquid biopsy into routine clinical practice and enabling more personalized and dynamic lymphoma care.

Prospects for liquid biopsy approaches in lymphomas

Analytes within liquid biopsies have emerged as promising alternatives to traditional tissue biopsies for various malignancies, including lymphomas. This review explores the clinical applications of one such liquid biopsy analyte, circulating tumor DNA (ctDNA) in different types of lymphoma, focusing on its role in diagnosis, disease monitoring, and relapse detection. Advancements in next-generation sequencing (NGS) and machine learning have enhanced ctDNA analysis, offering a multi-omic approach to understanding tumor genetics. In lymphoma, ctDNA provides insights into tumor heterogeneity, aids in genetic profiling, and predicts treatment response. Recent studies demonstrate the prognostic value of ctDNA and its potential to improve patient outcomes by facilitating early disease detection and personalized treatment strategies Despite these advancements, challenges remain in optimizing sample collection, processing, assay sensitivity, and overall consensus workflows in order to facilitate integration into routine clinical practice.

Benchmarking UMI-aware and standard variant callers for low frequency ctDNA variant detection

BACKGROUND

Circulating tumour DNA (ctDNA) is a subset of cell free DNA (cfDNA) released by tumour cells into the bloodstream. Circulating tumour DNA has shown great potential as a biomarker to inform treatment in cancer patients. Collecting ctDNA is minimally invasive and reflects the entire genetic makeup of a patient's cancer. ctDNA variants in NGS data can be difficult to distinguish from sequencing and PCR artefacts due to low abundance, particularly in the early stages of cancer. Unique Molecular Identifiers (UMIs) are short sequences ligated to the sequencing library before amplification. These sequences are useful for filtering out low frequency artefacts. The utility of ctDNA as a cancer biomarker depends on accurate detection of cancer variants.

RESULTS

In this study, we benchmarked six variant calling tools, including two UMI-aware callers for their ability to call ctDNA variants. The standard variant callers tested included Mutect2, bcftools, LoFreq and FreeBayes. The UMI-aware variant callers benchmarked were UMI-VarCal and UMIErrorCorrect. We used both datasets with known variants spiked in at low frequencies, and datasets containing ctDNA, and generated synthetic UMI sequences for these datasets. Variant callers displayed different preferences for sensitivity and specificity. Mutect2 showed high sensitivity, while returning more privately called variants than any other caller in data without synthetic UMIs - an indicator of false positive variant discovery. In data encoded with synthetic UMIs, UMI-VarCal detected fewer putative false positive variants than all other callers in synthetic datasets. Mutect2 showed a balance between high sensitivity and specificity in data encoded with synthetic UMIs.

CONCLUSIONS

Our results indicate UMI-aware variant callers have potential to improve sensitivity and specificity in calling low frequency ctDNA variants over standard variant calling tools. There is a growing need for further development of UMI-aware variant calling tools if effective early detection methods for cancer using ctDNA samples are to be realised.

Integrated approach to generate artificial samples with low tumor fraction for somatic variant calling benchmarking

BACKGROUND

High-throughput sequencing (HTS) has become the gold standard approach for variant analysis in cancer research. However, somatic variants may occur at low fractions due to contamination from normal cells or tumor heterogeneity; this poses a significant challenge for standard HTS analysis pipelines. The problem is exacerbated in scenarios with minimal tumor DNA, such as circulating tumor DNA in plasma. Assessing sensitivity and detection of HTS approaches in such cases is paramount, but time-consuming and expensive: specialized experimental protocols and a sufficient quantity of samples are required for processing and analysis. To overcome these limitations, we propose a new computational approach specifically designed for the generation of artificial datasets suitable for this task, simulating ultra-deep targeted sequencing data with low-fraction variants and demonstrating their effectiveness in benchmarking low-fraction variant calling.

RESULTS

Our approach enables the generation of artificial raw reads that mimic real data without relying on pre-existing data by using NEAT, a fine-grained read simulator that generates artificial datasets using models learned from multiple different datasets. Then, it incorporates low-fraction variants to simulate somatic mutations in samples with minimal tumor DNA content. To prove the suitability of the created artificial datasets for low-fraction variant calling benchmarking, we used them as ground truth to evaluate the performance of widely-used variant calling algorithms: they allowed us to define tuned parameter values of major variant callers, considerably improving their detection of very low-fraction variants.

CONCLUSIONS

Our findings highlight both the pivotal role of our approach in creating adequate artificial datasets with low tumor fraction, facilitating rapid prototyping and benchmarking of algorithms for such dataset type, as well as the important need of advancing low-fraction variant calling techniques.

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

Detection of low-frequency variants with high accuracy plays an important role in biomedical research and clinical practice. However, it is challenging to do so with next-generation sequencing (NGS) approaches due to the high error rates of NGS. To accurately distinguish low-level true variants from these errors, many statistical variants calling tools for calling low-frequency variants have been proposed, but a systematic performance comparison of these tools has not yet been performed. Here, we evaluated four raw-reads-based variant callers (SiNVICT, outLyzer, Pisces, and LoFreq) and four UMI-based variant callers (DeepSNVMiner, MAGERI, smCounter2, and UMI-VarCal) considering their capability to call single nucleotide variants (SNVs) with allelic frequency as low as 0.025% in deep sequencing data. We analyzed a total of 54 simulated data with various sequencing depths and variant allele frequencies (VAFs), two reference data, and Horizon Tru-Q sample data. The results showed that the UMI-based callers, except smCounter2, outperformed the raw-reads-based callers regarding detection limit. Sequencing depth had almost no effect on the UMI-based callers but significantly influenced on the raw-reads-based callers. Regardless of the sequencing depth, MAGERI showed the fastest analysis, while smCounter2 consistently took the longest to finish the variant calling process. Overall, DeepSNVMiner and UMI-VarCal performed the best with considerably good sensitivity and precision of 88%, 100%, and 84%, 100%, respectively. In conclusion, the UMI-based callers, except smCounter2, outperformed the raw-reads-based callers in terms of sensitivity and precision. We recommend using DeepSNVMiner and UMI-VarCal for low-frequency variant detection. The results provide important information regarding future directions for reliable low-frequency variant detection and algorithm development, which is critical in genetics-based medical research and clinical applications.

Cloud-native distributed genomic pileup operations

MOTIVATION

Pileup analysis is a building block of many bioinformatics pipelines, including variant calling and genotyping. This step tends to become a bottleneck of the entire assay since the straightforward pileup implementations involve processing of all base calls from all alignments sequentially. On the other hand, a distributed version of the algorithm faces the intrinsic challenge of splitting reads-oriented file formats into self-contained partitions to avoid costly data exchange between computational nodes.

RESULTS

Here, we present a scalable, distributed and efficient implementation of a pileup algorithm that is suitable for deploying in cloud computing environments. In particular, we implemented: (i) our custom data-partitioning algorithm optimized to work with the alignment reads, (ii) a novel and unique approach to process alignment events from sequencing reads using the MD tags, (iii) the source code micro-optimizations for recurrent operations, and (iv) a modular structure of the algorithm. We have proven that our novel approach consistently and significantly outperforms other state-of-the-art distributed tools in terms of execution time (up to 6.5× faster) and memory usage (up to 2× less), resulting in a substantial cloud cost reduction. SeQuiLa is a cloud-native solution that can be easily deployed using any managed Kubernetes and Hadoop services available in public clouds, like Microsoft Azure Cloud, Google Cloud Platform, or Amazon Web Services. Together with the already implemented distributed range join and coverage calculations, our package provides end-users with a unified SQL interface for convenient analyses of population-scale genomic data in an interactive way.

AVAILABILITY AND IMPLEMENTATION

https://biodatageeks.github.io/sequila/.

cfDNA Sequencing: Technological Approaches and Bioinformatic Issues

In the era of precision medicine, it is crucial to identify molecular alterations that will guide the therapeutic management of patients. In this context, circulating tumoral DNA (ctDNA) released by the tumor in body fluids, like blood, and carrying its molecular characteristics is becoming a powerful biomarker for non-invasive detection and monitoring of cancer. Major recent technological advances, especially in terms of sequencing, have made possible its analysis, the challenge still being its reliable early detection. Different parameters, from the pre-analytical phase to the choice of sequencing technology and bioinformatic tools can influence the sensitivity of ctDNA detection.