RNA-seq analysis

BaseClear offers a complete bioinformatics workflow for both prokaryotic and eukaryotic RNA-Seq projects. We use highly cited, continuously supported, and open-source computational tools for read quality control, reference alignment, and (differential) gene expression analysis. Our pipeline has been thoroughly tested using both publicly available RNA-Seq datasets as well as with those that we have generated in-house. In addition, we support you with project-based downstream analyses that may require a tailored approach.

Gene expression analysis

Paired-end reads from the Illumina platform are quality controlled using Trim Galore (Krueger, 2012). Next, QC’d reads are mapped to the reference genome using the STAR RNA-Seq aligner (Dobin et al., 2013). Based on alignments, gene counts are estimated by featureCounts (Liao et al., 2014). In addition to the alignments (BAM files), results per sample are provided as both raw and TPM-normalised counts.

Differential gene expression analysis

In projects that involve samples from different biological conditions, statistical analyses can be used to identify quantitative changes in gene expression between the different conditions. We perform this analysis using the DESeq2 framework (Love et al., 2014). The main output is a table that contains the average expression, fold-change, and associated statistics such as the P and corrected P values for each gene.

Request quotationContact specialist

Differential gene expression analysis

In projects with differential gene expression analysis, we provide you with several data analysis plots.

PCA plot

The PCA plot gives insight into how samples are associated based on their gene expression.

PCA plot

MA plot

The MA-plot gives an overview of the differential expression versus the average expression levels of genes.

MA Plot

Volcano plot

The volcano plot shows the relationship between the statistical significance (P value) and fold-change of gene expression. This allows you to see whether significantly differentially expressed genes also show large differences in expression.

Tailored (project-based) downstream analyses

Our main goal is helping you answer your research question as completely as possible. When results from individual gene expression changes become less straightforward to interpret, this may require the design and application of a tailored approach. In some project this is the simple inclusion of an additional functional database (e.g. eggNOG) in the analysis using which the genes are functionally annotated. In others, it may require a more complex approach that statistical analyses are performed not at gene level but at pathway level. Therefore, in addition to the above-mentioned framework of (differential) gene expression analysis, we are happy to help with you with more downstream analyses, where examples include (but are not limited to) clustering and heatmap generation, gene-set and Gene Ontology enrichment analysis, and pathway visualisation.

 

 

 

 

 

 

Figures: Downstream analyses where gene differences are summarised at higher functional levels can be performed based on the requirements of the project’s specific biological question.

References

  • Dobin, Alexander, et al. “STAR: ultrafast universal RNA-seq aligner.” Bioinformatics 29.1 (2013): 15-21.
  • Krueger, F. “Trim Galore: a wrapper tool around Cutadapt and FastQC to consistently apply quality and adapter trimming to FastQ files, with some extra functionality for MspI-digested RRBS-type (Reduced Representation Bisufite-Seq) libraries.” URL http://www. bioinformatics. babraham. ac. uk/projects/trim_galore/.(Date of access: 28/04/2016) (2012).
  • Liao, Yang, Gordon K. Smyth, and Wei Shi. “featureCounts: an efficient general purpose program for assigning sequence reads to genomic features.” Bioinformatics 30.7 (2013): 923-930.
  • Love, Michael I., Wolfgang Huber, and Simon Anders. “Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2.” Genome biology 15.12 (2014): 550.

 

Get in touch