Questions tagged [deseq2]
Deseq2 is an R package for analyzing RNA sequencing data
138
questions
0
votes
0
answers
16
views
How to use compareCluster in ClusterProfiler for Multiple GSEA Comparisons in a single GSEA analysis?
I am new to using compareCluster and just want to double check - the input for compareCluster can be ...
- 13
1
vote
1
answer
22
views
Normalizing AMR Gene Read Counts in Metagenomics Data Without Replicates: Is RPKM the Best Approach?
Following up from my previous question, I am currently conducting antimicrobial resistance (AMR) surveillance in water samples using long-read metagenomics (Nanopore sequencing). After analyzing the ...
- 157
0
votes
1
answer
43
views
Deseq2 model formula and longitudinal experimatal (or time-series) designs
At hand I have RNA expression data of 6 organs from mice with treated and untreated diets, collected at consecutive time points, and I wonder how to design an analysis accordingly in ...
- 3
1
vote
1
answer
21
views
Differential gene expression (aiming to create heatmap) for raw RNAseq count data from different publications - how to?
I'm new to bioinformatics and RNAseq analysis, I'll try my best to explain the question!
I have raw counts RNAseq data from 3 different publications (so 3 different datasets) that all had similar ...
1
vote
1
answer
33
views
Normalization and Comparison of AMR Gene Tables from Long-Read Metagenomics Data without Spike-In Controls
I am conducting antimicrobial resistance (AMR) surveillance in water samples using long-read metagenomics (Nanopore sequencing). After analyzing the data with KMA (a k-mer read alignment tool, Clausen ...
- 157
1
vote
1
answer
44
views
PCA of bulk RNA-seq doesn't show clustering and show large variation in general
I am very new to analyzing RNAseq and I am in a group with very little experience in this regard and I am looking for some advice.
My PCA after performing DESEQ2 analysis on my dataset doesn't show ...
- 11
2
votes
2
answers
69
views
Reasons for extremely low number of DESeq identified differentially expressed genes after RNAseq?
I am fairly new at RNA-sequencing analysis, but I have attempted to analyse the data I obtained from RNA sequencing on my own by following an online EdX class and by basing the majority of the ...
- 23
1
vote
0
answers
26
views
Is it normal to have low mapping to bacterial genome in a total RNA sample of plant root-colonized bacterial cells?
I had isolated the total RNA from a sample of plant roots with colonized bacteria by excising the roots (harboring the bacterial cells) and crushing them into a powder using liquid nitrogen. The RNA ...
- 149
2
votes
1
answer
37
views
rna seq matrix design hybrids
I am struggling to get my head around the correct design to use for my experiment.
I have RNAseq of 2 pure species (s1, s2) and of the hybrid (h). I aligned all three genotypes to the s1 genomes. I am ...
- 165
1
vote
1
answer
54
views
Multi Factor in Deseq2 Gene enrichment analysis
I want to see how the gene expression differs in breast cancer between three species, and I am using DESeqDataSetFromMatrix on my count table.
...
- 69
1
vote
0
answers
68
views
Continuous predictor variable to identify responder genes using DESeq2
I have samples from 45 individuals at 2 time points and these individuals were given a specific diet. The BMI of these individuals were recorded at both time points. The responders to the diet are the ...
- 237
2
votes
0
answers
65
views
How to incorporate negative controls in DESeq2
We are doing a comparison between two outcomes (positive and negative). We could not have any positive controls as we do not have any "control" data to set as baseline, either from ...
- 353
1
vote
0
answers
224
views
how is log2foldchange calculated in DESeq2
I am trying to calculate log2FoldChange in DESeq2 by hand because I am trying to reproduce my results and understand how DESeq2 works.
What I want to know is:
Are the counts first fitted to a ...
- 11
1
vote
0
answers
39
views
enrichment of DE genes within wide regions not pathways, hypergeometric test or GSEA?
I haven't found any relevant posts, so here is my question (it's more a question on statistics I believe):
I have the deseq2 output of DE genes( around 25000 genes) , and among them 1000 are ...
- 11
2
votes
2
answers
182
views
Understanding DESeq2 output for resultNames
I'm following along with some pre written code to learn how to do differential expression. I am using the same data set and have gone back through the code and copy and pasted what they used and the ...
- 73
3
votes
1
answer
205
views
Issue with Differential Gene expression analysis with Deseq2
I am using the Deseq2 package to perform a DE analysis on multiple factors.
So, basically to perform the analysis I consider patients characterized by the Brugada Syndrome type 1 and other patients ...
- 31
2
votes
1
answer
275
views
Correlation heatmap of RNA-seq clusters all samples together leading to very low no. of DEGs
I am writing to you to take an input or may be you can provide a different perspective. I am at wits end :(
So I have nearly 200 samples. These are separated into two groups (group I treated with ...
- 237
2
votes
2
answers
340
views
Integrate bulk RNA and ATAC-seq genes
I have up and down-regulated genes from bulk RNA result from DeSeq2. Genes that are differential in chromatin accessibility between two conditions from DiffBind. I would like to find the genes overlap ...
- 101
1
vote
1
answer
204
views
DESeq2 throwing error while normalizing raw microarray expression data due to presence of negative values
This question was also asked on Biostars
I am trying to download and analyze a miRNA expression dataset from NCBI GEO (GSE25631). I specifically want non-normalized ...
- 145
2
votes
2
answers
64
views
RNAseq alignment: best practices for aligning to multiple isoforms?
I have Illumina RNAseq data and would like to maximize my power to find candidate genes that are differentially expressed genes between experimental conditions.
Many of my (de novo assembled and ...
- 23
1
vote
1
answer
724
views
Differential gene expression: how do I correctly compare three groups to each other?
I'm trying to figure out the right way to do differential gene expression analysis with a discrete variable with three groups. The context is, I want to assess differential expression as a function of ...
- 11
0
votes
1
answer
142
views
Batch correction in differential expression analysis
I have sent two sets (two batches in matter of sending for sequencing) of different samples (plasma) to small RNA-seq to Qiagen company
This is how my meta data look
...
- 2,232
1
vote
0
answers
1k
views
Manually calculating log2 fold change values from DESeq2 normalized counts [closed]
This question was also asked on Biostars and Bioconductor Support
My aim is to perform gene clustering based on RNA seq datasets (raw expression values) from NCBI database. After clustering I realized ...
- 51
3
votes
2
answers
619
views
DESeq2 for large number of samples takes too much RAM
I am trying to run a very large number of transposase-accessible chromatin (ATAC)-seq samples through DESeq2 to find peaks of differential chromatin accessible across my study genome.
I have about 100 ...
1
vote
1
answer
101
views
DESeq2 and EdgeR
I am new to using both DESeq2 and EdgeR in Bioconductor used for transforming my RNA expression data.
However, I am struggling to understand their specific purpose, differences between them and ...
- 13
1
vote
1
answer
80
views
Filtering reads greater than 5 from HT-seq count files
I have some raw counts from HTSeq after aligning with hg38 human reference genome. I want to do filtering in a way that the filtered count files should have the same number of lines. The reason behind ...
- 143
1
vote
0
answers
42
views
Creating multiple phenotype datasets using bootstrap method "Bootstrap-samples-by-column-of-a-data-frame-in-r" for DEG analysis
I am working on a datasets and after some discussion with my group, we doubt that maybe one or more of our controls are different than the other controls. The motivation is to see if one or more ...
- 31
2
votes
4
answers
342
views
nested samples in design DEseq2
I'm running a DEG analysis with Deseq2 by specifing the following design that includes 12 samples from 6 pts in 4 different conditions:
...
- 21
1
vote
2
answers
107
views
Basic RNA Differential Expression in R
I have two matrices, one for individuals before treatment and one for the same individuals after treatment. Both matrices are raw read counts of RNA expression.
...
1
vote
1
answer
79
views
deseq2 single factor design output interpretation
This is regarding the single factor design For example if I have Age or other continuous numerical variable how to provide that into the design formula.
For this post do i need to 'You could ...
- 1,854
1
vote
1
answer
123
views
3
votes
1
answer
291
views
deseq2 makeContrasts function
This article talked about various deseq2 designs etc. One of the designs I would like to use is explained as this:
Control versus treatment average
...
- 1,854
1
vote
1
answer
909
views
deseq2 full and reduced model interpretation
I would like to know how to interpret the output of the formula although i been thorough loads of literature but I'm not confident yet.
So this is my full model I'm running
...
- 1,854
0
votes
0
answers
79
views
Different results of DESeq2 in different ways of comparison
I have different treatments (infected and control ones). When using DESeq2, I have different results in two ways of comparison.
first way:
factor: treatment, factor level1: Moribund Infected, factor ...
- 11
2
votes
1
answer
713
views
What is the difference between Normalized Expression in EdgeR vs DESeq2?
I am trying to access the normalized expression in both edgeR and DESeq2, yet the results are different. Does anyone know why?
How to get normalized expression using edgeR:
...
- 165
2
votes
1
answer
381
views
baseMean threshold
I have an RNA-seq dataset and I am using DEseq2 to find differentially expressed genes between the two groups. I used pre-filtering to remove any genes that have ...
- 21
1
vote
1
answer
163
views
Error in DESeq2 analysis with dimnames
I'm getting the following error in a DESeq2 analysis and I can't seem to figure out the issue:
...
- 165
3
votes
0
answers
80
views
LRT test in deseq2 on multiple condition
I have some doubts when i'm running LRT
This is my condition data-frame small subset
...
- 1,854
1
vote
2
answers
664
views
How to input data and metadata from NCBI for RNA-Seq analysis in R
So, I have downloaded some .txt and .csv files containing the raw count data (I am supposed to find out differentially expressed genes and see which genes and up- and down-regulated). I have also ...
- 11
0
votes
1
answer
769
views
How does DESeq2 "collapseReplicates()" function work on read counts data?
Comparing read counts from an RNA-seq experiment for two select genes before and after using DESeq2's collapseReplicates() and ...
- 335
1
vote
0
answers
60
views
Batch Effects that are Close to Collinear with Main Effect
Two years ago, I started working on a project uses both RNAseq and ATACseq. It's supposed to be a simple Healthy Control (HC) vs disease study. The sample collection was done in 2 phases, and it was ...
- 177
0
votes
1
answer
305
views
How to remedy a DESeq2 collapsing technical replicates error?
Goal: To ensure "the sum of the counts for [my samples] is the same as the counts in the [samples] columns in ddsColl" after collapsing replicates using ...
- 335
0
votes
1
answer
2k
views
Use of "baseMean" in DeSEQ2
The basemean is described as the "mean of normalized counts of all samples". My question is, how would we interpret differences in basemean between genes, pertaining to reliability and such? ...
-1
votes
1
answer
90
views
Visualizing top expressed genes
I have a spreadsheet of CPM values
I want to visualize the top 20 genes expressed across my samples
This is how my data look like
...
- 2,232
1
vote
1
answer
144
views
Best R program to use for RNA seq analysis
We have an RNA seq taken from a the same tissue in two different developmental stages. We have one gene, with the number of reads from stage A and stage B, and we want to see which genes are ...
- 11
1
vote
0
answers
51
views
The correct design for time series experiment
I have a treatment and control in two time points like this
...
- 2,232
1
vote
0
answers
40
views
How to simulate replicates for DGE analysis?
I am prototyping with data visualization of DGE results, and I would like to work on the analysis pipeline before the real data is available. Currently, I only have 3 samples for wild type and 1 ...
0
votes
2
answers
383
views
Help with performing PCA analysis on data from 3 different datasets of normalized counts data for RNA-seq experiment
So, I have limited knowledge of R but I need to do a PCA analysis of 3 different datasets of gene expression as a result of combined growth or mono-culture growth.
The 3 different datasets I performed ...
- 41
0
votes
1
answer
121
views
Understanding the different designs in DESeq2
I am using DEseq2 and trying to understand the results obtained using different models.
I have a data design with 2 genotypes and 2 time points.
...
1
vote
1
answer
233
views