Quantify QIAseq RNA

The Quantify QIAseq RNA tool can be found in the Toolbox here:

        Toolbox | Biomedical Genomics Analysis (Image biomedical_folder_closed_16_n_p) | QIAseq Tools (Image qiaseq_expert_folder_closed_16_n_p) | Quantify QIAseq RNA (Image rnaseqtrack_16_n_p)

The tool uses RNA-Seq reads as input, together with a reference sequence and target regions that are saved in the CLC_References folder of the Navigation Area when downloading the QIAseq RNA Panels hg38 reference data set.

We recommend using reads that have already been trimmed, as done for example in the Quantify QIAseq RNA Expression template workflow, see Quantify QIAseq RNA Expression. Trimming increases the results accuracy.

The tool performs multiple steps, as described below.

Mapping

The reads are mapped to the target regions using Map Reads to Reference (see http://resources.qiagenbioinformatics.com/manuals/clcgenomicsworkbench/current/index.php?manual=Map_Reads_Reference.html) with default settings. "Auto-detect paired distances" is however disabled, as the automatic estimation of paired distances is not appropriate for targeted data.

The reference sequences of the target regions are created as follows:

Note that all other regions of the genome are ignored.

Filtering

Reads that successfully mapped to the target regions are removed if the UMI does not have the expected length of 12.

Merging

PCR and sequencing errors also happen in the UMIs. To account for this, UMIs can be merged as described in [Peng et al., 2015]. Briefly, UMIs that have a count of at most 2 or lower than 5% of the maximum observed UMI count, are considered for merging. For a UMI to be merged into another UMI, it must differ by at most one base pair and have a count that is at least 6 fold smaller.