BOL: Related items

RNA-Seq Data Pathway and Gene-set Analysis Workflows

Jit — Fri, 25 Oct 2013 08:00:48 -0500

It describe the GAGE (Luo et al., 2009) /Pahview (Luo and Brouwer, 2013) workflows on RNA-Seq data pathway analysis and gene-set analysis. The gage package (2.12.0) now includes a new tutorial, “RNA-Seq Data Pathway and Gene-set Analysis Workflows“.

First cover a full workflow from preparation, reads counting, data preprocessing, gene set test, to pathway visualization in about 40 lines of codes. The same workflow can be used for GO analysis or other types of gene set analysis too. We also describe joint workflows, i.e. to do gene-level analysis using one of the major RNA-Seq analysis tools, DEseq/DEseq2, edgeR, limma and Cufflinks, and feed the results into GAGE/Pahview for pathway analysis or visualization. All these workflows are implemented in R/Bioconductor.

The work ows cover the most common situations and issues for RNA-Seq data pathway analysis. Issues like data quality assessment are relevant for data analysis in general yet out the scope of this tutorial. Although we focus on RNA-Seq data here, but pathway analysis work ow remains similar for microarray, particularly step 3-4 would be the same. Please check gage and pathview vigenttes for details.

Note: You need to update to current release versions of R(3.0.2)/ Bioconductor(2.13) to use all the features.

Reference:

Please check it out:
http://bioconductor.org/packages/release/bioc/html/gage.html
http://bioconductor.org/packages/release/bioc/vignettes/gage/inst/doc/RNA-seqWorkflow.pdf

Scallop: reference-based transcriptome assembler for RNA-seq

Rahul Nayak — Tue, 08 May 2018 04:23:27 -0500

Scallop is an accurate reference-based transcript assembler. Scallop features its high accuracy in assembling multi-exon transcripts as well as lowly expressed transcripts. Scallop achieves this improvement through a novel algorithm that can be proved preserving all phasing paths from reads and paired-end reads, while also achieves both transcripts parsimony and coverage deviation minimization.

Scallop paper has been published at Nature Biotechnology. The datasets and scripts used in this paper to compare the performance of Scallop and other assemblers are available at scalloptest.

Please also checkout the podcast about Scallop (thanks Roman Cheplyaka for the interview). It is available at both the bioinformatics chat and iTunes.

https://github.com/Kingsford-Group/scallop

Address of the bookmark: https://github.com/Kingsford-Group/scallop

BioJupies: Automatically Generates RNA-seq Data Analysis Notebooks

Rahul Nayak — Sun, 20 Dec 2020 11:43:45 -0600

With BioJupies you can produce in seconds a customized, reusable, and interactive report from your own raw or processed RNA-seq data through a simple user interface

BioJupies now supports user accounts! Sign in from the top right corner of the page for access to unlimited private notebooks, RNA-seq datasets and alignment jobs.

Address of the bookmark: https://amp.pharm.mssm.edu/biojupies/

RNA-seq workflow: gene-level exploratory analysis and differential expression

Jit — Sat, 09 Oct 2021 07:59:23 -0500

Here we walk through an end-to-end gene-level RNA-seq differential expression workflow using Bioconductor packages. We will start from the FASTQ files, show how these were quantified to the reference transcripts, and prepare gene-level count datasets for downstream analysis. We will perform exploratory data analysis (EDA) for quality assessment and to explore the relationship between samples, perform differential gene expression analysis, and visually explore the results.

Address of the bookmark: http://master.bioconductor.org/packages/release/workflows/vignettes/rnaseqGene/inst/doc/rnaseqGene.html

Pathway Analysis

Rahul Agarwal — Fri, 03 Oct 2014 08:51:13 -0500

Pathway Analysis is usually performed with aim to enrich the genes with their functional information and reveal the underlying biological mechanisms pursue by genes. Pathway Analysis is not only limited to what biological pathways a particular set of expressed genes follow but also to disclose the relationships between these genes. With availability of more genomics, transcriptomics and proteomics data, interactions between genes involve in multiple pathways become more clear and also relationships between the genes, their transcripts, and their gene products. However, existing tools and dbs mainly based on knowledge driven approach in which pathways will be identified by finding the correlation between the information in one of the pathway knowledge databases (KEGG,Reactome,Panther,BioCarta, Panther,GO,NCI,WikiPathways,etc) and gene expression result for a specific conditions for instance tumor, obesity , cold resistant crops/plants, etc.

Introductory Articles/ppt/sources:

http://www.ploscompbiol.org/article/info%3Adoi%2F10.1371%2Fjournal.pcbi.1002375

http://bioinformatics.mdanderson.org/MicroarrayCourse/Lectures09/Pathway%20Analysis.pdf

http://gettinggeneticsdone.blogspot.de/2012/03/pathway-analysis-for-high-throughput.html

http://davetang.org/muse/tag/pathway/

https://www.biostars.org/p/42219/

http://bioinformatics.ca//files/public/Pathways_2014_Module4_v2.pdf

http://bioinformatics.ca//files/public/Pathways_2014_Module2.pdf

Impotant Database and Tools:

GeneMANIA, Cytoscape, IPA and Metacore (Commerical ), Pathway Commons, Reactome ,Panther, BioCyc, WikiPathways, Pathvisio, KEGG, NCI, Stringdb, Amigo, WebGestalt ,ConsensusPathDB ,GSEA,Blast2go

Popular R based tools:

Reactome.db, ReactomePA, ClusterProfiler, Gage, SPIA, topGO, Pathview,DOSE,GOStat

More:

http://www.bioconductor.org/help/search/index.html?q=Enrichment+analysis+

RNA-Seq differential expression work flow using DESeq2

BioStar — Mon, 23 Aug 2021 10:57:14 -0500

One of the aim of RNAseq data analysis is the detection of differentially expressed genes. The package DESeq2 provides methods to test for differential expression analysis.

Address of the bookmark: http://www.sthda.com/english/wiki/rna-seq-differential-expression-work-flow-using-deseq2

Kallisto vs Salmon: Choosing the Right Tool for RNA-Seq Quantification

BioStar — Fri, 02 May 2025 06:28:46 -0500

In the world of transcriptomics, quantifying gene and transcript expression accurately and efficiently is crucial. With the explosion of RNA-Seq data, researchers have turned to fast, alignment-free tools that streamline the quantification process without compromising accuracy. Two leading tools in this space are Kallisto and Salmon. Both tools are highly efficient and widely used in the bioinformatics community, but they differ in subtle yet important ways. If you're unsure which one to use for your next RNA-Seq project, this post is for you.

What Are Kallisto and Salmon?

At their core, both Kallisto and Salmon are tools for quantifying transcript abundance from RNA-Seq reads. They bypass traditional alignment-based methods, replacing them with pseudoalignment or quasi-mapping, which drastically speeds up the process.

Kallisto was developed by Lior Pachter’s lab and introduced the concept of pseudoalignment using a de Bruijn graph.
Salmon, developed by Rob Patro’s group, builds on this idea with quasi-mapping and offers additional features like advanced bias correction.

Head-to-Head Comparison

1. Algorithm

Kallisto uses pseudoalignment, focusing on matching k-mers from reads to a transcriptome index.
Salmon uses quasi-mapping, which adds more flexibility and can also work with aligned reads (BAM files).

2. Input and Flexibility

Kallisto works with raw FASTQ reads and requires a custom transcriptome index.
Salmon accepts FASTQ or pre-aligned BAM files, giving you more workflow options.

3. Bias Correction

One of Salmon’s major advantages is its sophisticated bias correction system. It corrects for:

Sequence-specific bias
Positional bias
GC-content bias

Kallisto offers basic sequence bias correction but lacks the comprehensive models found in Salmon.

4. Speed and Resources

Kallisto is blazing fast and slightly more memory-efficient.
Salmon is still very fast, but the added features can come at a small computational cost.

5. Output and Downstream Analysis

Both tools provide transcript-level quantifications and support bootstrapping for variance estimation.
Salmon can also summarize counts at the gene level if provided with a mapping file (--geneMap).
Kallisto integrates seamlessly with Sleuth for differential expression analysis.
Salmon works well with tximport, DESeq2, edgeR, and other Bioconductor tools.

Choosing the Right Tool

Goal	Recommended Tool
Maximum speed	Kallisto
Advanced bias correction	Salmon
Use BAM files	Salmon
Transcript-level quantification with Sleuth	Kallisto
Integration with DESeq2/edgeR	Salmon

Example Command Lines

Kallisto (paired-end):

kallisto quant -i transcriptome.idx -o output -b 100 sample_R1.fastq sample_R2.fastq

Salmon (paired-end, bias correction):

salmon quant -i salmon_index -l A -1 sample_R1.fastq -2 sample_R2.fastq \
  -p 8 --validateMappings --seqBias --gcBias -o output

Conclusion

Both Kallisto and Salmon are exceptional tools that have transformed RNA-Seq analysis. Your choice largely depends on your priorities—whether it's speed, accuracy, flexibility, or compatibility with downstream tools.

For many users, Salmon offers a more complete and flexible solution, especially when bias correction and gene-level outputs are essential. However, Kallisto remains a favorite for quick, accurate quantification, especially when paired with the Sleuth pipeline.

Will MinION Nanopore sequencing increase the number of Next Generation Sequencing projects?

Strand — Tue, 04 Aug 2015 05:14:07 -0500

Will MinION Nanopore sequencing increase the number of Next Generation Sequencing projects?

Webinar: Wednesday 21 August 2013 at Noon EDT

Jitendra Narayan — Sun, 11 Aug 2013 19:31:56 -0500

This webinar will describe the use of combinatorial pooling to reconstruct gene sequences within BACs. Recent work in barley has shown that this level of sequence knowledge is sufficient to support critical end-point objectives such as map-based cloning and marker-assisted breeding.

http://www.extension.org/pages/67926/upcoming-webinar:-selective-sequencing-through-combinatorial-pooling#.UggsVuHyPqU

Surrogate Variable Analysis (SVA)

Jit — Thu, 30 Oct 2014 08:01:58 -0500

The sva package contains functions for removing batch effects and other unwanted variation in high-throughput experiment. Specifically, the sva package contains functions for the identifying and building surrogate variables for high-dimensional data sets. Surrogate variables are covariates constructed directly from high-dimensional data (like gene expression/RNA sequencing/methylation/brain imaging data) that can be used in subsequent analyses to adjust for unknown, unmodeled, or latent sources of noise. The sva package can be used to remove artifacts in three ways:

(1) identifying and estimating surrogate variables for unknown sources of variation in high-throughput experiments (Leek and Storey 2007 PLoS Genetics,2008 PNAS),

(2) directly removing known batch effects using ComBat (Johnson et al. 2007 Biostatistics) and

(3) removing batch effects with known control probes (Leek 2014 biorXiv).

Removing batch effects and using surrogate variables in differential expression analysis have been shown to reduce dependence, stabilize error rate estimates, and improve reproducibility, see (Leek and Storey 2007 PLoS Genetics, 2008 PNAS or Leek et al. 2011 Nat. Reviews Genetics).

More at http://www.bioconductor.org/packages/release/bioc/html/sva.html