Identify Variants and Add Expression Values
The Identify Variants and Add Expression Values workflow can be used to identify novel and known mutations in RNA-Seq data, automatically map, quantify, and annotate the transcriptomes, and compare the mutational patterns in the samples with the expression values of the corresponding transcripts and genes.
To run the Template workflow:
Toolbox | Template Workflows | Biomedical Workflows () | Whole Transcriptome Sequencing () | (Human, Mouse or Rat) | Identify Variants and Add Expression Values ()
- Double-click on the Identify Variants and Add Expression Values tool to start the analysis. If you are connected to a server, you will first be asked where you would like to run the analysis.
- Specify the RNA-Seq reads to analyze. The reads can be selected by double-clicking on the reads file name or clicking once on the file and then clicking on the arrow pointing to the right side in the middle of the wizard (figure 23.26).
Figure 23.26: Select the sequencing reads to analyze.
- Select the Reference Data Set that is relevant to your study (figure 23.27).
Figure 23.27: Select the relevant data set for the samples being studied.
- Configure the parameters for the RNA-Seq Analysis (figure 23.28).
Figure 23.28: Configure the RNA-Seq Analysis.
If you wish to use spike-in controls, add the relevant file in the "Spike-in controls" field.
You can also specify that the reads should be mapped only in their forward or reverse orientation (it is by default set to both). Choosing to restrict mapping to one direction is typically appropriate when a strand specific protocol for read generation has been used, as it allows assignment of the reads to the right gene in cases where overlapping genes are located on different strands. Also, applying the 'strand specific' 'reverse' option in an RNA-Seq run could allow the user to assess the degree of antisense transcription. Note that mate pairs are not supported when choosing the forward only or reverse only option.
- Specify a target region for the Indels and Structural Variants tool (figure 23.29).
Figure 23.29: Specify the target region for the Indels and Structural Variants tool.
The targeted region file is a file that specifies which regions have been sequenced. This file is something that you must provide yourself, as this file depends on the technology used for sequencing. You can obtain the targeted regions file from the vendor of your targeted sequencing reagents. Remember that you have a hg38-specific BED file when using hg38 as reference, and hg19-specific BED file when using hg19 as reference.
- Set the parameters for the Low Frequency Variant Detection step (see figure 23.30). For a description of the different parameters that can be adjusted in the variant detection step, see http://resources.qiagenbioinformatics.com/manuals/clcgenomicsworkbench/current/index.php?manual=Low_Frequency_Variant_Detection.html.
Figure 23.30: Specify the parametes for transcriptomic variant detection.
- If you are working with the workflow for Human, specify here the relevant 1000 Genomes population (and HapMap populations at the next step) (see figure 23.31). Choose the population that matches best the population your samples are derived from.
Figure 23.31: Select the relevant population from the drop-down list for Hapmap databases.
Under "Locked settings" you can see that "Automatically join adjacent MNVs and SNVs" has been selected. The reason for this is that many databases do not report a succession of SNVs as one MNV as is the case for CLC Workbench, and as a consequence it is not possible to directly compare variants called with CLC Workbench with these databases. In order to support filtering against these databases anyway, the option to Automatically join adjacent MNVs and SNVs is enabled. This means that an MNV in the experimental data will get an exact match, if a set of SNVs and MNVs in the database can be combined to provide the same allele. This assumes that SNVs and MNVs in the track of known variants represent the same allele, although there is no evidence for this in the track of known variants.
- 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.
- Choose to Save your results and click on the button labeled Finish.
The following outputs are generated:
- Gene expression () A track showing gene expression annotations. Hold the mouse over or right-click on the track: a tooltip will appear with information about e.g. gene name and expression values.
- Transcript expression () A track showing transcript expression annotations.
- RNA-Seq Mapping Report () This report contains information about the reads, reference, transcripts, and statistics. This is explained in more details here: (http://resources.qiagenbioinformatics.com/manuals/clcgenomicsworkbench/current/index.php?manual=RNA_Seq_report.html).
- Read Mapping () The mapped RNA-Seq reads. The RNA-Seq 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).
- Annotated Variants with Expression Values () Annotation track showing the variants. Hold the mouse over one of the variants or right-clicking on the variant. A tooltip will appear with detailed information about the variant.
- RNA-Seq Track List () A collection of tracks presented together. Shows the annotated variant track together with the human reference sequence, genes, transcripts, coding regions, and variants detected in ClinVar and dbSNP (see figure 23.16).
- Log () A log of the workflow execution.