A review of somatic single nucleotide variant calling algorithms for next-generation sequencing data

Recommendations for Tumor Mutational Burden Assay Validation and Reporting: A Joint Consensus Recommendation of the Association for Molecular Pathology, College of American Pathologists, and Society for Immunotherapy of Cancer

Tumor mutational burden (TMB) has been recognized as a predictive biomarker for immunotherapy response in several tumor types. Several laboratories offer TMB testing, but there is significant variation in how TMB is calculated, reported, and interpreted among laboratories. TMB standardization efforts are underway, but no published guidance for TMB validation and reporting is currently available. Recognizing the current challenges of clinical TMB testing, the Association for Molecular Pathology convened a multidisciplinary collaborative working group with representation from the American Society of Clinical Oncology, the College of American Pathologists, and the Society for the Immunotherapy of Cancer to review the laboratory practices surrounding TMB and develop recommendations for the analytical validation and reporting of TMB testing based on survey data, literature review, and expert consensus. These recommendations encompass pre-analytical, analytical, and postanalytical factors of TMB analysis, and they emphasize the relevance of comprehensive methodological descriptions to allow comparability between assays.

A Novel Affordable and Reliable Framework for Accurate Detection and Comprehensive Analysis of Somatic Mutations in Cancer

Accurate detection and analysis of somatic variants in cancer involve multiple third-party tools with complex dependencies and configurations, leading to laborious, error-prone, and time-consuming data conversions. This approach lacks accuracy, reproducibility, and portability, limiting clinical application. Musta was developed to address these issues as an end-to-end pipeline for detecting, classifying, and interpreting cancer mutations. Musta is based on a Python command-line tool designed to manage tumor-normal samples for precise somatic mutation analysis. The core is a Snakemake-based workflow that covers all key cancer genomics steps, including variant calling, mutational signature deconvolution, variant annotation, driver gene detection, pathway analysis, and tumor heterogeneity estimation. Musta is easy to install on any system via Docker, with a Makefile handling installation, configuration, and execution, allowing for full or partial pipeline runs. Musta has been validated at the CRS4-NGS Core facility and tested on large datasets from The Cancer Genome Atlas and the Beijing Institute of Genomics. Musta has proven robust and flexible for somatic variant analysis in cancer. It is user-friendly, requiring no specialized programming skills, and enables data processing with a single command line. Its reproducibility ensures consistent results across users following the same protocol.

Unlocking the potential of Molecular Tumor Boards: from cutting-edge data interpretation to innovative clinical pathways

The emerging era of precision medicine is characterized by an increasing availability of targeted anticancer therapies and by the parallel development of techniques to obtain more refined molecular data, whose interpretation may not always be straightforward. Molecular tumor boards gather various professional figures, in order to leverage the analysis of molecular data and provide prognostic and predictive insights for clinicians. In addition to healthcare development, they could also become a tool to promote knowledge and research spreading. A growing body of evidence on the application of molecular tumor boards to clinical practice is forming and positive signals are emerging, although a certain degree of heterogeneity exists. This work analyzes molecular tumor boards' potential workflows, figures involved, data sources, sample matrices and eligible patients, as well as available evidence and learning examples. The emerging concept of multi-institutional, disease-specific molecular tumor boards is also considered by presenting two ongoing nationwide experiences.

Pan-cancer discovery of somatic mutations from RNA sequencing data

Identification of somatic mutations (SMs) is essential for characterizing cancer genomes. While DNA-seq is the prevalent method for identifying SMs, RNA-seq provides an alternative strategy to discover tumor mutations in the transcribed genome. Here, we have developed a machine learning based pipeline to discover SMs based on RNA-seq data (designated as RNA-SMs). Subsequently, we have conducted a pan-cancer analysis to systematically identify RNA-SMs from over 8,000 tumors in The Cancer Genome Atlas (TCGA). In this way, we have identified over 105,000 novel SMs that had not been reported in previous TCGA studies. These novel SMs have significant clinical implications in designing targeted therapy for improved patient outcomes. Further, we have combined the SMs identified by both RNA-seq and DNA-seq analyses to depict an updated mutational landscape across 32 cancer types. This new online SM atlas, OncoDB ( https://oncodb.org ), offers a more complete view of gene mutations that underline the development and progression of various cancers.

Accelerated somatic mutation calling for whole-genome and whole-exome sequencing data from heterogenous tumor samples

Accurate detection of somatic mutations in DNA sequencing data is a fundamental prerequisite for cancer research. Previous analytical challenges were overcome by consensus mutation calling from four to five popular callers. This, however, increases the already nontrivial computing time from individual callers. Here, we launch MuSE 2, powered by multistep parallelization and efficient memory allocation, to resolve the computing time bottleneck. MuSE 2 speeds up 50 times more than MuSE 1 and eight to 80 times more than other popular callers. Our benchmark study suggests combining MuSE 2 and the recently accelerated Strelka2 achieves high efficiency and accuracy in analyzing large cancer genomic data sets.

Detecting Somatic Mutations Without Matched Normal Samples Using Long Reads

DNA sequencing of tumours to identify somatic mutations has become a critical tool to guide the type of treatment given to cancer patients. The gold standard for mutation calling is comparing sequencing data from the tumour to a matched normal sample to avoid mis-classifying inherited SNPs as mutations. This procedure works extremely well, but in certain situations only a tumour sample is available. While approaches have been developed to find mutations without a matched normal, they have limited accuracy or require specific types of input data (e.g. ultra-deep sequencing). Here we explore the application of single molecule long read sequencing to calling somatic mutations without matched normal samples. We develop a simple theoretical framework to show how haplotype phasing is an important source of information for determining whether a variant is a somatic mutation. We then use simulations to assess the range of experimental parameters (tumour purity, sequencing depth) where this approach is effective. These ideas are developed into a prototype somatic mutation caller, smrest, and its use is demonstrated on two highly mutated cancer cell lines. Finally, we argue that this approach has potential to measure clinically important biomarkers that are based on the genome-wide distribution of mutations: tumour mutation burden and mutation signatures.

Whole-Genome Sequencing Reveals Rare Off-Target Mutations in MC1R-Edited Pigs Generated by Using CRISPR-Cas9 and Somatic Cell Nuclear Transfer

The clustered regularly interspaced short palindromic repeats (CRISPR)-Cas9 system has been widely used to create animal models for biomedical and agricultural use owing to its low cost and easy handling. However, the occurrence of erroneous cleavage (off-targeting) may raise certain concerns for the practical application of the CRISPR-Cas9 system. In this study, we created a melanocortin 1 receptor (MC1R)-edited pig model through somatic cell nuclear transfer (SCNT) by using porcine kidney cells modified by the CRISPR-Cas9 system. We then carried out whole-genome sequencing of two MC1R-edited pigs and two cloned wild-type siblings, together with the donor cells, to assess the genome-wide presence of single-nucleotide variants and small insertions and deletions (indels) and found only one candidate off-target indel in both MC1R-edited pigs. In summary, our study indicates that the minimal off-targeting effect induced by CRISPR-Cas9 may not be a major concern in gene-edited pigs created by SCNT.

Empirical Bayes single nucleotide variant-calling for next-generation sequencing data

One of the fundamental computational problems in cancer genomics is the identification of single nucleotide variants (SNVs) from DNA sequencing data. Many statistical models and software implementations for SNV calling have been developed in the literature, yet, they still disagree widely on real datasets. Based on an empirical Bayesian approach, we introduce a local false discovery rate (LFDR) estimator for germline SNV calling. Our approach learns model parameters without prior information, and simultaneously accounts for information across all sites in the genomic regions of interest. We also propose another LFDR-based algorithm that reliably prioritizes a given list of mutations called by any other variant-calling algorithm. We use a suite of gold-standard cell line data to compare our LFDR approach against a collection of widely used, state of the art programs. We find that our LFDR approach approximately matches or exceeds the performance of all of these programs, despite some very large differences among them. Furthermore, when prioritizing other algorithms' calls by our LFDR score, we find that by manipulating the type I-type II tradeoff we can select subsets of variant calls with minimal loss of sensitivity but dramatic increases in precision.

Advancements and hurdles in the development of a vaccine for triple-negative breast cancer: A comprehensive review of multi-omics and immunomics strategies

Triple-Negative Breast Cancer (TNBC) presents a significant challenge in oncology due to its aggressive behavior and limited therapeutic options. This review explores the potential of immunotherapy, particularly vaccine-based approaches, in addressing TNBC. It delves into the role of immunoinformatics in creating effective vaccines against TNBC. The review first underscores the distinct attributes of TNBC and the importance of tumor antigens in vaccine development. It then elaborates on antigen detection techniques such as exome sequencing, HLA typing, and RNA sequencing, which are instrumental in identifying TNBC-specific antigens and selecting vaccine candidates. The discussion then shifts to the in-silico vaccine development process, encompassing antigen selection, epitope prediction, and rational vaccine design. This process merges computational simulations with immunological insights. The role of Artificial Intelligence (AI) in expediting the prediction of antigens and epitopes is also emphasized. The review concludes by encapsulating how Immunoinformatics can augment the design of TNBC vaccines, integrating tumor antigens, advanced detection methods, in-silico strategies, and AI-driven insights to advance TNBC immunotherapy. This could potentially pave the way for more targeted and efficacious treatments.

Performance analysis of conventional and AI-based variant callers using short and long reads

BACKGROUND

The accurate detection of variants is essential for genomics-based studies. Currently, there are various tools designed to detect genomic variants, however, it has always been a challenge to decide which tool to use, especially when various major genome projects have chosen to use different tools. Thus far, most of the existing tools were mainly developed to work on short-read data (i.e., Illumina); however, other sequencing technologies (e.g. PacBio, and Oxford Nanopore) have recently shown that they can also be used for variant calling. In addition, with the emergence of artificial intelligence (AI)-based variant calling tools, there is a pressing need to compare these tools in terms of efficiency, accuracy, computational power, and ease of use.

RESULTS

In this study, we evaluated five of the most widely used conventional and AI-based variant calling tools (BCFTools, GATK4, Platypus, DNAscope, and DeepVariant) in terms of accuracy and computational cost using both short-read and long-read data derived from three different sequencing technologies (Illumina, PacBio HiFi, and ONT) for the same set of samples from the Genome In A Bottle project. The analysis showed that AI-based variant calling tools supersede conventional ones for calling SNVs and INDELs using both long and short reads in most aspects. In addition, we demonstrate the advantages and drawbacks of each tool while ranking them in each aspect of these comparisons.

CONCLUSION

This study provides best practices for variant calling using AI-based and conventional variant callers with different types of sequencing data.

A Bioinformatics Toolkit for Next-Generation Sequencing in Clinical Oncology

Next-generation sequencing (NGS) has taken on major importance in clinical oncology practice. With the advent of targeted therapies capable of effectively targeting specific genomic alterations in cancer patients, the development of bioinformatics processes has become crucial. Thus, bioinformatics pipelines play an essential role not only in the detection and in identification of molecular alterations obtained from NGS data but also in the analysis and interpretation of variants, making it possible to transform raw sequencing data into meaningful and clinically useful information. In this review, we aim to examine the multiple steps of a bioinformatics pipeline as used in current clinical practice, and we also provide an updated list of the necessary bioinformatics tools. This resource is intended to assist researchers and clinicians in their genetic data analyses, improving the precision and efficiency of these processes in clinical research and patient care.

How low can you go?: Methodologic considerations in clonal hematopoiesis variant calling

Clonal hematopoiesis (CH) is defined by the presence of an expanded clonal hematopoietic cell population due to an acquired mutation conferring a selective growth advantage and is known to predispose to hematologic malignancy. In this review, we discuss sequencing methods for CH detection in bulk sequencing data and corresponding bioinformatic approaches for variant calling, filtering, and curation. We detail practical recommendations for CH calling. Finally, we discuss how improvements in CH sequencing and bioinformatic approaches will enable the characterization of CH trajectories, its impact on human health, and therapeutic approaches to mitigate its adverse effects.

Validation and benchmarking of targeted panel sequencing for cancer genomic profiling

OBJECTIVES

To validate a large next-generation sequencing (NGS) panel for comprehensive genomic profiling and improve patient access to more effective precision oncology treatment strategies.

METHODS

OncoPanScan was designed by targeting 825 cancer-related genes to detect a broad range of genomic alterations. A practical validation strategy was used to evaluate the assay's analytical performance, involving 97 tumor specimens with 25 paired blood specimens, 10 engineered cell lines, and 121 artificial reference DNA samples.

RESULTS

Overall, 1107 libraries were prepared and the sequencing failure rate was 0.18%. Across alteration classes, sensitivity ranged from 0.938 to more than 0.999, specificity ranged from 0.889 to more than 0.999, positive predictive value ranged from 0.867 to more than 0.999, repeatability ranged from 0.908 to more than 0.999, and reproducibility ranged from 0.832 to more than 0.999. The limit of detection for variants was established based on variant frequency, while for tumor mutation burden and microsatellite instability, it was based on tumor content, resulting in a minimum requirement of 20% tumor content. Benchmarking variant calls against validated NGS assays revealed that variations in the dry-bench processes were the primary cause of discordances.

CONCLUSIONS

This study presents a detailed validation framework and empirical recommendations for large panel validation and elucidates the sources of discordant alteration calls by comparing with "gold standard measures."

Genomic variant benchmark: if you cannot measure it, you cannot improve it

Genomic benchmark datasets are essential to driving the field of genomics and bioinformatics. They provide a snapshot of the performances of sequencing technologies and analytical methods and highlight future challenges. However, they depend on sequencing technology, reference genome, and available benchmarking methods. Thus, creating a genomic benchmark dataset is laborious and highly challenging, often involving multiple sequencing technologies, different variant calling tools, and laborious manual curation. In this review, we discuss the available benchmark datasets and their utility. Additionally, we focus on the most recent benchmark of genes with medical relevance and challenging genomic complexity.

UniAligner: a parameter-free framework for fast sequence alignment

Even though the recent advances in 'complete genomics' revealed the previously inaccessible genomic regions, analysis of variations in centromeres and other extra-long tandem repeats (ETRs) faces an algorithmic challenge since there are currently no tools for accurate sequence comparison of ETRs. Counterintuitively, the classical alignment approaches, such as the Smith-Waterman algorithm, fail to construct biologically adequate alignments of ETRs. We present UniAligner-the parameter-free sequence alignment algorithm with sequence-dependent alignment scoring that automatically changes for any pair of compared sequences. UniAligner prioritizes matches of rare substrings that are more likely to be relevant to the evolutionary relationship between two sequences. We apply UniAligner to estimate the mutation rates in human centromeres, and quantify the extremely high rate of large duplications and deletions in centromeres. This high rate suggests that centromeres may represent some of the most rapidly evolving regions of the human genome with respect to their structural organization.

Diagnostic Potential of Endometrial Cancer DNA from Pipelle, Pap-Brush, and Swab Sampling

Endometrial cancer (EC) is a major gynecological malignancy with rising morbidity and mortality worldwide. The aim of this study was to explore a safe and readily available sample and a sensitive and effective detection method and its biomarkers for early diagnosis of EC, which is critical for patient prognosis. This study designed a panel targeting variants for EC-related genes, assessed its technical performance by comparing it with whole-exon sequencing, and explored the diagnostic potential of endometrial biopsies using the Pipelle aspirator, cervical samples using the Pap brush, and vaginal specimens using the swab from 38 EC patients and 208 women with risk factors for EC by applying targeted panel sequencing (TPS). TPS produced high-quality data (Q30 > 85% and mapping ratios > 99.35%) and was found to have strong consistency with whole-exome sequencing (WES) in detecting pathogenic mutations (92.11%), calculating homologous recombination deficiency (HRD) scores (r = 0.65), and assessing the microsatellite instability (MSI) status of EC (100%). The sensitivity of TPS in detection of EC is slightly better than that of WES (86.84% vs. 84.21%). Of the three types of samples detected using TPS, endometrial biopsy using the Pipelle aspirator had the highest sensitivity in detection of pathogenic mutations (81.87%) and the best consistency with surgical tumor specimens in MSI (85.16%). About 84% of EC patients contained pathogenic mutations in PIK3CA, PTEN, TP53, ARID1A, CTNNB1, KRAS, and MTOR, suggesting that this small gene set can achieve an excellent pathogenic mutation detection rate in Chinese EC patients. The custom panel combined with ultra-deep sequencing serves as a sensitive method for detecting genetic lesions from endometrial biopsy using the Pipelle aspirator.

Accelerated somatic mutation calling for whole-genome and whole-exome sequencing data from heterogenous tumor samples

Accurate detection of somatic mutations in DNA sequencing data is a fundamental prerequisite for cancer research. Previous analytical challenge was overcome by consensus mutation calling from four to five popular callers. This, however, increases the already nontrivial computing time from individual callers. Here, we launch MuSE2.0, powered by multi-step parallelization and efficient memory allocation, to resolve the computing time bottleneck. MuSE2.0 speeds up 50 times than MuSE1.0 and 8-80 times than other popular callers. Our benchmark study suggests combining MuSE2.0 and the recently expedited Strelka2 can achieve high efficiency and accuracy in analyzing large cancer genomic datasets.

Simple combination of multiple somatic variant callers to increase accuracy

Publications comparing variant caller algorithms present discordant results with contradictory rankings. Caller performances are inconsistent and wide ranging, and dependent upon input data, application, parameter settings, and evaluation metric. With no single variant caller emerging as a superior standard, combinations or ensembles of variant callers have appeared in the literature. In this study, a whole genome somatic reference standard was used to derive principles to guide strategies for combining variant calls. Then, manually annotated variants called from the whole exome sequencing of a tumor were used to corroborate these general principles. Finally, we examined the ability of these principles to reduce noise in targeted sequencing.

What makes TMB an ambivalent biomarker for immunotherapy? A subtle mismatch between the sample-based design of variant callers and real clinical cohort

Tumor mutation burden (TMB) is a widely recognized biomarker for predicting the efficacy of immunotherapy. However, its use still remains highly controversial. In this study, we examine the underlying causes of this controversy based on clinical needs. By tracing the source of the TMB errors and analyzing the design philosophy behind variant callers, we identify the conflict between the incompleteness of biostatistics rules and the variety of clinical samples as the critical issue that renders TMB an ambivalent biomarker. A series of experiments were conducted to illustrate the challenges of mutation detection in clinical practice. Additionally, we also discuss potential strategies for overcoming these conflict issues to enable the application of TMB in guiding decision-making in real clinical settings.

cfSNV: a software tool for the sensitive detection of somatic mutations from cell-free DNA

Cell-free DNA (cfDNA) in blood, viewed as a surrogate for tumor biopsy, has many clinical applications, including diagnosing cancer, guiding cancer treatment and monitoring treatment response. All these applications depend on an indispensable, yet underdeveloped task: detecting somatic mutations from cfDNA. The task is challenging because of the low tumor fraction in cfDNA. Recently, we developed the computational method cfSNV, the first method that comprehensively considers the properties of cfDNA for the sensitive detection of mutations from cfDNA. cfSNV vastly outperformed the conventional methods that were developed primarily for calling mutations from solid tumor tissues. cfSNV can accurately detect mutations in cfDNA even with medium-coverage (e.g., ≥200×) sequencing, which makes whole-exome sequencing (WES) of cfDNA a viable option for various clinical utilities. Here, we present a user-friendly cfSNV package that exhibits fast computation and convenient user options. We also built a Docker image of it, which is designed to enable researchers and clinicians with a limited computational background to easily carry out analyses on both high-performance computing platforms and local computers. Mutation calling from a standard preprocessed WES dataset (~250× and ~70 million base pair target size) can be carried out in 3 h on a server with eight virtual CPUs and 32 GB of random access memory.

Validation of genetic variants from NGS data using deep convolutional neural networks

Accurate somatic variant calling from next-generation sequencing data is one most important tasks in personalised cancer therapy. The sophistication of the available technologies is ever-increasing, yet, manual candidate refinement is still a necessary step in state-of-the-art processing pipelines. This limits reproducibility and introduces a bottleneck with respect to scalability. We demonstrate that the validation of genetic variants can be improved using a machine learning approach resting on a Convolutional Neural Network, trained using existing human annotation. In contrast to existing approaches, we introduce a way in which contextual data from sequencing tracks can be included into the automated assessment. A rigorous evaluation shows that the resulting model is robust and performs on par with trained researchers following published standard operating procedure.

Computational prediction of MHC anchor locations guides neoantigen identification and prioritization

Neoantigens are tumor-specific peptide sequences resulting from sources such as somatic DNA mutations. Upon loading onto major histocompatibility complex (MHC) molecules, they can trigger recognition by T cells. Accurate neoantigen identification is thus critical for both designing cancer vaccines and predicting response to immunotherapies. Neoantigen identification and prioritization relies on correctly predicting whether the presenting peptide sequence can successfully induce an immune response. Because most somatic mutations are single-nucleotide variants, changes between wild-type and mutated peptides are typically subtle and require cautious interpretation. A potentially underappreciated variable in neoantigen prediction pipelines is the mutation position within the peptide relative to its anchor positions for the patient's specific MHC molecules. Whereas a subset of peptide positions are presented to the T cell receptor for recognition, others are responsible for anchoring to the MHC, making these positional considerations critical for predicting T cell responses. We computationally predicted anchor positions for different peptide lengths for 328 common HLA alleles and identified unique anchoring patterns among them. Analysis of 923 tumor samples shows that 6 to 38% of neoantigen candidates are potentially misclassified and can be rescued using allele-specific knowledge of anchor positions. A subset of anchor results were orthogonally validated using protein crystallography structures. Representative anchor trends were experimentally validated using peptide-MHC stability assays and competition binding assays. By incorporating our anchor prediction results into neoantigen prediction pipelines, we hope to formalize, streamline, and improve the identification process for relevant clinical studies.

In Silico Pipeline to Identify Tumor-Specific Antigens for Cancer Immunotherapy Using Exome Sequencing Data

Tumor-specific antigens or neoantigens are peptides that are expressed only in cancer cells and not in healthy cells. Some of these molecules can induce an immune response, and therefore, their use in immunotherapeutic strategies based on cancer vaccines has been extensively explored. Studies based on these approaches have been triggered by the current high-throughput DNA sequencing technologies. However, there is no universal nor straightforward bioinformatic protocol to discover neoantigens using DNA sequencing data. Thus, we propose a bioinformatic protocol to detect tumor-specific antigens associated with single nucleotide variants (SNVs) or "mutations" in tumoral tissues. For this purpose, we used publicly available data to build our model, including exome sequencing data from colorectal cancer and healthy cells obtained from a single case, as well as frequent human leukocyte antigen (HLA) class I alleles in a specific population. HLA data from Costa Rican Central Valley population was selected as an example. The strategy included three main steps: (1) pre-processing of sequencing data; (2) variant calling analysis to detect tumor-specific SNVs in comparison with healthy tissue; and (3) prediction and characterization of peptides (protein fragments, the tumor-specific antigens) derived from the variants, in the context of their affinity with frequent alleles of the selected population. In our model data, we found 28 non-silent SNVs, present in 17 genes in chromosome one. The protocol yielded 23 strong binders peptides derived from the SNVs for frequent HLA class I alleles for the Costa Rican population. Although the analyses were performed as an example to implement the pipeline, to our knowledge, this is the first study of an in silico cancer vaccine using DNA sequencing data in the context of the HLA alleles. It is concluded that the standardized protocol was not only able to identify neoantigens in a specific but also provides a complete pipeline for the eventual design of cancer vaccines using the best bioinformatic practices.

Supplementary Information

The online version contains supplementary material available at 10.1007/s43657-022-00084-9.

Evaluation of ctDNA in the Prediction of Response to Neoadjuvant Therapy and Prognosis in Locally Advanced Rectal Cancer Patients: A Prospective Study

“Watch and wait” is becoming a common treatment option for patients with locally advanced rectal cancer (LARC) submitted to neoadjuvant treatment. However, currently, no clinical modality has an acceptable accuracy for predicting pathological complete response (pCR). The aim of this study was to assess the clinical utility of circulating tumor DNA (ctDNA) in predicting the response and prognosis in these patients. We prospectively enrolled a cohort of three Iberian centers between January 2020 and December 2021 and performed an analysis on the association of ctDNA with the main response outcomes and disease-free survival (DFS). The rate of pCR in the total sample was 15.3%. A total of 24 plasma samples from 18 patients were analyzed by next-generation sequencing. At baseline, mutations were detected in 38.9%, with the most common being TP53 and KRAS. Combination of either positive magnetic resonance imaging (MRI) extramural venous invasion (mrEMVI) and ctDNA increased the risk of poor response (p = 0.021). Also, patients with two mutations vs. those with fewer than two mutations had a worse DFS (p = 0.005). Although these results should be read carefully due to sample size, this study suggests that baseline ctDNA combined with mrEMVI could potentially help to predict the response and baseline ctDNA number of mutations might allow the discrimination of groups with different DFS. Further studies are needed to clarify the role of ctDNA as an independent tool in the selection and management of LARC patients.

A Clinician's Guide to Bioinformatics for Next-Generation Sequencing

Next-generation sequencing (NGS) technologies are high-throughput methods for DNA sequencing and have become a widely adopted tool in cancer research. The sheer amount and variety of data generated by NGS assays require sophisticated computational methods and bioinformatics expertise. In this review, we provide background details of NGS technology and basic bioinformatics concepts for the clinician investigator interested in cancer research applications, with a focus on DNA-based approaches. We introduce the general principles of presequencing library preparation, postsequencing alignment, and variant calling. We also highlight the common variant annotations and NGS applications for other molecular data types. Finally, we briefly discuss the revealed utility of NGS methods in NSCLC research and study design considerations for research studies that aim to leverage NGS technologies for clinical care.

Neoepitopes prediction strategies: an integration of cancer genomics and immunoinformatics approaches

A major near-term medical impact of the genomic technology revolution will be the elucidation of mechanisms of cancer pathogenesis, leading to improvements in the diagnosis of cancer and the selection of cancer treatment. Next-generation sequencing technologies have accelerated the characterization of a tumor, leading to the comprehensive discovery of all the major alterations in a given cancer genome, followed by the translation of this information using computational and immunoinformatics approaches to cancer diagnostics and therapeutic efforts. In the current article, we review various components of cancer immunoinformatics applied to a series of fields of cancer research, including computational tools for cancer mutation detection, cancer mutation and immunological databases, and computational vaccinology.

DeepSom: a CNN-based approach to somatic variant calling in WGS samples without a matched normal

MOTIVATION

Somatic mutations are usually called by analyzing the DNA sequence of a tumor sample in conjunction with a matched normal. However, a matched normal is not always available, for instance, in retrospective analysis or diagnostic settings. For such cases, tumor-only somatic variant calling tools need to be designed. Previously proposed approaches demonstrate inferior performance on whole-genome sequencing (WGS) samples.

RESULTS

We present the convolutional neural network-based approach called DeepSom for detecting somatic single nucleotide polymorphism and short insertion and deletion variants in tumor WGS samples without a matched normal. We validate DeepSom by reporting its performance on five different cancer datasets. We also demonstrate that on WGS samples DeepSom outperforms previously proposed methods for tumor-only somatic variant calling.

AVAILABILITY AND IMPLEMENTATION

DeepSom is available as a GitHub repository at https://github.com/heiniglab/DeepSom.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Reconstructing Phylogenetic Relationship in Bladder Cancer: A Methodological Overview

Bladder cancer (BC) expresses itself as a highly heterogeneous disease both at the histological and molecular level, often occurring as synchronous or metachronous multifocal disease with high risk of recurrence and potential to metastasize. Multiple sequencing studies focusing on both non-muscle-invasive bladder cancer (NMIBC) and muscle-invasive bladder cancer (MIBC) gave insights into the extent of both inter- and intrapatient heterogeneity, but many questions on clonal evolution in BC remain unanswered. In this review article, we provide an overview over the technical and theoretical concepts linked to reconstructing evolutionary trajectories in BC and propose a set of tools and established software for phylogenetic analysis.

Optimizing Insertion and Deletion Detection Using Next-Generation Sequencing in the Clinical Laboratory

Detection of insertions and deletions (InDels) by short-read next-generation sequencing (NGS) technology can be challenging because of frequent misaligned reads. A systematic analysis of short InDels (1 to 30 bases) and fms-related receptor tyrosine kinase 3 (FLT3) internal tandem duplications (ITDs; 6 to 183 bases) from 46 clinical cases of solid or hematologic malignancy processed with a clinical NGS assay identified misaligned reads in every case, ranging from 3% to 100% of reads with the InDel showing mismapped bases. Mismaps also increased with InDel size. As a consequence, the clinical NGS bioinformatics pipeline undercalled the variant allele frequency by 1% to 84%, incorrectly called simultaneous single-base substitutions along with InDels, or did not report an FLT3 ITD that had been detected by capillary electrophoresis. To improve the ability of the pipeline to better detect and quantify InDels, we utilized a software program called Assembly-Based ReAligner (ABRA2) to more accurately remap reads. ABRA2 was able to correct 41% to 100% of the reads with mismapped bases and led to absolute increases in the variant allele frequency from 1% to 61% along with correction of all of the single-base substitutions except for two cases. ABRA2 could also detect multiple FLT3 ITD clones except for one 183-base ITD. Our analysis has found that ABRA2 performs well on short InDels as well as FLT3 ITDs that are <100 bases.

From Samples to Germline and Somatic Sequence Variation: A Focus on Next-Generation Sequencing in Melanoma Research

Next-generation sequencing (NGS) applications have flourished in the last decade, permitting the identification of cancer driver genes and profoundly expanding the possibilities of genomic studies of cancer, including melanoma. Here we aimed to present a technical review across many of the methodological approaches brought by the use of NGS applications with a focus on assessing germline and somatic sequence variation. We provide cautionary notes and discuss key technical details involved in library preparation, the most common problems with the samples, and guidance to circumvent them. We also provide an overview of the sequence-based methods for cancer genomics, exposing the pros and cons of targeted sequencing vs. exome or whole-genome sequencing (WGS), the fundamentals of the most common commercial platforms, and a comparison of throughputs and key applications. Details of the steps and the main software involved in the bioinformatics processing of the sequencing results, from preprocessing to variant prioritization and filtering, are also provided in the context of the full spectrum of genetic variation (SNVs, indels, CNVs, structural variation, and gene fusions). Finally, we put the emphasis on selected bioinformatic pipelines behind (a) short-read WGS identification of small germline and somatic variants, (b) detection of gene fusions from transcriptomes, and (c) de novo assembly of genomes from long-read WGS data. Overall, we provide comprehensive guidance across the main methodological procedures involved in obtaining sequencing results for the most common short- and long-read NGS platforms, highlighting key applications in melanoma research.

Detection and evolutionary dynamics of somatic FAS variants in autoimmune lymphoproliferative syndrome: Diagnostic implications

Autoimmune lymphoproliferative syndrome (ALPS) is a rare primary immune disorder characterized by impaired apoptotic homeostasis. The clinical characteristics include lymphoproliferation, autoimmunity (mainly cytopenia), and an increased risk of lymphoma. A distinctive biological feature is accumulation (&gt;2.5%) of an abnormal cell subset composed of TCRαβ+ CD4-CD8- T cells (DNTs). The most common genetic causes of ALPS are monoallelic pathogenic variants in the FAS gene followed by somatic FAS variants, mainly restricted to DNTs. Identification of somatic FAS variants has been typically addressed by Sanger sequencing in isolated DNTs. However, this approach can be costly and technically challenging, and may not be successful in patients with normal DNT counts receiving immunosuppressive treatment. In this study, we identified a novel somatic mutation in FAS (c.718_719insGTCG) by Sanger sequencing on purified CD3+ cells. We then followed the evolutionary dynamics of the variant along time with an NGS-based approach involving deep amplicon sequencing (DAS) at high coverage (20,000-30,000x). Over five years of clinical follow-up, we obtained six blood samples for molecular study from the pre-treatment (DNTs&gt;7%) and treatment (DNTs&lt;2%) periods. DAS enabled detection of the somatic variant in all samples, even the one obtained after five years of immunosuppressive treatment (DNTs: 0.89%). The variant allele frequency (VAF) range was 4%-5% in pre-treatment samples and &lt;1.5% in treatment samples, and there was a strong positive correlation between DNT counts and VAF (Pearson’s R: 0.98, p=0.0003). We then explored whether the same approach could be used in a discovery setting. In the last follow-up sample (DNT: 0.89%) we performed somatic variant calling on the FAS exon 9 DAS data from whole blood and purified CD3+ cells using VarScan 2. The c.718_719insGTCG variant was identified in both samples and showed the highest VAF (0.67% blood, 1.58% CD3+ cells) among &gt;400 variants called. In summary, our study illustrates the evolutionary dynamics of a somatic FAS mutation before and during immunosuppressive treatment. The results show that pathogenic somatic FAS variants can be identified with the use of DAS in whole blood of ALPS patients regardless of their DNT counts.

Mycobacterium tuberculosis genetic features associated with pulmonary tuberculosis severity

OBJECTIVES

Mycobacterium tuberculosis (Mtb) infections result in a wide spectrum of clinical presentations but without proven Mtb genetic determinants. Herein, we hypothesized that the genetic features of Mtb clinical isolates, such as specific polymorphisms or microdiversity, may be linked to tuberculosis (TB) severity.

METHODS

A total of 234 patients with pulmonary TB (including 193 drug-susceptible and 14 monoresistant cases diagnosed between 2017 and 2020 and 27 multidrug-resistant cases diagnosed between 2010 and 2020) were stratified according to TB disease severity, and Mtb genetic features were explored using whole genome sequencing, including heterologous single-nucleotide polymorphism (SNP), calling to explore microdiversity. Finally, we performed a structural equation modeling analysis to relate TB severity to Mtb genetic features.

RESULTS

The clinical isolates from patients with mild TB carried mutations in genes associated with host-pathogen interaction, whereas those from patients with moderate/severe TB carried mutations associated with regulatory mechanisms. Genome-wide association study identified an SNP in the promoter of the gene coding for the virulence regulator espR, statistically associated with moderate/severe disease. Structural equation modeling and model comparisons indicated that TB severity was associated with the detection of Mtb microdiversity within clinical isolates and to the espR SNP.

CONCLUSION

Taken together, these results provide a new insight to better understand TB pathophysiology and could provide a new prognosis tool for pulmonary TB severity.

grenepipe: a flexible, scalable and reproducible pipeline to automate variant calling from sequence reads

SUMMARY

We developed grenepipe, an all-in-one Snakemake workflow to streamline the data processing from raw high-throughput sequencing data of individuals or populations to genotype variant calls. Our pipeline offers a range of popular software tools within a single configuration file, automatically installs software dependencies, is highly optimized for scalability in cluster environments and runs with a single command.

AVAILABILITY AND IMPLEMENTATION

grenepipe is published under the GPLv3 and freely available at github.com/moiexpositoalonsolab/grenepipe.

Neoantigens in precision cancer immunotherapy: from identification to clinical applications

Immunotherapies targeting cancer neoantigens are safe, effective, and precise. Neoantigens can be identified mainly by genomic techniques such as next-generation sequencing and high-throughput single-cell sequencing; proteomic techniques such as mass spectrometry; and bioinformatics tools based on high-throughput sequencing data, mass spectrometry data, and biological databases. Neoantigen-related therapies are widely used in clinical practice and include neoantigen vaccines, neoantigen-specific CD8+ and CD4+ T cells, and neoantigen-pulsed dendritic cells. In addition, neoantigens can be used as biomarkers to assess immunotherapy response, resistance, and prognosis. Therapies based on neoantigens are an important and promising branch of cancer immunotherapy. Unremitting efforts are needed to unravel the comprehensive role of neoantigens in anti-tumor immunity and to extend their clinical application. This review aimed to summarize the progress in neoantigen research and to discuss its opportunities and challenges in precision cancer immunotherapy.

RNA-SSNV: A Reliable Somatic Single Nucleotide Variant Identification Framework for Bulk RNA-Seq Data

The usage of expressed somatic mutations may have a unique advantage in identifying active cancer driver mutations. However, accurately calling mutations from RNA-seq data is difficult due to confounding factors such as RNA-editing, reverse transcription, and gap alignment. In the present study, we proposed a framework (named RNA-SSNV, https://github.com/pmglab/RNA-SSNV) to call somatic single nucleotide variants (SSNV) from tumor bulk RNA-seq data. Based on a comprehensive multi-filtering strategy and a machine-learning classification model trained with comprehensively curated features, RNA-SSNV achieved the best precision–recall rate (0.880–0.884) in a testing dataset and robustly retained 0.94 AUC for the precision–recall curve in three validation adult-based TCGA (The Cancer Genome Atlas) datasets. We further showed that the somatic mutations called by RNA-SSNV tended to have a higher functional impact and therapeutic power in known driver genes. Furthermore, VAF (variant allele fraction) analysis revealed that subclonal harboring expressed mutations had evolutional selection advantage and RNA had higher detection power to rescue DNA-omitted mutations. In sum, RNA-SSNV will be a useful approach to accurately call expressed somatic mutations for a more insightful analysis of cancer drive genes and carcinogenic mechanisms.

Recent Progress on Therapeutic Vaccines for Breast Cancer

Breast cancer remains the most frequently diagnosed malignancy worldwide. Advanced breast cancer is still an incurable disease mainly because of its heterogeneity and limited immunogenicity. The great success of cancer immunotherapy is paving the way for a new era in cancer treatment, and therapeutic cancer vaccination is an area of interest. Vaccine targets include tumor-associated antigens and tumor-specific antigens. Immune responses differ in different vaccine delivery platforms. Next-generation sequencing technologies and computational analysis have recently made personalized vaccination possible. However, only a few cases benefiting from neoantigen-based treatment have been reported in breast cancer, and more attention has been given to overexpressed antigen-based treatment, especially human epidermal growth factor 2-derived peptide vaccines. Here, we discuss recent advancements in therapeutic vaccines for breast cancer and highlight near-term opportunities for moving forward.

Artificial intelligence and machine learning in precision and genomic medicine

The advancement of precision medicine in medical care has led behind the conventional symptom-driven treatment process by allowing early risk prediction of disease through improved diagnostics and customization of more effective treatments. It is necessary to scrutinize overall patient data alongside broad factors to observe and differentiate between ill and relatively healthy people to take the most appropriate path toward precision medicine, resulting in an improved vision of biological indicators that can signal health changes. Precision and genomic medicine combined with artificial intelligence have the potential to improve patient healthcare. Patients with less common therapeutic responses or unique healthcare demands are using genomic medicine technologies. AI provides insights through advanced computation and inference, enabling the system to reason and learn while enhancing physician decision making. Many cell characteristics, including gene up-regulation, proteins binding to nucleic acids, and splicing, can be measured at high throughput and used as training objectives for predictive models. Researchers can create a new era of effective genomic medicine with the improved availability of a broad range of datasets and modern computer techniques such as machine learning. This review article has elucidated the contributions of ML algorithms in precision and genome medicine.

Artificial intelligence and machine learning in precision and genomic medicine

The advancement of precision medicine in medical care has led behind the conventional symptom-driven treatment process by allowing early risk prediction of disease through improved diagnostics and customization of more effective treatments. It is necessary to scrutinize overall patient data alongside broad factors to observe and differentiate between ill and relatively healthy people to take the most appropriate path toward precision medicine, resulting in an improved vision of biological indicators that can signal health changes. Precision and genomic medicine combined with artificial intelligence have the potential to improve patient healthcare. Patients with less common therapeutic responses or unique healthcare demands are using genomic medicine technologies. AI provides insights through advanced computation and inference, enabling the system to reason and learn while enhancing physician decision making. Many cell characteristics, including gene up-regulation, proteins binding to nucleic acids, and splicing, can be measured at high throughput and used as training objectives for predictive models. Researchers can create a new era of effective genomic medicine with the improved availability of a broad range of datasets and modern computer techniques such as machine learning. This review article has elucidated the contributions of ML algorithms in precision and genome medicine.

Combinatorial and Machine Learning Approaches for Improved Somatic Variant Calling From Formalin-Fixed Paraffin-Embedded Genome Sequence Data

Formalin fixation of paraffin-embedded tissue samples is a well-established method for preserving tissue and is routinely used in clinical settings. Although formalin-fixed, paraffin-embedded (FFPE) tissues are deemed crucial for research and clinical applications, the fixation process results in molecular damage to nucleic acids, thus confounding their use in genome sequence analysis. Methods to improve genomic data quality from FFPE tissues have emerged, but there remains significant room for improvement. Here, we use whole-genome sequencing (WGS) data from matched Fresh Frozen (FF) and FFPE tissue samples to optimize a sensitive and precise FFPE single nucleotide variant (SNV) calling approach. We present methods to reduce the prevalence of false-positive SNVs by applying combinatorial techniques to five publicly available variant callers. We also introduce FFPolish, a novel variant classification method that efficiently classifies FFPE-specific false-positive variants. Our combinatorial and statistical techniques improve precision and F1 scores compared to the results of publicly available tools when tested individually.

Detection of oncogenic and clinically actionable mutations in cancer genomes critically depends on variant calling tools

Motivation

The analysis of cancer genomes provides fundamental information about its etiology, the processes driving cell transformation or potential treatments. While researchers and clinicians are often only interested in the identification of oncogenic mutations, actionable variants or mutational signatures, the first crucial step in the analysis of any tumor genome is the identification of somatic variants in cancer cells (i.e. those that have been acquired during their evolution). For that purpose, a wide range of computational tools have been developed in recent years to detect somatic mutations in sequencing data from tumor samples. While there have been some efforts to benchmark somatic variant calling tools and strategies, the extent to which variant calling decisions impact the results of downstream analyses of tumor genomes remains unknown.

Results

Here, we quantify the impact of variant calling decisions by comparing the results obtained in three important analyses of cancer genomics data (identification of cancer driver genes, quantification of mutational signatures and detection of clinically actionable variants) when changing the somatic variant caller (MuSE, MuTect2, SomaticSniper and VarScan2) or the strategy to combine them (Consensus of two, Consensus of three and Union) across all 33 cancer types from The Cancer Genome Atlas. Our results show that variant calling decisions have a significant impact on these analyses, creating important differences that could even impact treatment decisions for some patients. Moreover, the Consensus of three calling strategy to combine the output of multiple variant calling tools, a very widely used strategy by the research community, can lead to the loss of some cancer driver genes and actionable mutations. Overall, our results highlight the limitations of widespread practices within the cancer genomics community and point to important differences in critical analyses of tumor sequencing data depending on variant calling, affecting even the identification of clinically actionable variants.

Availability and implementation

Code is available at https://github.com/carlosgarciaprieto/VariantCallingClinicalBenchmark.

Supplementary information

Supplementary data are available at Bioinformatics online.

Computational analysis of cancer genome sequencing data

Distilling biologically meaningful information from cancer genome sequencing data requires comprehensive identification of somatic alterations using rigorous computational methods. As the amount and complexity of sequencing data have increased, so has the number of tools for analysing them. Here, we describe the main steps involved in the bioinformatic analysis of cancer genomes, review key algorithmic developments and highlight popular tools and emerging technologies. These tools include those that identify point mutations, copy number alterations, structural variations and mutational signatures in cancer genomes. We also discuss issues in experimental design, the strengths and limitations of sequencing modalities and methodological challenges for the future.

Somatic Mutations Detected in Parkinson Disease Could Affect Genes With a Role in Synaptic and Neuronal Processes

The role of somatic mutations in complex diseases, including neurodevelopmental and neurodegenerative disorders, is becoming increasingly clear. However, to date, no study has shown their relation to Parkinson disease’s phenotype. To explore the relevance of embryonic somatic mutations in sporadic Parkinson disease, we performed whole-exome sequencing in blood and four brain regions of ten patients. We identified 59 candidate somatic single nucleotide variants (sSNVs) through sensitive calling and a careful filtering strategy (COSMOS). We validated 27 of them with amplicon-based ultra-deep sequencing, with a 70% validation rate for the highest-confidence variants. The identified sSNVs are in genes with synaptic functions that are co-expressed with genes previously associated with Parkinson disease. Most of the sSNVs were only called in blood but were also found in the brain tissues with ultra-deep amplicon sequencing, demonstrating the strength of multi-tissue sampling designs.

Network Approaches for Charting the Transcriptomic and Epigenetic Landscape of the Developmental Origins of Health and Disease

The early developmental phase is of critical importance for human health and disease later in life. To decipher the molecular mechanisms at play, current biomedical research is increasingly relying on large quantities of diverse omics data. The integration and interpretation of the different datasets pose a critical challenge towards the holistic understanding of the complex biological processes that are involved in early development. In this review, we outline the major transcriptomic and epigenetic processes and the respective datasets that are most relevant for studying the periconceptional period. We cover both basic data processing and analysis steps, as well as more advanced data integration methods. A particular focus is given to network-based methods. Finally, we review the medical applications of such integrative analyses.

Consistency and reproducibility of large panel next-generation sequencing: Multi-laboratory assessment of somatic mutation detection on reference materials with mismatch repair and proofreading deficiency

INTRODUCTION

Clinical precision oncology increasingly relies on accurate genome-wide profiling using large panel next generation sequencing; however, difficulties in accurate and consistent detection of somatic mutation from individual platforms and pipelines remain an open question.

OBJECTIVES

To obtain paired tumor-normal reference materials that can be effectively constructed and interchangeable with clinical samples, and evaluate the performance of 56 panels under routine testing conditions based on the reference samples.

METHODS

Genes involved in mismatch repair and DNA proofreading were knocked down using the CRISPR-Cas9 technology to accumulate somatic mutations in a defined GM12878 cell line. They were used as reference materials to comprehensively evaluate the reproducibility and accuracy of detection results of oncopanels and explore the potential influencing factors.

RESULTS

In total, 14 paired tumor-normal reference DNA samples from engineered cell lines were prepared, and a reference dataset comprising 168 somatic mutations in a high-confidence region of 1.8 Mb were generated. For mutations with an allele frequency (AF) of more than 5% in reference samples, 56 panels collectively reported 1306 errors, including 729 false negatives (FNs), 179 false positives (FPs) and 398 reproducibility errors. The performance metric varied among panels with precision and recall ranging from 0.773 to 1 and 0.683 to 1, respectively. Incorrect and inadequate filtering accounted for a large proportion of false discovery (including FNs and FPs), while low-quality detection, cross-contamination and other sequencing errors during the wet bench process were other sources of FNs and FPs. In addition, low AF (<5%) considerably influenced the reproducibility and comparability among panels.

CONCLUSIONS

This study provided an integrated practice for developing reference standard to assess oncopanels in detecting somatic mutations and quantitatively revealed the source of detection errors. It will promote optimization, validation, and quality control among laboratories with potential applicability in clinical use.

Early Breast Cancer Evolution by Autosomal Broad Copy Number Alterations

The availability of comprehensive genomic datasets across patient populations enables the application of novel methods for reconstructing tumor evolution within individual patients. To this end, we propose studying autosomal broad copy number alterations (CNAs) as a framework to better understand early tumor evolution. We compared the broad CNAs and somatic mutations of patients with 1 to 10 autosomal broad CNAs against the full set of patients, using data from The Cancer Genome Atlas breast cancer project. We reveal here that the frequency of a chromosome arm obtaining a broad CNA and a genome acquiring somatic mutations changes as autosomal broad CNAs accumulate. Therefore, we propose that the number of autosomal broad CNAs is an important characteristic of breast tumors that needs to be taken into consideration when studying breast tumors. To investigate this idea more in-depth, we next studied the frequency that specific chromosome arms acquire broad CNAs in patients with 1 to 10 broad CNAs. With this process, we identified the broad CNAs that exhibit the fastest rates of accumulation across all patients. This finding suggests a likely order of occurrence of these alterations in patients, which is apparent when we consider a subset of patients with few broad CNAs. Here, we lay the foundation for future studies to build upon our findings and use autosomal broad CNAs as a method to monitor breast tumor progression in vivo to further our understanding of how early tumor evolution unfolds.

Optimized analysis parameters of variant calling for whole genome-based phylogeny of Mycobacteroides abscessus

Whole-genome sequence (WGS) analysis provides the best resolution for reconstructing bacterial phylogeny. However, the resulting tree could vary according to parameters used in the WGS pipeline, making it difficult to compare results across multiple studies. This study compares effects on phylogenies when applying different parameter stringencies. We used as the study model to optimize parameters strains of Mycobacteroides abscessus serially isolated at various intervals, isolates known to represent persistent infection (PI) cases or re-infection (RI) cases and isolates from different subspecies. Un-optimized parameters with low stringency provided an excessive number of SNPs (823) compared to the optimized setting (3 SNPs) between paired strains isolated 1 day apart from PI cases, discordant tree topology and misclassification of subspecies and of instances of RI. We demonstrated that using high-quality variants provides more accuracy for recognizing serial isolates of the same clone versus different clones and for phylogenetic analysis of M. abscessus. Our approach might be used as a model for analyses requiring phylogenetic reconstruction of other bacteria.

SmMIP-tools: a computational toolset for processing and analysis of single-molecule molecular inversion probes-derived data

MOTIVATION

Single-molecule molecular inversion probes (smMIPs) provide an exceptionally cost-effective and modular approach for routine or large-cohort next-generation sequencing. However, processing the derived raw data to generate highly accurate variants calls remains challenging.

RESULTS

We introduce SmMIP-tools, a comprehensive computational method that promotes the detection of single nucleotide variants and short insertions and deletions from smMIP-based sequencing. Our approach delivered near-perfect performance when benchmarked against a set of known mutations in controlled experiments involving DNA dilutions and outperformed other commonly used computational methods for mutation detection. Comparison against clinically approved diagnostic testing of leukaemia patients demonstrated the ability to detect both previously reported variants and a set of pathogenic mutations that did not pass detection by clinical testing. Collectively, our results indicate that increased performance can be achieved when tailoring data processing and analysis to its related technology. The feasibility of using our method in research and clinical settings to benefit from low-cost smMIP technology is demonstrated.

AVAILABILITY AND IMPLEMENTATION

The source code for SmMIP-tools, its manual and additional scripts aimed to foster large-scale data processing and analysis are all available on github (https://github.com/abelson-lab/smMIP-tools). Raw sequencing data generated in this study have been submitted to the European Genome-Phenome Archive (EGA; https://ega-archive.org) and can be accessed under accession number EGAS00001005359.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Clonal and subclonal TP53 molecular impairment is associated with prognosis and progression in multiple myeloma

Aberrations on TP53, either as deletions of chromosome 17p (del17p) or mutations, are associated with poor outcome in multiple myeloma (MM), but conventional detection methods currently in use underestimate their incidence, hindering an optimal risk assessment and prognostication of MM patients. We have investigated the altered status of TP53 gene by SNPs array and sequencing techniques in a homogenous cohort of 143 newly diagnosed MM patients, evaluated both at diagnosis and at first relapse: single-hit on TP53 gene, either deletion or mutation, detected both at clonal and sub-clonal level, had a minor effect on outcomes. Conversely, the coexistence of both TP53 deletion and mutation, which defined the so-called double-hit patients, was associated with the worst clinical outcome (PFS: HR 3.34 [95% CI: 1.37–8.12] p = 0.008; OS: HR 3.47 [95% CI: 1.18–10.24] p = 0.02). Moreover, the analysis of longitudinal samples pointed out that TP53 allelic status might increase during the disease course. Notably, the acquisition of TP53 alterations at relapse dramatically worsened the clinical course of patients. Overall, our analyses showed these techniques to be highly sensitive to identify TP53 aberrations at sub-clonal level, emphasizing the poor prognosis associated with double-hit MM patients.

Halvade somatic: Somatic variant calling with Apache Spark

Background

The accurate detection of somatic variants from sequencing data is of key importance for cancer treatment and research. Somatic variant calling requires a high sequencing depth of the tumor sample, especially when the detection of low-frequency variants is also desired. In turn, this leads to large volumes of raw sequencing data to process and hence, large computational requirements. For example, calling the somatic variants according to the GATK best practices guidelines requires days of computing time for a typical whole-genome sequencing sample.

Findings

We introduce Halvade Somatic, a framework for somatic variant calling from DNA sequencing data that takes advantage of multi-node and/or multi-core compute platforms to reduce runtime. It relies on Apache Spark to provide scalable I/O and to create and manage data streams that are processed on different CPU cores in parallel. Halvade Somatic contains all required steps to process the tumor and matched normal sample according to the GATK best practices recommendations: read alignment (BWA), sorting of reads, preprocessing steps such as marking duplicate reads and base quality score recalibration (GATK), and, finally, calling the somatic variants (Mutect2). Our approach reduces the runtime on a single 36-core node to 19.5 h compared to a runtime of 84.5 h for the original pipeline, a speedup of 4.3 times. Runtime can be further decreased by scaling to multiple nodes, e.g., we observe a runtime of 1.36 h using 16 nodes, an additional speedup of 14.4 times. Halvade Somatic supports variant calling from both whole-genome sequencing and whole-exome sequencing data and also supports Strelka2 as an alternative or complementary variant calling tool. We provide a Docker image to facilitate single-node deployment. Halvade Somatic can be executed on a variety of compute platforms, including Amazon EC2 and Google Cloud.

Conclusions

To our knowledge, Halvade Somatic is the first somatic variant calling pipeline that leverages Big Data processing platforms and provides reliable, scalable performance. Source code is freely available.

Achieving robust somatic mutation detection with deep learning models derived from reference data sets of a cancer sample

BACKGROUND

Accurate detection of somatic mutations is challenging but critical in understanding cancer formation, progression, and treatment. We recently proposed NeuSomatic, the first deep convolutional neural network-based somatic mutation detection approach, and demonstrated performance advantages on in silico data.

RESULTS

In this study, we use the first comprehensive and well-characterized somatic reference data sets from the SEQC2 consortium to investigate best practices for using a deep learning framework in cancer mutation detection. Using the high-confidence somatic mutations established for a cancer cell line by the consortium, we identify the best strategy for building robust models on multiple data sets derived from samples representing real scenarios, for example, a model trained on a combination of real and spike-in mutations had the highest average performance.

CONCLUSIONS

The strategy identified in our study achieved high robustness across multiple sequencing technologies for fresh and FFPE DNA input, varying tumor/normal purities, and different coverages, with significant superiority over conventional detection approaches in general, as well as in challenging situations such as low coverage, low variant allele frequency, DNA damage, and difficult genomic regions.