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
1
vote
1answer
21 views

Error to running "CIRIquant" function in CIRIquant package at python Language

"CIRIquant" is a Python package for detecting Circular RNA. I installed "CIRIquant" through the below command in terminal: pip install CIRIquant My Linux distribution is Ubuntu ...
1
vote
1answer
45 views

Different contrasts in DESeq2

I have a treatment and controls in two time points like this ...
0
votes
1answer
29 views

Convert VCF file to mpileup.txt

I am working on an iterative analysis that uses orthologous pipelines that require mpileup.txt files as input for a visualization step. This requires me to convert VCF files to mpileup.txt. This ...
2
votes
0answers
12 views

Alternatives to CAP-miRSeq for microRNA processing

I am involved in a Drosophila transcriptomics project in which we have been obtaining RNAseq data for microRNAs as well as mRNAs. In initial work we used CAP-miRSeq for the latter as the Tuxedo ...
0
votes
0answers
39 views

Why is bcl2fastq2 taking so long to calculate stats?

Our lab has been using bcl2fastq v2.20.0.422 to demultiplex RNA-seq data sequenced on an Illumina Novaseq machine on a beefy EC2 instance and we've run into the strange problem: namely that while ...
0
votes
0answers
33 views

How to perform a meta-analysis using data consisting of paired-end and single-end reads generated from Illumina and Ion Torrent?

So basically I have RNA-seq reads that were generated from Illumina and Ion Torrent platforms for yeast species. I have seen an article where they compared liver cells of a rat that were sequenced ...
1
vote
2answers
52 views

Is it possible to do DEG analysis without replicates?

I have three bulk RNA-seq samples (test1, test2, test3) without replicates. And I noticed that DEG analysis tools such as DESeq2/edgeR/etc cannot be applied for data with no replicates. So, I just ...
1
vote
1answer
35 views

Assembling all transcripts for an individual gene? (using single sequence to seed the assembly)

Let's say I have a candidate gene and I believe that in an individual sample, the genome sequence differs from the reference which then interferes with alignment. Is there a way for me to do a "...
0
votes
0answers
12 views

java error when running CIRI-full.jar Merge module

I have an issue running a module trying to detect circRNA from RNAseq data. The error is a java error and I have not knowledge of java so I don't understand the issues and googling hasn't helped. the ...
0
votes
0answers
18 views

A very large number of clones in BCR reportorie

I am using mixcr to convert fq.gz raw data file (single cell BCR sequencing) to txt files with the names such as JX01_d3-B.clonotypes.IGH.txt , Then I use the immunarch to load the file to explore the ...
1
vote
2answers
44 views

Counting the number of gene isoforms in a GFF3, is this method correct?

Recently I've been tasked to count the numbers of gene isoforms for each locus in a .gff3 file. I'm still doing my first steps in biology and bioinformatics, so I ...
0
votes
1answer
32 views

Error when using awk command to trim sequence files

I am attempting to edit my fastq sequence files with an awk command from a published article. The sequence files are all ...
0
votes
0answers
50 views

Splitting .bam files into separate samples for tophat2

I have RNA seq data of 12 samples (paired end, so 24 files) and I ran tophat2 for alignment on all of them which gave me one bam file for all the samples. Now I need separate bam files for all the ...
0
votes
1answer
36 views

Limma decideTests function: what kind of multiple hypothesis testing correction does parameter "method" involve?

What kind of multiple hypothesis testing correction does method="global" do in Limma's decideTests function? According ...
0
votes
0answers
23 views

Gene co-expression network analysis

I am working on a genome-wide study and struggling with gene co-expression network analysis. I am working on Tomato(Solanum lycopersicum). Can someone suggest me some tools or methods for this ...
0
votes
1answer
28 views

Keeping rows with a cut-off in all columns

I have a data frame of raw read counts, genes in rows and samples in columns ...
1
vote
1answer
45 views

shiny dashboard to visualise RNA seq data

I am trying to build a shiny dashboard for the visualization of different RNAseq data sets. However, I am encountering problem with reactive input, the code is always using the second dataset that I ...
0
votes
1answer
58 views

Beginner with DESeq2 having issue with analysis

Hi there I am brand new to using RStudio and trying to run DESeq2 analysis on my featureCounts output counts table. I ran the following code: ...
0
votes
1answer
30 views

Parallelizing microRNA targets

I am trying to look for miRNA targets using a file called Conservedfamily.txt from the Zebrafish target scan fish website. I have written a python program to ...
2
votes
1answer
47 views

DESeq2. Which is better: using the altHypothesis argument in the results function or manually filtering for your desired P-value and Log2FoldChange?

In DESeq2 you can identify your differentially expressed (DE) genes using the results function. I noticed people in my lab using the results() function with the minimum number of arguments supplied (...
0
votes
3answers
91 views

The confusion of using TPM (transcripts per million)

It is shown that TPM values are not suitable for DEG analysis but good for within-sample comparison since TPM normalized the gene length. My question is first: if TPM is not suitable for across ...
1
vote
0answers
9 views

Taiji doesn't run Pagerank

I'm trying to use Taiji to run a pagerank on our combined bulk RNAseq and ATACseq dataset. Here's how we run it: taiji run --config config.yml -n 10 Here's ...
1
vote
1answer
49 views

how to get normcounts for singlecellexperiment object?

I need to perform differential expression analysis using the scDD package from R, but I am not able to since I miss the normcounts assay in my SCE object (of course in the example they show, the assay ...
1
vote
1answer
48 views

Understanding "Measurement of mRNA abundance using RNA-seq data" by Wagner et al 2012

I'm trying to understand why RPKM is not an appropriate way to normalize for RNA seq (I understand the general idea, but I'd like to gain a deeper understanding). So I'm reading the original paper by ...
0
votes
0answers
40 views

Correcting for batch effect in bulk RNA seq datasets

I have samples from many different bulk RNA seq studies that exhibit an evident batch effect (samples from different labs cluster together regardless of treatment). I've found a lot of papers (here, ...
0
votes
1answer
17 views

Combining read counts from three separate GEO studies

I want to do differential expression analysis with DESEQ2. I have three read counts files downloaded from GEO (small RNAseq based) where the number of miRNAs and id is nearly the same. These studies ...
2
votes
1answer
60 views

Batch effect correction using a subset of samples - using DE genes

I have some RNA-seq data with two very obvious batches as you can see in the PCA: The samples of interest (A - H) are from tumor tissues. In addition, there are data for two cell lines I and J. The ...
0
votes
0answers
15 views

Installation of software Capmirseq

I am trying to install a software called Capmirseq. After installing all relevant modules I followed the installation instruction code as follows ./install.sh -p /path/to/software. Which is expected ...
0
votes
1answer
43 views

RNA_Seq data aligned used uniquely or multi mapped reads impact on result interpretation

I have some transcriptomic (Whole) sequencing data that I should analyse. I would like to do raw data alignment to a reference genome taking into account the multi mapped reads and uniquely mapped ...
1
vote
0answers
12 views

Ccle miRna data units

Does anybody know what units the mirna expression data on ccle is in? I saw the pipeline in the original paper but it doesn't mention units. Original paper: https://www.nature.com/articles/s41586-019-...
1
vote
0answers
32 views

How to create a DESeqDataSet and define experiment design before varianceStabilizingTransformation?

I have an RNAseq count matrix consisting of 2 groups (high, low) with 6 timepoints per group (T1,T2,...,T6) and 3 replicates per timepoint (rep1, rep2, rep3). So a 2-factor design with 36 samples in ...
2
votes
1answer
48 views

How does FDRtool work?

I have a question about using FDRtool. In the below code (on RNA seq data whose p values were acquired using Deseq2), the FDRtool was first used and thereafter p.adjust using the benjamini hochberg ...
3
votes
2answers
163 views

Calculating most abundant transcript from RNA-Seq data

vcf2maf uses VEP to annotate variants, and I believe selects the default Ensembl transcript to use for annotation. Sometimes the ...
2
votes
1answer
48 views

Get Salmon mapping/alignment summary

With HISAT2, after the alignment of fastq files you get an alignment summary like this: ...
0
votes
1answer
46 views

Base abbreviations other than ACTGU in sequence file

Short question: where can I find what F and J and < and ...
1
vote
1answer
73 views

What is the meaning of split read?

I want to use rna seq data to later perform functional tests on fusion genes. so before that I need to filter the "best results" (of rnaseq) for deciding which candidates I actually want to ...
1
vote
1answer
153 views

What are common ways to calculate gene length for TPM calculations? [duplicate]

I am aware of this similar question. But the accepted answer there answers how to calculate TPM given a mapping from gene name to gene length. My question is, given an annotation file of genes (for ...
0
votes
2answers
42 views

RNAseq DE comparison with samples of different read length [duplicate]

I have RNA samples but with different read lengths (Eg, HiSeq 2x125 and NovaSeq 2x150bp data). I would like to do DE analysis on these samples. What do you recommend? Do you recommend to trim the ...
1
vote
5answers
137 views

Normalize RNA seq data from multiple runs for expression analysis

I have RNA samples sequenced with TruSeq Stranded Total RNA kit protocol in Illumina HiSeq (2x125bp) and NovaSeq platforms (2x150bp) - almost 100 samples altogether. I have to use the samples data for ...
1
vote
1answer
47 views

correlation coefficient versus DEGs analysis: what's the best approach for low expressed genes?

I have 6 folders. Each one contains 7 datasets of a specific type of cancer (RNA-Seq) and 7 datasets of normal tissue (healthy control for that type of cancer). A total of 84 datasets. I want to ...
1
vote
0answers
26 views

Data format for pathway based clustering of samples

I came across this paper as one of the examples from this paper, this one Figure 2. Host Protein Alterations in Infected iAT2s ...
1
vote
1answer
57 views

Differential gene expression analysis of time series with replicates

I have a dataset that has two groups, perturbed vs control. Each group has 3 replicates. Each replicate has 8 time points. How do we do Differential gene expression analysis to find significantly ...
3
votes
1answer
241 views

Why did expression based subtypng of breast cancer gain much more acceptance than others

This is may not be entirely technical question but rather a academic question. But the technique behind the application is within the scope of bioinformatics. So I would try to ask here that: In each ...
2
votes
1answer
46 views

Toptable error, wont recognize condition

I am getting a rather strange error from topTable. When I run my code I get Error in fit$coefficients[, coef] : subscript out of bounds from topTable as if it is ...
1
vote
1answer
45 views

Challenging benchmarks for supervised learning on sparse scRNA-seq data

One challenging aspect of modeling scRNA-seq data is data sparsity, that is, scRNA-seq measurements typically suffer from large fractions of observed zeros (i.e. dropouts), where a given gene in a ...
0
votes
1answer
55 views

Demultiplexing FASTQ file without index information

I am trying to understand how data that was uploaded to SRA (https://www.ncbi.nlm.nih.gov/sra?LinkName=biosample_sra&from_uid=4510743) can be analyzed with the assumption that the FASTQ file ...
0
votes
1answer
81 views

RNASeq analysis using featureCount and EdgeR

I am using a pipeline (bam -> featurecount-> EdgeR) to do some RNASeq analysis of several groups and sub-groups. For example, I have the following dataset with two types (T1 and T2) and T1 has ...
1
vote
2answers
95 views

Un-normalize DESeq2 counts

I obtained a matrix of RNA-seq count data that has been normalized by DESeq2's median of ratio method. I know that DESeq2 wants to take in un-normalized counts, but I do not have access to those data. ...
2
votes
2answers
91 views

Calculate genome coverage and depth from alignment

I have a .bam alignment file and a genome reference .fasta file. I am looking for a easy to use tool (that I can reference in a publication) to calculate the percentage coverage of the reference by ...
0
votes
1answer
45 views

Common gene naming conventions and how to convert between them

I have two questions, one is a direct question about some RNA expression data. The other is more broad, but motivated by the first. I have downloaded an expression set for a class project and I am ...

1
2 3 4 5
11