DeepVariant is an analysis pipeline that uses a deep neural network to call genetic variants from next-generation DNA sequencing data.
- DeepVariant is highly accurate. In 2016 DeepVariant won PrecisionFDA Truth Challenge in the best SNP Performance category.
- DeepVariant is robust. DeepVariant maintains high accuracy even for error-prone sequencing conditions, including PCR-positive samples and low quality sequencing runs.
- DeepVariant is fast. Using Google Cloud Platform, a whole human genome(1) analysis completes in as little as ~70 minutes(2) and an exome in ~25 minutes.
- DeepVariant is cost-efficient. Using Google Cloud Platform, calling a whole genome costs ~$2-3 and an exome $0.20.
- DeepVariant is flexible. DeepVariant can be easily adjusted or used out of the box for different sequencing technologies and even for non-human species.
- DeepVariant is easy to use. No filtering is needed beyond setting your preferred minimum quality threshold.
- DeepVariant is where you need it. It builds and runs on-premises and in the cloud, as a stand-alone binary or a Docker image, and makes use of hardware accelerators like GPUs and TPUs.
(1): Examples assume 30x coverage. (2): Time estimates do not include mapping.
DeepVariant is a suite of Python/C++ programs that run on any Unix-like operating system. For convenience the documentation refers to building and running DeepVariant on Google Cloud Platform, but the tools themselves can be built and run on any standard Linux computer, including on-premise machines. Note that DeepVariant currently requires Python 2.7 and does not yet work with Python 3.
Pre-built binaries are available at gs://deepvariant/. These are compiled to use SSE4 and AVX instructions, so you'll need a CPU (such as Intel Sandy Bridge) that supports them. (The file /proc/cpuinfo lists these features under "flags".)
Alternatively, see Building and testing DeepVariant for more information on building DeepVariant from sources for your platform.
For managed pipeline execution of DeepVariant see the cost- and speed-optimized, Docker-based pipelines created for Google Cloud Platform.
- Advanced Case Study: Train a customized SNP and small indel variant caller for BGISEQ-500 data
- DeepVariant training data
- DeepVariant usage guide
- Building and testing DeepVariant
- DeepVariant Genomic VCF (gVCF) support
- Getting Started with GCP (It is not required to run DeepVariant on GCP.)
- Colab example: visualizing pileup images/tensors
- Google Developer Codelab: Variant Calling on a Rice genome with DeepVariant
For technical details describing how DeepVariant works please see our Nature Biotechnology publication.
Briefly, we started with some of the reference genomes from Genome in a Bottle, for which there is high-quality ground truth available (or the closest approximation currently possible). Using multiple replicates of these genomes, we produced approximately one hundred million training examples in the form of multi-channel tensors encoding the sequencing instrument data, and then trained a TensorFlow-based image classification model (inception-v3) to assign genotype likelihoods from the experimental data produced by the instrument. Read additional information on the Google Research blog.
Under the hood, DeepVariant relies on Nucleus, a library of Python and C++ code for reading and writing data in common genomics file formats (like SAM and VCF) designed for painless integration with the TensorFlow machine learning framework.
We are delighted to see several external evaluations of the DeepVariant method.
The 2016 PrecisionFDA Truth Challenge, administered by the FDA, assessed several community-submitted variant callsets on the (at the time) blinded evaluation sample, HG002. DeepVariant won the Highest SNP Performance award in the challenge.
DNAnexus posted an extensive evaluation of several variant calling methods, including DeepVariant, using a variety of read sets from HG001, HG002, and HG005. They have also evaluated DeepVariant under a variety of noisy sequencing conditions.
Independent evaluations of DeepVariant v0.6 from both DNAnexus and bcbio are also available. Their analyses support our findings of improved indel accuracy, and also include comparisons to other variant calling tools.
The Genomics team in Google Brain actively supports DeepVariant and are always interested in improving the quality of DeepVariant. If you run into an issue, please report the problem on our Issue tracker. Make sure to add enough detail to your report that we can reproduce the problem and fix it. We encourage including links to snippets of BAM/VCF/etc. files that provoke the bug, if possible. Depending on the severity of the issue we may patch DeepVariant immediately with the fix or roll it into the next release.
If you have questions about next-generation sequencing, bioinformatics, or other general topics not specific to DeepVariant we recommend you post your question to a community discussion forum such as BioStars.
Interested in contributing? See CONTRIBUTING.
DeepVariant is licensed under the terms of the BSD-3-Clause license.
DeepVariant happily makes use of many open source packages. We'd like to specifically call out a few key ones:
-
abseil-cpp and abseil-py
We thank all of the developers and contributors to these packages for their work.
- This is not an official Google product.