mRNA Differential Gene Expression workflow enables faster and quantitative analysis of differentially expressed genes between your samples from your RNA-seq data. Downstream analysis such as GO and KEGG-pathway analysis can be performed on the set of interesting genes. In this workflow, parameters can also be optimized depending on datasets and analysis needs. State-of-the-art publication-quality visualization and reports of your analysis results can be easily downloaded in the form of tables, images and html.

Try it in Genome Explorer®

Pipeline

Import Disease/Normal Sample (FASTQ | SRA) Cleaned Disease/Normal (FASTQ) QC Reports Map to Reference (Disease/Normal BAM) QC Reports BAM Processing (Cleaned Disease/Normal BAM) QC Reports Estimate Transcript GTF Compare Transcript GTF Merge Transcript GTF Differential Expression Project Reports

Import Samples

Samples are downloaded from amazon s3 bucket.

Pre-Alignment QC

Pre-alignment quality control step is performed to ensure the quality of raw data and it also provides a QC report which can spot biases in the starting library material or during sequencing. We use the FASTQC program which compiles important statistics on NGS data, such as per-base and per-sequence quality scores, GC content and checks for duplicated and overrepresented sequences within the file.

Alignment

In this step sequencing reads are aligned to the reference genome. The mapping is performed using HISAT (hierarchical indexing for spliced alignment of transcripts) which uses an indexing scheme based on the Burrows-Wheeler transform and the Ferragina-Manzini (FM) index. Read alignment information is reported in SAM/BAM format that is used for calculating depth of coverage and calling variants.

Alignment QC

Alignment QC step examines the sequencing alignment data in SAM/BAM file according to the features of the mapped reads and their genomic properties and provides an overall view of the data that helps to detect biases in the sequencing and/or mapping of the data which you should be aware of before doing any further analysis. All statistics about sufficient coverage, depth of sequencing and read counts are provided at this step.

Estimate transcript

The objective is to assemble reads that were previously aligned to the reference genome into transcript structure (transcript assembly process) and to estimate the abundance levels of genes and transcripts. The expression levels are reported in terms of Fragments Per Kilobase of Exon per Million mapped fragments (FPKM) units.

Compare transcript

The expression values of genes and transcripts can be compared across samples post transcript-assembly process.

Merge transcript

In this objective, GTF files from each sample are merged together into a single GTF file, which is used as an input to Differential Expression objective.

Differential expression

The objective is to find significant changes in gene/transcript expression, splicing, and promoters between two samples in a group.

Supported Organisms

1

Homo sapiens

2

Mus musculus

Reports

Turn your complex genomic data into a meaningful report of actionable information and insightful visualizations.
Reports are highly customizable and exportable in industry-class formats.

Splice Junction

Displaying Putative Splice Sites Around Donor

Junction Analysis

Displaying AT-AS Events

FPKM Distribution

Gene density plot

volcano plot

Read Duplication Summary

Reads Genomic Origin

Pre-filter Coverage Profile (Low)