ALEXEY KOLESNIKOV
Research Areas
Authored Publications
Sort By
Accurate somatic small variant discovery for multiple sequencing technologies with DeepSomatic
Jimin Park
Daniel E. Cook
Lucas Brambrink
Joshua Gardner
Brandy McNulty
Samuel Sacco
Ayse G. Keskus
Asher Bryant
Tanveer Ahmad
Jyoti Shetty
Yongmei Zhao
Bao Tran
Giuseppe Narzisi
Adrienne Helland
Byunggil Yoo
Irina Pushel
Lisa A. Lansdon
Chengpeng Bi
Adam Walter
Margaret Gibson
Tomi Pastinen
Rebecca Reiman
Sharvari Mankame
T. Rhyker Ranallo-Benavidez
Christine Brown
Nicolas Robine
Floris P. Barthel
Midhat S. Farooqi
Karen H. Miga
Andrew Carroll
Mikhail Kolmogorov
Benedict Paten
Kishwar Shafin
Nature Biotechnology (2025)
Preview abstract
Somatic variant detection is an integral part of cancer genomics analysis. While most methods have focused on short-read sequencing, long-read technologies offer potential advantages in repeat mapping and variant phasing. We present DeepSomatic, a deep-learning method for detecting somatic small nucleotide variations and insertions and deletions from both short-read and long-read data. The method has modes for whole-genome and whole-exome sequencing and can run on tumor–normal, tumor-only and formalin-fixed paraffin-embedded samples. To train DeepSomatic and help address the dearth of publicly available training and benchmarking data for somatic variant detection, we generated and make openly available the Cancer Standards Long-read Evaluation (CASTLE) dataset of six matched tumor–normal cell line pairs whole-genome sequenced with Illumina, PacBio HiFi and Oxford Nanopore Technologies, along with benchmark variant sets. Across samples, both cell line and patient-derived, and across short-read and long-read sequencing technologies, DeepSomatic consistently outperforms existing callers.
View details
Accurate human genome analysis with Element Avidity sequencing
Andrew Carroll
Daniel Cook
Lucas Brambrink
Bryan Lajoie
Kelly N. Wiseman
Sophie Billings
Semyon Kruglyak
Bryan R. Lajoie
Junhua Zhao
Shawn E. Levy
Kishwar Shafin
Maria Nattestad
BMC Bioinformatics (2025)
Preview abstract
We investigate the new sequencing technology Avidity from Element Biosciences. We show that Avidity whole genome sequencing matches mapping and variant calling accuracy with Illumina at high coverages (30x-50x) and is noticeably more accurate at lower coverages (20x-30x). We quantify base error rates of Element reads, finding lower error rates, especially in homopolymer and tandem repeat regions. We use Element’s ability to generate paired end sequencing with longer insert sizes than typical short–read sequencing. We show that longer insert sizes result in even higher accuracy, with long insert Element sequencing giving noticeably more accurate genome analyses at all coverages.
View details
Accurate somatic small variant discovery for multiple sequencing technologies with DeepSomatic
Jimin Park
Daniel Cook
Lucas Brambrink
Joshua Gardner
Brandy McNulty
Samuel Sacco
Ayse Keskus
Asher Bryant
Tanveer Ahmad
Jyoti Shetty
Yongmei Zhao
Bao Tran
Giuseppe Narzisi
Adrienne Helland
Byunggil Yoo
Irina Pushel
Lisa Lansdon
Chengpeng Bi
Adam Walter
Margaret Gibson
Tomi Pastinen
Rebecca Reiman
Sharvari Mankame
Rhyker Ranallo-Benavidez
Christine Brown
Nicolas Robine
Floris Barthel
Midhat S. Farooqi
Karen Miga
Andrew Carroll
Mikhail Kolmogorov
Benedict Paten
Kishwar Shafin
Nature Biotechnology (2025)
Preview abstract
Somatic variant detection is an integral part of cancer genomics analysis. While most methods have focused on short-read sequencing, long-read technologies offer potential advantages in repeat mapping and variant phasing. We present DeepSomatic, a deep-learning method for detecting somatic small nucleotide variations and insertions and deletions from both short-read and long-read data. The method has modes for whole-genome and whole-exome sequencing and can run on tumor–normal, tumor-only and formalin-fixed paraffin-embedded samples. To train DeepSomatic and help address the dearth of publicly available training and benchmarking data for somatic variant detection, we generated and make openly available the Cancer Standards Long-read Evaluation (CASTLE) dataset of six matched tumor–normal cell line pairs whole-genome sequenced with Illumina, PacBio HiFi and Oxford Nanopore Technologies, along with benchmark variant sets. Across samples, both cell line and patient-derived, and across short-read and long-read sequencing technologies, DeepSomatic consistently outperforms existing callers.
View details
Highly accurate assembly polishing with DeepPolisher
Mira Mastoras
Mobin Asri
Lucas Brambrink
Prajna Hebbar
Daniel Cook
Maria Nattestad
Julian K. Lucas
Taylor S. Won
Andrew Carroll
Benedict Paten
Kishwar Shafin
Genome Research (2025)
Preview abstract
Accurate genome assemblies are essential for biological research, but even the highest-quality assemblies retain errors caused by the technologies used to construct them. Base-level errors are typically fixed with an additional polishing step that uses reads aligned to the draft assembly to identify necessary edits. However, current methods struggle to find a balance between over- and underpolishing. Here, we present an encoder-only transformer model for assembly polishing called DeepPolisher, which predicts corrections to the underlying sequence using Pacific Biosciences (PacBio) HiFi read alignments to a diploid assembly. Our pipeline introduces a method, PHAsing Reads in Areas Of Homozygosity (PHARAOH), which uses ultralong Oxford Nanopore Technologies (ONT) data to ensure alignments are accurately phased and to correctly introduce heterozygous edits in falsely homozygous regions. We demonstrate that the DeepPolisher pipeline can reduce assembly errors by approximately half, mostly driven by reductions in indel errors. We have applied our DeepPolisher-based pipeline to 180 assemblies from the next Human Pangenome Reference Consortium (HPRC) data release, producing an average predicted quality value (QV) improvement of 3.4 (54% error reduction) for the majority of the genome.
View details
Local read haplotagging enables accurate long-read small variant calling
Daniel Cook
Maria Nattestad
Lucas Brambrink
Brandy McNulty
John E. Gorzynski
Sneha D. Goenka
Euan Ashley
Miten Jain
Karen Miga
Benedict Paten
Andrew Carroll
Kishwar Shafin
Nature Communications (2024)
Preview abstract
Long-read sequencing technology has enabled variant detection in difficult-to-map regions of the genome and enabled rapid genetic diagnosis in clinical settings. Rapidly evolving third-generation sequencing platforms like Pacific Biosciences (PacBio) and Oxford Nanopore Technologies (ONT) are introducing newer platforms and data types. It has been demonstrated that variant calling methods based on deep neural networks can use local haplotyping information with long-reads to improve the genotyping accuracy. However, using local haplotype information creates an overhead as variant calling needs to be performed multiple times which ultimately makes it difficult to extend to new data types and platforms as they get introduced. In this work, we have developed a local haplotype approximate method that enables state-of-the-art variant calling performance with multiple sequencing platforms including PacBio Revio system, ONT R10.4 simplex and duplex data. This addition of local haplotype approximation simplifies long-read variant calling with DeepVariant.
View details
A draft human pangenome reference
Wen-Wei Liao
Mobin Asri
Jana Ebler
Daniel Doerr
Marina Haukness
Shuangjia Lu
Julian K. Lucas
Jean Monlong
Haley J. Abel
Silvia Buonaiuto
Xian Chang
Haoyu Cheng
Justin Chu
Vincenza Colonna
Jordan M. Eizenga
Xiaowen Feng
Christian Fischer
Robert S. Fulton
Shilpa Garg
Cristian Groza
Andrea Guarracino
William T. Harvey
Simon Heumos
Kerstin Howe
Miten Jain
Tsung-Yu Lu
Charles Markello
Fergal J. Martin
Matthew W. Mitchell
Katherine M. Munson
Moses Njagi Mwaniki
Adam M. Novak
Hugh E. Olsen
Trevor Pesout
David Porubsky
Pjotr Prins
Jonas A. Sibbesen
Jouni Sirén
Chad Tomlinson
Flavia Villani
Mitchell R. Vollger
Lucinda L Antonacci-Fulton
Gunjan Baid
Carl A. Baker
Anastasiya Belyaeva
Konstantinos Billis
Andrew Carroll
Sarah Cody
Daniel Cook
Robert M. Cook-Deegan
Omar E. Cornejo
Mark Diekhans
Peter Ebert
Susan Fairley
Olivier Fedrigo
Adam L. Felsenfeld
Giulio Formenti
Adam Frankish
Yan Gao
Nanibaa’ A. Garrison
Carlos Garcia Giron
Richard E. Green
Leanne Haggerty
Kendra Hoekzema
Thibaut Hourlier
Hanlee P. Ji
Eimear E. Kenny
Barbara A. Koenig
Jan O. Korbel
Jennifer Kordosky
Sergey Koren
HoJoon Lee
Alexandra P. Lewis
Hugo Magalhães
Santiago Marco-Sola
Pierre Marijon
Ann McCartney
Jennifer McDaniel
Jacquelyn Mountcastle
Maria Nattestad
Sergey Nurk
Nathan D. Olson
Alice B. Popejoy
Daniela Puiu
Mikko Rautiainen
Allison A. Regier
Arang Rhie
Samuel Sacco
Ashley D. Sanders
Valerie A. Schneider
Baergen I. Schultz
Kishwar Shafin
Michael W. Smith
Heidi J. Sofia
Ahmad N. Abou Tayoun
Francoise Thibauld-Nissen
Francesa Floriana Tricomi
Justin Wagner
Brian Walenz
Jonathan M. D. Wood
Aleksey V. Zimin
Guillaume Borque
Mark J. P. Chaisson
Paul Flicek
Adam M. Phillippy
Justin Zook
Evan E. Eichler
David Haussler
Ting Wang
Erich D. Jarvis
Karen H. Miga
Glenn Hickey
Erik Garrison
Tobias Marschall
Ira M. Hall
Heng Li
Benedict Paten
Nature (2023)
Preview abstract
Here the Human Pangenome Reference Consortium presents a first draft of the human pangenome reference. The pangenome contains 47 phased, diploid assemblies from a cohort of genetically diverse individuals. These assemblies cover more than 99% of the expected sequence in each genome and are more than 99% accurate at the structural and base pair levels. Based on alignments of the assemblies, we generate a draft pangenome that captures known variants and haplotypes and reveals new alleles at structurally complex loci. We also add 119 million base pairs of euchromatic polymorphic sequences and 1,115 gene duplications relative to the existing reference GRCh38. Roughly 90 million of the additional base pairs are derived from structural variation. Using our draft pangenome to analyse short-read data reduced small variant discovery errors by 34% and increased the number of structural variants detected per haplotype by 104% compared with GRCh38-based workflows, which enabled the typing of the vast majority of structural variant alleles per sample.
View details
Improving variant calling using population data and deep learning
Nae-Chyun Chen
Sidharth Goel
Andrew Carroll
BMC Bioinformatics (2023)
Preview abstract
Large-scale population variant data is often used to filter and aid interpretation of variant calls in a single sample. These approaches do not incorporate population information directly into the process of variant calling, and are often limited to filtering which trades recall for precision. In this study, we develop population-aware DeepVariant models with a new channel encoding allele frequencies from the 1000 Genomes Project. This model reduces variant calling errors, improving both precision and recall in single samples, and reduces rare homozygous and pathogenic clinvar calls cohort-wide. We assess the use of population-specific or diverse reference panels, finding the greatest accuracy with diverse panels, suggesting that large, diverse panels are preferable to individual populations, even when the population matches sample ancestry. Finally, we show that this benefit generalizes to samples with different ancestry from the training data even when the ancestry is also excluded from the reference panel.
View details
precisionFDA Truth Challenge V2: Calling variants from short- and long-reads in difficult-to-map regions.
Andrew Carroll
Gunjan Baid
Howard Yang
Maria Nattestad
Sidharth Goel
Cell Genomics (2022)
Preview abstract
The precisionFDA Truth Challenge V2 aimed to assess the state-of-the-art of variant calling in difficult-to-map regions and the Major Histocompatibility Complex (MHC). Starting with fastq files, 20 challenge participants applied their variant calling pipeline and submitted 64 variant callsets for one or more sequencing technologies (~35X Illumina, ~35X PacBio HiFi, and ~50X Oxford Nanopore Technologies). Submissions were evaluated following best practices for benchmarking small variants with the new GIAB benchmark sets and genome stratifications. Challenge submissions included a number of innovative methods for all three technologies, with graph-based methods and machine-learning methods scoring best for short-reads and long-read datasets, respectively. New methods out-performed the winners of the 2016 Truth Challenge across technologies, and new machine-learning approaches combining multiple sequencing technologies performed particularly well. Recent developments in sequencing and variant calling have enabled benchmarking variants in challenging genomic regions, paving the way for the identification of previously unknown clinically relevant variants.
View details
DeepConsensus improves the accuracy of sequences with a gap-aware sequence transformer
Gunjan Baid
Daniel Cook
Kishwar Shafin
Felipe Llinares
Anastasiya Belyaeva
Armin Töpfer
Aaron Wenger
William J. Rowell
Howard Cheng-Hao Yang
Waleed Ammar
Jean-Philippe Vert
Ashish Teku Vaswani
Maria Nattestad
Andrew Walker Carroll
Nature Biotechnology (2022)
Preview abstract
Genomic analysis requires accurate sequencing in sufficient coverage and over difficult genome regions. Through repeated sampling of a circular template, Pacific Biosciences developed long (10-25kb) reads with high overall accuracy, but lower homopolymer accuracy. Here, we introduce DeepConsensus, a transformer-based approach which leverages a unique alignment loss to correct sequencing errors. DeepConsensus reduces errors in PacBio HiFi reads by 42%, compared to the current approach. We show this increases the yield of PacBio HiFi reads at Q20 by 9%, at Q30 by 27%, and at Q40 by 90%. With two SMRT cells of HG003, reads from DeepConsensus improve hifiasm assembly contiguity (NG50 4.9Mb to 17.2Mb), increase gene completeness (94% to 97%), reduce false gene duplication rate (1.1% to 0.5%), and improve assembly base accuracy (QV43 to QV45), and also reduce variant calling errors by 24%.
View details
Ultra-rapid whole genome nanopore sequencing in a critical care setting
Andrew Carroll
Ankit Sethia
Benedict Paten
Christopher Wright
Courtney J. Wusthoff
Daniel R Garalde
Dianna G. Fisk
Elizabeth Spiteri
Euan Ashley
Fritz J. Sedlazeck
Gunjan Baid
Henry Chubb
Jeffrey W Christle
Jeffrey W. Christle
John E. Gorzynski
Jonathan A Bernstein
Joseph Guillory
Joshua W. Knowles
Katherine Xiong
Kishwar Shafin
Kyla Dunn
Marco Perez
Maria Nattestad
Maura RZ Ruzhnikov
Megan E. Grove
Mehrzad Samadi
Michael Ma
Miten Jain
Scott R. Ceresnak
Sneha D. Goenka
Tanner D. Jensen
Tia Moscarello
Tong Zhu
Trevor Pesout
New England Journal of Medicine (2022)
Preview abstract
Background
Genetic disease is a major contributor to critical care hospitalization, especially in younger patients. While early genetic diagnosis can guide clinical management, the turnaround time for whole genome based diagnostic testing has traditionally been measured in months. Recent programs in neonatal populations have reduced turnaround time into the range of days and shown that rapid genetic diagnosis enhances patient care and reduces healthcare costs. Yet, most decisions in critical care need to be made on hourly timescales.
Methods
We developed a whole genome sequencing approach designed to provide a genetic diagnosis within hours. Optimized highly parallel nanopore sequencing was coupled to a high-performance cloud compute system to implement near real-time basecalling and alignment followed by accelerated central and graphics processor unit variant calling. A custom scheme for variant prioritization took only minutes to rank variants most likely to be deleterious allowing efficient manual review and classification according to American College of Medical Genetics and Genomics guidelines.
Results
We performed whole genome sequencing on 12 patients from the critical care units of Stanford hospitals. In 10 cases, the pipeline produced diagnostic results faster than all previously published clinical genome analyses. Per patient, DNA extraction, library preparation, and nanopore sequencing across 48 flow cells generated 173–236 GigaBases of sequencing data in as little as 1:50 hours. After optimization, the average turnaround time was 7:58 hours (range 7:18–9:0 hours). A pathogenic or likely pathogenic variant was identified in five out of 12 patients (42%). After Sanger or short read sequencing confirmation in a CLIA-approved laboratory, this validated diagnosis altered clinical management in every case.
Conclusions
We developed an approach to make a genetic diagnosis from whole genome sequencing in hours, returning actionable, cost-saving diagnostic information on critical care timescales.
View details
