Webinar: Trio Analysis with Sentieon Rapid DNAseq on DNAnexus

Please join us Thursday, May 18 at 11am PST (2pm EST) for an educational webinar on trio analysis and cohort joint genotyping using Rapid DNAseq, Sentieon’s variant discovery tool.

Sentieon and DNAnexus are collaborating to make tools available to accelerate genome analysis, while reducing costs for researchers. By leveraging the power and scalability of the cloud,  rapid runtimes for joint genotyping easily scale to thousands of samples. Additionally, using the variant discovery pipelines on DNAnexus, you can obtain identical results to GATK, MuTect, or MuTect2 at a fraction of the cost. Our customers are already experiencing a more than 10-fold increase in processing speed and cost savings of 5-7x.

This webinar will demonstrate the power of running Sentieon Rapid DNAseq on DNAnexus for variant calling of a 30x whole-genome trio in just about one hour. Using the single multi-sample VCF generated from this data, it is possible to identify de novo variants, the parental origin of interesting inherited mutations, and examine the carrier status of individuals for rare recessive mutations.

Webinar Details
Title: Rapid Trio Analysis and Cohort Joint Genotyping with Sentieon on DNAnexus
Speaker: Don Freed, PhD, Bioinformatics Scientist, Sentieon
Date: Thursday, May 18, 2017
Time: 11:00 AM PT, 2:00 PM ET

Attendees will learn how to:

  • Rapidly run joint genotyping and identify de novo mutations from a trio in about one hour
  • Achieve deterministic variant calls, with no run-to-run differences and no down-sampling in high coverage regions
  • Perform accelerated analysis on cohorts of thousands of samples

Upgrading to TLS 1.2

At DNAnexus, security and compliance are paramount. In order to maintain the highest security standards, we occasionally need to update our system protocols. We are in the process of deprecating support for TLS 1.0 and TLS 1.1; starting July 1st 2017, the DNAnexus Platform will only support TLS 1.2. After this date, customers using TLS 1.0 or TLS 1.1 will no longer be able to access the DNAnexus Platform. You may need to upgrade your system to maintain continued access to DNAnexus.

What is TLS?

TLS stands for “Transport Layer Security” and is a protocol that ensures connections made to a remote endpoint are the intended destination through encryption and endpoint identity verification. DNAnexus web and API connections use TLS as a key security component, thus it is important that the latest version of TLS is supported.

What do I need to do? 

For macOS/OS X Users:  

If you are using the DNAnexus Platform SDK (a.k.a. dx-toolkit) with the Python version provided with your operating system, you will need to install an alternative Python version with TLS 1.2 support. We suggest using the Homebrew package manager to install Python version 2.7.13:

  1. Install Homebrew using the instructions at https://brew.sh/.
  2. Once Homebrew is installed, run the following command in your terminal prompt:
    • brew install python
  3. Once the Python version is installed follow the instructions at https://wiki.dnanexus.com/Downloads to download, unpack, and activate the latest dx-toolkit release.

For Internet Explorer Web Browser Users:

If you use Internet Explorer version 10 or earlier, you will need to upgrade your web browser to Internet Explorer 11 or later.

For All Other Users:

If you use a PC, or if your web browser is Internet Explorer 11 or later, Chrome, Firefox, or Safari, we do not expect that your access to  DNAnexus will be impacted by this upgrade.

If you are impacted, please make the necessary modifications by June 30th, 2017 in order to maintain continued access to DNAnexus. Please do not hesitate to contact support@dnanexus.com with any questions or concerns.

Case Study: Calling Somatic Variants with Sentieon on DNAnexus

 Editor’s Note: This blog post was written by Don Freed, Bioinformatics Scientist at Sentieon. Email him at don.freed@sentieon.com.

Sentieon on DNAnexus provides clinicians and researchers an easy-to-use interface to process genomics data with state-of-the art tools in the cloud. In my previous post, we showed how the Sentieon Rapid DNAseq app is capable of processing an entire trio from FASTQ to VCF in about an hour. Today we’ll demonstrate the processing of targeted and exome-capture sequencing of tumor-normal paired samples.

Cancer is a genomic disease, arising as the result of sequential mutations in somatic cells. These mutations change the properties of cancer cells, leading to uncontrolled growth, cell death resistance, and eventually tumor formation and metastasis. Clinicians and researchers are very interested in digging into the genetic changes that give rise to cancer with next-generation sequence (NGS) data. NGS data allows researchers to better understanding the process of cancer and tumor formation and allows clinicians to determine the best course of treatment for cancer patients. As a result, the sequencing of paired tumor-normal samples is a common study design.

Somatic variant calling is a difficult task due to tumor heterogeneity and low abundance of tumor cells.  Two of the more popular tools to address this task are MuTect and MuTect2 developed by the Broad Institute. High coverage across the protein coding sequences of the tumor genome is essential for the sensitive detection of low frequency somatic mosaic mutations. In targeted sequencing experiments, coverage of the targeted regions can exceed 1,000 or even 10,000 reads. Processing such high-coverage datasets is computationally intensive and by default MuTect (1.1.5) and MuTect2 (3.7) downsample the data to 1,000x and 500x, respectively, for improved runtime performance. Unfortunately, downsampling may also result in missed variant calls.

The Sentieon TNseq app on DNAnexus does not perform downsampling, enabling you to achieve the highest sensitivity for low frequency variants present in your sample. TNseq is also deterministic, given the same input data, the same variants will be called. In addition, the app produces results that match MuTect and MuTect2, but with a drastically improved runtime and cost savings.

Running TNseq Tumor-Normal Analysis on DNAnexus

For this case study, we used publicly available tumor-normal paired samples from SRP079186 and whole-genome sequence data from the Texas Cancer Research Biobank (TCRB). According to the project description in SRA, the tumor samples in SRP079186 consist of sequence data from small cell carcinoma mixed with squamous cell carcinoma of the esophagus, while the TCRB data are described in detail in this peer-reviewed publication. We obtained the FASTQ files from the European Nucleotide Archive and directly from the TCRB website and uploaded them into DNAnexus. To process the data, we used the Sentieon TNseq FASTQ to VCF app on DNAnexus. The pipelines produced VCF files of somatic variants called in the exome capture and targeted sequencing data. Alignment, duplicate marking, indel realignment, base-quality score recalibration, indel co-realignment and variant calling were performed with the DNAnexus workflows shown below.

Somatic variant calling on the targeted sequencing sample took only 37 minutes, while somatic variant calling from the paired whole-exome sequenced samples took 10 hours and 18 minutes. Variant calling in the 60x tumor and the 30x normal TCRB whole-genome sequence sample took 16 hours and 45 minutes.

Using Sentieon’s TNseq FASTQ to VCF app, we were able to identify six and ten somatic mosaic mutations in the targeted sequence data, 381 and 739 somatic mosaic mutations in the whole-exome sequence data, and 5,812 and 2,171 somatic mutations in the whole-genome sequence data with TNsnv and TNhaplotyper, respectively. The number of non-silent mutations was 72 and 251 from the whole-exome sequence data as determined by TNsnv and TNhaplotyper, respectively. Other whole-exome sequencing studies of similar cancers have shown that 82 non-silent somatic mutations occur on average in these cancers. This is consistent with our findings given the differences in somatic variant calling pipelines used. Coverage of the targeted regions in the tumor sample exceeded 6,500x at some positions. MuTect2, with default settings, would have discarded more than 92% of the data at these sites, possibly resulting in false-negative or false-positive variant calls.

Digging deeper into the whole-exome data, we find p.R248Q and p.R158H somatic mutations in the gene TP53 encoding the protein p53, perhaps the best studied tumor suppressor protein. p53 functions as a key cell cycle regulator, stopping cell growth in the case of genomic damage. In many cancers p53 is inactivated as a method of bypassing cellular senesce. Germline inactivating mutations in TP53 cause a cancer predisposition syndrome called Li-Fraumeni syndrome. Somatic inactivating mutations in TP53 are expected in this sample; TP53 is the most frequent target of somatic mutation in esophageal squamous cell carcinomas and about half of all cancers have TP53 mutations. If you are interested in learning more about our analysis, you can see our workflow in the “Sentieon FASTQ to VCF (Tumor/Normal)” project on DNAnexus.

With the secure and collaborative cloud-based DNAnexus Platform and Sentieon TNseq FASTQ to VCF app, researchers and clinicians can rapidly and confidently identify somatic mutations from thousands of paired samples. Register here for a license-free trial of Sentieon tools on DNAnexus, available through April 7th.