Subsections

Identify Somatic Variants from Tumor Normal Pair (WES)

The Identify Somatic Variants from Tumor Normal Pair (WES) ready-to-use workflow can be used to identify potential somatic variants in a tumor sample when you also have a normal/control sample from the same individual.

When running this workflow the reads are mapped and the variants identified. An internal workflow removes germline variants that are found in the mapped reads of the normal/control sample and variants outside the target region are removed as they are likely to be false positives due to non-specific mapping of sequencing reads. Next, remaining variants are annotated with gene names, amino acid changes, conservation scores and information from relevant databases like ClinVar (variants with clinically relevant association). Finally, information from dbSNP is added to see which of the detected variants have been observed before and which are completely new.

Before starting the workflow, you will need to import in the workbench a file with the genomic regions targeted by the amplicon or hybridization kit. Such a file (a BED or GFF file) is usually available from the vendor of the enrichment kit and sequencing machine. Use the Import | Tracks tool to import it in your Navigation Area.

Run the Identify Somatic Variants from Tumor Normal Pair (WES) workflow

To run the Identify Somatic Variants from Tumor Normal Pair (WES) tool, go to:

        Toolbox | Ready-to-Use Workflows | Whole Exome Sequencing (Image exome_sequencing_closed_16_n_p) | Somatic Cancer (Image somatic_folder_closed_16_n_p) | Identify Somatic Variants from Tumor Normal Pair (WES) (Image filter_somatic_var_wes_16_n_p)

  1. Go to the toolbox and double-click on the Identify Somatic Variants from Tumor Normal Pair (WES) ready-to-use workflow.

  2. First (figure 17.21), select the tumor sample reads.

    Image filter_somatic_variants_from_tumor_normal_step1_wes
    Figure 17.21: Select the tumor sample reads.

  3. In the next wizard step, specify the normal sample reads.

  4. In the next dialog, select which data set should be used to identify variants (figure 17.22).

    Image identify_somatic_variants_wes
    Figure 17.22: Choose the relevant reference Data Set to identify variants.

  5. The following 2 steps allow you to restrict the calling of Indels and Structural Variants to the targeted regions, both for tumor and normal reads (figure 17.23).

    Image filter_somatic_variants_from_tumor_normal_step3_wes
    Figure 17.23: Specify the target regions track.

  6. Set the parameters for the Low Frequency Variant Detection step (figure 17.24).

    Image filter_somatic_variants_from_tumor_normal_step5_wes
    Figure 17.24: Specify the settings for the variant detection.

    For a description of the different parameters that can be adjusted, see http://resources.qiagenbioinformatics.com/manuals/clcgenomicsworkbench/current/index.php?manual=Low_Frequency_Variant_Detection.html. If you click on "Locked Settings", you will be able to see all parameters used for variant detection in the ready-to-use workflow.

  7. In the following 2 wizard steps, you can select your target regions track to be used for reporting the performance of the targeted re-sequencing experiment for the tumor and normal samples successively (figure 17.25). The targeted region track should be the same as the track you selected in the previous wizard steps. Variants found outside the targeted regions will not be included in the output that is generated with the ready-to-use workflow.

    Image filter_somatic_variants_from_tumor_normal_step4_wes
    Figure 17.25: Select your target region track.

    For a description of the different parameters that can be adjusted, see http://resources.qiagenbioinformatics.com/manuals/clcgenomicsworkbench/current/index.php?manual=QC_Targeted_Sequencing.html. If you click on "Locked Settings", you will be able to see all parameters used for the QC for Targeted Sequencing tool in the ready-to-use workflow.

  8. In the Remove Variants Present in Control Reads step, you can adjust the settings for removal of germline variants (figure 17.26).

    Image filter_somatic_variants_from_tumor_normal_step6_wes
    Figure 17.26: Specify setting for removal of germline variants.

  9. In the last wizard step you can check the selected settings by clicking on the button labeled Preview All Parameters. In the Preview All Parameters wizard you can only check the settings, and if you wish to make changes you have to use the Previous button from the wizard to edit parameters in the relevant windows.

  10. Choose to Save your results and click Finish.

Output from the Identify Somatic Variants from Tumor Normal Pair (WES) workflow

The following outputs are generated:

  1. Read Mapping Normal (Image read_track_16_n_p) The mapped sequencing reads for the normal sample. The reads are shown in different colors depending on their orientation, whether they are single reads or paired reads, and whether they map unambiguously (see http://resources.qiagenbioinformatics.com/manuals/clcgenomicsworkbench/current/index.php?manual=Coloring_mapped_reads.html).
  2. Read Mapping Tumor (Image read_track_16_n_p) The mapped sequencing reads for the tumor sample.
  3. Target Region Coverage Report Normal (Image proteinreport_16_n_p) The report consists of a number of tables and graphs that in different ways provide information about the mapped reads from the normal sample.
  4. Target Region Coverage Tumor (Image annotation_track_16_n_p) A track showing the targeted regions. The table view provides information about the targeted regions such as target region length, coverage, regions without coverage, and GC content.
  5. Target Region Coverage Report Tumor (Image proteinreport_16_n_p) The report consists of a number of tables and graphs that in different ways provide information about the mapped reads from the tumor sample.
  6. Amino Acids Changes Track that shows the consequences of the variants at the amino acid level in the context of the original amino acid sequence. A variant introducing a stop mutation is illustrated with a red amino acid.
  7. Annotated Somatic Variants (Image variant_track_16_n_p) A variant track holding the identified and annotated somatic variants. The variants can be shown in track format or in table format. When holding the mouse over the detected variants in the Track List, a tooltip appears with information about the individual variants. You will have to zoom in on the variants to be able to see the detailed tooltip.
  8. Track List Tumor Normal Comparison (Image trackset_16_n_p) A collection of tracks presented together. Shows the annotated variant track together with the human reference sequence, genes, transcripts, coding regions, the mapped reads for both normal and tumor, the annotated somatic variants, information from the ClinVar database, and finally a track showing the conservation score (see figure 17.27).

Image identify_somatic_variants_genomebrowserview_wes
Figure 17.27: The Track List presents all the different data tracks together and makes it easy to compare different tracks.