Nicolas L Bray, Harold Pimentel, Páll Melsted and Lior Pachter, Near-optimal probabilistic RNA-seq quantification, Nature Biotechnology 34, 525–527 (2016), doi:10.1038/nbt.3519

Address of the bookmark: https://pachterlab.github.io/kallisto/about

What is RNA-Seq?

RNA-Seq is a next-generation sequencing (NGS) technology used to study the transcriptome—the complete set of RNA molecules in a cell. It quantifies gene expression, detects novel transcripts, and captures alternative splicing events with high sensitivity and resolution.

Workflow for RNA-Seq Analysis

RNA-Seq analysis involves several stages, each requiring computational tools and expertise.

1. Experimental Design and Data Acquisition

Before diving into analysis, bioinformaticians should consider:

• Biological Replicates: Ensure statistical power to detect meaningful differences.
• Sequencing Depth: Align sequencing depth to study objectives (e.g., higher depth for low-abundance transcripts).
• Paired-End vs. Single-End: Paired-end sequencing provides more detailed information on transcript structure.

Once sequencing is complete, raw data is provided in FASTQ format, containing sequence reads and quality scores.

2. Quality Control and Preprocessing

Quality control (QC) ensures data integrity. Tools such as FastQC evaluate metrics like base quality, GC content, and adapter contamination.

Preprocessing Steps:

• Trimming: Tools like Trimmomatic or Cutadapt remove low-quality bases and adapter sequences.
• Filtering: Discard reads below a certain quality threshold or length.

3. Read Alignment

Reads are mapped to a reference genome or transcriptome to determine their origin. Alignment tools include:

• HISAT2: Handles large genomes efficiently and supports spliced alignments.
• STAR: High-speed aligner optimized for RNA-Seq.
• Bowtie2: Suitable for short-read alignment.

Output: A SAM/BAM file containing aligned reads.

4. Transcript Assembly and Quantification

This step involves identifying transcripts and quantifying their expression levels. Tools used include:

• StringTie: Assembles and quantifies transcripts from aligned reads.
• Salmon/Kallisto: Perform pseudo-alignment for rapid and accurate quantification.

Expression levels are typically measured as TPM (transcripts per million) or FPKM (fragments per kilobase of transcript per million mapped reads).

5. Differential Expression Analysis

To identify genes with altered expression between conditions, bioinformaticians use tools such as:

• DESeq2: Accounts for data normalization and variability.
• edgeR: Handles overdispersed count data efficiently.
• Limma-voom: Combines linear modeling with RNA-Seq count data.

The output includes a list of differentially expressed genes (DEGs) with statistical significance and fold-change values.

6. Functional Annotation and Pathway Analysis

Understanding the biological significance of DEGs involves:

• Gene Ontology (GO) Analysis: Tools like DAVID or clusterProfiler categorize genes based on their biological functions.
• Pathway Enrichment Analysis: Identifies pathways enriched in DEGs using tools like KEGG, Reactome, or GSEA.

7. Visualization

Visualizing results enhances interpretability. Common visualizations include:

• Heatmaps: Show expression patterns across samples (e.g., pheatmap).
• Volcano Plots: Highlight significant DEGs (e.g., ggplot2).
• PCA/UMAP: Assess sample clustering and variability (e.g., Seurat).

Challenges in RNA-Seq Analysis

1. Batch Effects: Technical variability can confound biological signals. Combat this with normalization techniques or batch-correction tools like ComBat.
2. Low-Quality Samples: Poor-quality RNA impacts downstream analyses.
3. Computational Complexity: RNA-Seq generates massive datasets, requiring robust computing resources and optimized pipelines.

Key Tools and Resources

• Bioconductor: A treasure trove of R packages for RNA-Seq analysis.
• Galaxy: A web-based platform for running RNA-Seq workflows.
• Nextflow/Snakemake: Workflow management tools to streamline analyses.

Applications of RNA-Seq

RNA-Seq is used in diverse research areas, including:

• Cancer Transcriptomics: Identifying tumor-specific expression profiles.
• Developmental Biology: Studying dynamic transcriptome changes.
• Drug Discovery: Screening genes modulated by therapeutic compounds.

Conclusion

RNA-Seq analysis is a cornerstone of modern transcriptomics, offering bioinformaticians a versatile toolkit for unraveling gene expression and regulation. Mastering RNA-Seq workflows and tools empowers researchers to transform raw sequencing data into biological discoveries.

Whether you’re investigating disease mechanisms, exploring cellular pathways, or developing new therapeutics, RNA-Seq is a powerful ally in your bioinformatics arsenal.

High affinity antisense oligonucleotides (ASOs) containing bicylic modifications (BNA) such as locked nucleic acid (LNA) or constrained ethyl (cEt) designed to induce target RNA cleavage have been shown to have enhanced potency along with a higher propensity to cause hepatotoxicity. In order to unravel the mechanism of this hepatotoxicity, we leveraged GeneSpring GX analysis software to analyze transcriptional profiles from the livers of mice treated with a panel of highly efficacious hepatotoxic or non-hepatotoxic LNA ASOs.

Speaker:
Sebastien A. Burel, PhD
Director, Nonclinical Development, Ionis Pharmaceuticals, California

Details:
14 June 2017, 8 AM PST

Register for this Webinar

“As we roll out genomic medicine we are fighting against this society-wide misconception that having the bad gene means you’re going to get the disease. That’s only true in a very few cases.”

Source:Today Health

Address of the bookmark: http://www.today.com/health/should-you-get-sequenced-not-all-bad-genes-predict-disease-8C11017154

Title: Comparison of Short Read De Novo Alignment Algorithms 

Author: Nikhil Gopal

Address of the bookmark: http://biochem218.stanford.edu/Projects%202011/Gopal%202011.pdf

http://bmi.osu.edu/hpc/software/benchmark/

http://bmi.osu.edu/hpc/software/pmap/pmap.html

Other similiar papers:

http://online.liebertpub.com/doi/pdf/10.1089/cmb.2012.0022

http://bioinformatics.oxfordjournals.org/content/28/24/3169

Some new Mapping tool links:

GSNAP

http://research-pub.gene.com/gmap/

RMAP

http://rulai.cshl.edu/rmap/

mrsFAST

http://mrsfast.sourceforge.net/Home

http://sourceforge.net/projects/mrsfast/files/mrsfast-ultra-3.1.0/

BFAST

http://sourceforge.net/apps/mediawiki/bfast/index.php?title=Main_Page

SHRiMP (for AB SOLiD color-space reads)

http://compbio.cs.toronto.edu/shrimp/

RazerA 3

http://www.seqan.de/projects/razers/

Address of the bookmark: http://www.biomedcentral.com/1471-2105/14/184

Source:

http://news.ufl.edu/2014/04/28/real-time-genome-sequencing-at-sea/

miRNA Sequencing data analysis : http://tools.genxpro.net/omiras/

miRNApath( R based tool) :  http://www.bioconductor.org/packages/release/bioc/html/miRNApath.html

miRWalk DB : http://www.umm.uni-heidelberg.de/apps/zmf/mirwalk/

TargetScanHuman : http://www.targetscan.org/

RNAhybrid : http://bibiserv.techfak.uni-bielefeld.de/rnahybrid/welcome.html

RNA22 predictor : http://cbcsrv.watson.ibm.com/rna22.html

miRNA predictor : http://www.microrna.org/microrna/home.do

Plant miRNA DB :http://bioinformatics.cau.edu.cn/PMRD/

miRBASE DB: http://www.mirbase.org/

Plant RNA predictor : http://plantgrn.noble.org/psRNATarget/

miRNA Interaction DB : http://starbase.sysu.edu.cn/

Sequencing based miRNA DB : http://mirgator.kobic.re.kr/

predicted A-to-I edited miRNA DB : http://microrna.osumc.edu/mireditar/

Animal, plant and virus miRNA DB : http://lemur.amu.edu.pl/share/php/mirnest/

Atlantic Salmon miRNAs DB : http://www.molgenv.com/ssa_mirnas_db_home.php

miRNA prediction on UTRs : http://genie.weizmann.ac.il/pubs/mir07/mir07_prediction.html

Idea of analysing miRNA Sequencing data :

http://www.illumina.com/applications/epigenetics/small_rna_analysis.ilmn

More:

http://www.bioconductor.org/help/search/index.html?q=miRNA+target