Questions tagged [rna-seq]

Questions should include this tag if they pertain to issues related to bioinformatics analysis of RNA-seq data, e.g. normalization, differential expression analysis, sequencing and experimental design or transcriptome assembly...

Filter by
Sorted by
Tagged with
22 votes
3 answers
2k views

What is the actual cause of excessive zeroes in single cell RNA-seq data? Is it PCR?

First, sorry if I am missing something basic - I am a programmer recently turned bioinformatician so I still don't know a lot of stuff. This is a cross post with a Biostars question hope that's not ...
Martin Modrák's user avatar
20 votes
3 answers
8k views

How exactly is "effective length" used in FPKM calculated?

According to this famous blog post, the effective transcript length is: $\tilde{l}_i = l_i - \mu$ where $l_i$ is the length of transcript and $\mu$ is the average fragment length. However, typically ...
user172818's user avatar
  • 6,505
19 votes
2 answers
535 views

Confirm success or failure of RNA-Seq normalization

I am working with a set of (bulk) RNA-Seq data collected across multiple runs, run at different times of the year. I have normalized my data using library size / quantile / RUV normalization, and ...
Scott Gigante's user avatar
18 votes
2 answers
3k views

How can we distinguish between true zero and dropout-zero counts in single-cell RNA-seq?

In single-cell RNA-seq data we have an inflated number of 0 (or near-zero) counts due to low mRNA capture rate and other inefficiencies. How can we decide which genes are 0 due to gene dropout (lack ...
Peter's user avatar
  • 2,634
17 votes
4 answers
11k views

How to compute RPKM in R?

I have the following data of fragment counts for each gene in 16 samples: ...
Iakov Davydov's user avatar
16 votes
2 answers
24k views

Difference between CPM and TPM and which one for downstream analysis?

What the difference between TPM and CPM when dealing with RNA seq data? What metrics would you use if you have to perform some down stream analysis other than Differential expression for eg. ...
novicebioinforesearcher's user avatar
15 votes
2 answers
2k views

Alignment based vs reference-free (transcriptome analysis)?

I want to focus on transcriptome analysis. We know it's possible to analyze RNA-Seq experiment based on alignment or k-mers. Possible alignment workflow: Align sequence reads with TopHat2 Quantify ...
SmallChess's user avatar
  • 2,699
13 votes
4 answers
4k views

What methods are available to find a cutoff value for non-expressed genes in RNA-seq?

I have a gene expression count matrix produced from bulk RNA-seq data. I'd like to find genes that were not expressed in a group of samples and were expressed in another group. The problem of course ...
Peter's user avatar
  • 2,634
13 votes
2 answers
4k views

Normalization methods with RNA-Seq ERCC spike in?

ERCC spike-in is a set of synthetic controls developed for RNA-Seq. I'm interested in using it to normalize my RNA-Seq samples. In particular, I'd like to use the spike-ins to remove technical bias ...
SmallChess's user avatar
  • 2,699
12 votes
5 answers
256 views

How do I efficiently perform a metagenome screen of “all” species?

I’ve got an RNA-seq dataset with a large proportion of environmental RNA “contamination”. BLASTing random reads reveals that much of the data comes from bacterial, plant and viral RNA. My target ...
Konrad Rudolph's user avatar
12 votes
2 answers
4k views

determining doublets in single-cell RNA-seq

Doublets are a known problem with scRNA-seq experiments, where 2 or more cells are sometimes captured instead. To determine their presence, there are studies that mix multiple species (such as human ...
burger's user avatar
  • 2,179
11 votes
2 answers
1k views

Can I model technical replicates in DESeq2?

I’d normally use collapseReplicates (or do the collapsing upstream) to handle technical replicates. However, in my current RNA-seq experimental design, samples ...
Konrad Rudolph's user avatar
11 votes
1 answer
1k views

Quantifying reads mapping to multiple loci

I have been using STAR for our RNA-Seq samples. The final.out log file reports percentage of uniquely mapped reads along with percentage of reads that map to ...
rightskewed's user avatar
11 votes
1 answer
311 views

Classifying samples based on marker gene expression

I have a few sets of marker genes that I can classify RNA-seq samples using semi-supervised clustering. I would like to automate the process, however, I am struggling to find the ideal algorithm that ...
GWW's user avatar
  • 752
10 votes
4 answers
961 views

What methods exist to calculate RNA expression profile similarity

Some of the work in our lab requires a comparison of a strain across several experimental conditions. We are looking to identify most similar experimental conditions based on the gene transcription ...
chiffa's user avatar
  • 201
10 votes
2 answers
892 views

Building STAR Genome Index for nanopore RNA sequencing

I am aligning a dataset of 1,000,000 reads oh human mRNA sequenced on Oxford Nanopore Technologies' MinION, and would like to use the STAR aligner, using the parameters recommended by Pacific ...
Scott Gigante's user avatar
10 votes
1 answer
5k views

How to apply upperquartile normalization on RSEM expected counts?

I see that TCGA RNASeq V2 RSEM data is normalized with upper-quartile normalization. After doing Quantification with RSEM with the samples I have, I got "genes.results" as output which has gene id, ...
stack_learner's user avatar
10 votes
2 answers
2k views

*very* unbalanced group sizes for DE

I downloaded some publicly available RNA-seq data and want to compare those samples carrying a mutation (~4) against the rest (~800!). I ran both EdgeR and DESeq2, and the first results in an ...
Kraken's user avatar
  • 405
9 votes
3 answers
3k views

How to identify gene expression signatures from gene expression data?

I have TCGA gene expression data. I'm interested in identifying gene expression signatures using the data. I would like to know whether there are any tools or R packages for identifying gene ...
stack_learner's user avatar
9 votes
2 answers
4k views

Run kallisto iteratively across many samples

I am on a Mac using UNIX. I am trying to use the kallisto quant command on all files in a directory (instead of manually entering them). Because I'm running the analysis against the same index file, I ...
ZincFingers's user avatar
9 votes
3 answers
940 views

Visualisation of long read RNA-Seq splicing

I have a dataset of Oxford Nanopore cDNA reads. Many of my reads are full-length or close to full-length transcripts, and I and am interested in examining alternative splicing. For this, I would like ...
Scott Gigante's user avatar
9 votes
3 answers
4k views

Convert R RNA-seq data object to a Python object

I have done some work in R and would like to try a Python tool. What is a good way to import the data (and its annotations etc) as a Python object? I am particularly interested in converting a ...
Peter's user avatar
  • 2,634
9 votes
2 answers
2k views

Duplicate genes with RSEM counts: Which one to choose?

I have Ensembl ids in the first column and samples with RSEM counts data in other columns. I converted Ensembl ids to gene symbols. Now I see there are three genes repeated twice. ...
stack_learner's user avatar
9 votes
0 answers
145 views

Run cuffcompare in strand-agnostic mode

Is there a way to run Cufflinks' cuffcompare in a strand-agnostic mode? I would like to do this because I have some RNA-seq datasets derived from an unstranded run, that should be compared to a ...
aechchiki's user avatar
  • 2,676
8 votes
4 answers
320 views

How can the cell line contribution be estimated from RNASeq data?

Using a laser-capture microdissection of cells a group of cells stained with the marker of interest was sequenced. In another cohort of patients (this is all human liver tissue) the whole tissue was ...
llrs's user avatar
  • 4,693
8 votes
2 answers
4k views

Which measure should be used in a PCA or RNA-seq data? TPM or counts?

I'm trying to understand the magnitude of batch effects in my RNA-seq samples, and I was wondering which expression units are more suitable to draw a PCA. I'm thinking of either ...
mgalardini's user avatar
8 votes
2 answers
10k views

Are mitochondrial genes to exclude in scRNA-seq such as ribosomal genes?

In this answer, it is stated that ribosomal genes should be excluded prior to normalization in scRNA-seq as contaminants. Do mitochondrial genes have to be excluded as well? I plotted the top 50 ...
gc5's user avatar
  • 1,783
8 votes
4 answers
104 views

Introduce errors in reference transcripts according to external dataset error model

I would like to modify some reference transcripts from Ensembl (D. melanogaster) to introduce a controlled rate of random errors in the sequences. The idea would be to introduce random base ...
aechchiki's user avatar
  • 2,676
8 votes
2 answers
981 views

Getting a "system is computationally singular" error in sleuth

I am analysing 142 samples belonging to 6 batches. Additionally, those samples belong to 72 strains, which means that for most of the strains there are two samples. I could fit simple models (for ...
mgalardini's user avatar
8 votes
1 answer
2k views

Interpreting Intergrative Genomic Viewer (IGV)

I was following a tutorial on "Tuxedo Genome Guided Transcriptome Assembly Workshop" and was wondering how to interpret the following: From what I understand from 'Color Legends', the color blue ...
AlwaysTrying44's user avatar
8 votes
1 answer
418 views

When performing differential expression analysis, should genes with low read counts be removed before or after normalization?

I have RNA seq data which I've quantified using Kallisto. I'd like to use tximport to transform the read count data into input for EdgeR, following the R code supplied in the tximport documentation: ...
J0HN_TIT0R's user avatar
8 votes
1 answer
589 views

lower mapping rates in salmon v0.13 compared to previous versions

Hi there :) Thanks for the tool! I recently updated to the new salmon (from 0.8... its been a couple years) and I noticed that my mapping percentages change dramatically between the two versions. For ...
Courtney Stairs's user avatar
8 votes
1 answer
707 views

Sleuth: transcripts with beta close to 0 are considered differentially expressed in a likelihood-ratio test

I'm comparing the results that I obtain when doing a DE analysis with the Wald test and the likelihood-ratio test. One the thing that I've noticed is that there are many genes with 'beta' close to ...
elsoja's user avatar
  • 241
7 votes
2 answers
559 views

differential gene expression complex design no replicates

We have an experimental design as seen below Where we administered drug at 0 min for each mouse genotype and took them down at given below intervals. wt and ko mouse models were administered only ...
novicebioinforesearcher's user avatar
7 votes
2 answers
2k views

Stranded vs. unstranded library preparation protocols in RNAseq

I've been reading this paper lately: Sailfish Enables Alignment-Free Isoform Quantification from RNA-Seq Reads Using Lightweight Algorithms I don't really understand the second paragraph under ...
Paghillect's user avatar
7 votes
4 answers
10k views

How to convert gff to gtf?

My annotation file is in .gff format. I would like to convert it to .gtf format or to know if there is a way to directly download the annotation file in .gtf format? I am working on sequences from ...
Diango's user avatar
  • 161
7 votes
2 answers
4k views

Filtering step for read counts data

I have around 1200 samples as columns and 60,000 genes with Htseq-Counts data. Before normalization with voom function I want to do filtering step. I want to remove genes whose expression is == 0 in ...
stack_learner's user avatar
7 votes
3 answers
2k views

RIP-seq analysis?

Given an experiment consisting of an input (baseline RNA) and IP (pulldown to find RNAs bound to certain protein of interest)... Is a DE analysis performed over the RNA-seq data from the samples (lets ...
Kraken's user avatar
  • 405
7 votes
2 answers
562 views

Missing genes and normalisation of RSEM output using EBSeq

Without going into too much background, I just joined up with a lab as a bioinformatics intern while I'm completing my masters degree in the field. The lab has data from an RNA-seq they outsourced, ...
J0HN_TIT0R's user avatar
7 votes
2 answers
364 views

Which tools can detect chimeric RNA (fusion genes) from WGS or RNA-Seq data?

Given WGS data or RNA-seq data, which tools can I use to detect gene fusions?
L42's user avatar
  • 281
7 votes
1 answer
1k views

Coverage calculation: long reads (RNA-seq)

Say your aim is to calculate the coverage of an RNA-seq experiment generated with long-read sequencing (so, uneven read length). Up to now, I relied on the Lander/Waterman equation: $$C = L*N / G$$...
aechchiki's user avatar
  • 2,676
7 votes
1 answer
158 views

How can HISAT2/StringTie report decimal coverage values

I have performed RNA-seq analysis using HISAT2 & StringTie workflow suggested in: Transcript-level expression analysis of RNA-seq experiments with HISAT, StringTie and Ballgown. Some of the ...
pogibas's user avatar
  • 203
7 votes
2 answers
2k views

Spliced vs. unspliced ratios for transcripts in RNA-seq data

Is there a computational tool for measuring what percentage of RNA is spliced in an RNAseq experiment? I'm not particularly interested in complicated analyses that give ratios for all possible ...
Jessime Kirk's user avatar
7 votes
1 answer
169 views

Network comparison of single cells (from sequencing data)

What methods can we use to compare networks (PPI, gene regulatory etc) created from single-cell gene expression (or proteomics) data? There are methods that construct networks by comparing two ...
Peter's user avatar
  • 2,634
7 votes
1 answer
672 views

insert size pre and post trimming

I have a problem here with my rna seq data: Sequencing details: rRNA was removed, followed by cDNA preparation and generation of stranded libraries using the TruSeq Stranded Total RNA Sample Prep Kit. ...
novicebioinforesearcher's user avatar
7 votes
2 answers
597 views

Why are TPMs per 10k or 100k in many scRNA-seq studies?

I noticed that many scRNA-seq papers normalize TPMs to 10k or 100k as opposed to 1M (as the abbreviation defines them). It doesn't really matter since you are just moving the decimal point, so why ...
burger's user avatar
  • 2,179
7 votes
2 answers
489 views

Correct for gene length or read counts in GO enrichment analysis

It is a well reported fact that GO analysis of RNAseq results is affected by a number of biases, including length bias and expression level bias. The bioconductor ...
Ian Sudbery's user avatar
  • 3,311
7 votes
1 answer
728 views

Why are my kallisto and salmon results differing so much just for lncRNA transcripts?

I am running some analysis on an RNA-seq dataset. I have a list of transcripts that are potential lncRNA for which I ran both Kallisto and Salmon aligners. The input data for index building and ...
Marta Silva's user avatar
6 votes
3 answers
8k views

Volcano plot in R

This question has also been asked on biostars How can I reproduce this volcano plot? I'm only able to do the traditional one, I'm kind knew too these field.
Sofia's user avatar
  • 351
6 votes
2 answers
294 views

Clarification on Gene Enrichment

When I run a GSEA analysis on two conditions from the same RNaseq (negative control PBS injection VS positive control CpG injection) from the same dataset/same gene list, I get results that look ...
julianstanley's user avatar

1
2 3 4 5
14