Questions tagged [differential-expression]
Use it when comparing a measure of the expression between two (or more) samples, where the units are relative to the comparison.
105
questions
0
votes
0
answers
9
views
Unclear variation in level 5 LINCS/L1000 gene perturbation Tx data and best practice to process
This question was also asked on Biostars
I am trying to use level 5 gene KO/KD/OE & compound perturbation data (moderated z-scores for the 978 landmark genes) from LINCS L1000 (source: clue.io) to ...
- 13
1
vote
0
answers
12
views
What is the difference between AnnotationDbi/org.Mm.eg.db and biomaRt/Mus.musculus for converting to gene symbols?
I am interested in the differences between AnnotationDbi/org.Mm.eg.db and biomaRt/Mus.musculus. They yield the same results as seen in the code below.
...
- 153
2
votes
1
answer
51
views
scRNA-Seq pseudobulk for marker detection
A standard approach for scRNA-Seq is to partition the single cells into individual clusters, then use a Wilcoxon test to find markers that characterize each cluster (or other statistical methods that ...
- 202
1
vote
1
answer
46
views
RNAseq small sample size in control group
I am doing an RNAseq experiment with affected (n=6) and control groups (n=6) in cow. However it turned out that 4 of my control samples have very low quality so I have to discard them.
My problem is ...
- 45
2
votes
0
answers
64
views
How can I separate upregulated genes from downregulated on a boxplot?
I have a table of fold change shifts for two different groups. Some of them are upregulated, others are downregulated. I wonder if it is ok to subset only positive values and to boxplot them to show ...
- 103
2
votes
1
answer
30
views
How to compare differential expression of same genes between two groups? [duplicate]
I need a statistically significant "decision" that DEGs from group A are different (or not) from group B (after drug treatment).
For that, I have a mean fold change delta, calculated as the ...
- 103
1
vote
2
answers
64
views
What is the best way to compare differences in DEG expression between two populations?
I am interested in differential response to a drug in the "RED" group versus the "GRN" group. I mean, I need to compare differences in DEG expression between two populations (RED ...
- 103
0
votes
1
answer
58
views
What is the best way to compare treated/control samples? [closed]
I have two populations of mice, say RED and GREEN; there are 5 pairs treatment/control in each.
The table is the following:
...
- 103
0
votes
0
answers
13
views
Why is effective length of transcript better than normal length in Diferential expression gene analysis? [duplicate]
I want to ask, what is a benefit of effective length in RSEM analysis?
Thank you.
0
votes
1
answer
43
views
Comparing differential expression across samples - is batch effect correction needed?
I have a bulk RNA-seq dataset made up of control and treatment conditions for a range of cell lines. This dataset was generated in two batches, such that the cell lines are split between batches but ...
1
vote
0
answers
32
views
How can I get a column with real gene names in my ballgown analysis?
I am doing RNAseq analysis and I am using ballgown procedure for it.
There is an option in R to calculate differentially expressed genes and use FPKM in calculating differential gene expression, which ...
- 103
-1
votes
1
answer
47
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
...
- 1,767
0
votes
1
answer
49
views
Gene expression changes in one condition
I want to look at an RNA-seq data-set on 1 condition to look for expressed genes. As there is only one condition, I can’t perform differential analysis.
The only ...
- 1,767
1
vote
0
answers
37
views
The correct design for time series experiment
I have a treatment and control in two time points like this
...
- 1,767
0
votes
1
answer
41
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
128
views
2
votes
1
answer
50
views
Help with Limma-model
I'm starting to use Limma to find differential expressed sites within my data and would like to ask if my approach was correct.
The data has a control and a treatment set, each containing data of 48 ...
- 23
0
votes
0
answers
39
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 ...
- 41
0
votes
1
answer
45
views
Different results in differential expression analysis of microarray data
I am performing differential gene expression analysis to microarray data for type 2 diabetes donors and nondiabetic donors. When I run the code I get some different results in each time (about 50 or ...
- 1
0
votes
0
answers
62
views
Differential miRNA expression using RPM
I have microRNA (miRNA) expression data in RPM. I would like to do differential gene expression between two groups.
How can I do this?
I have considered edgeR and DESeq2 in R, but it looks like they ...
- 205
3
votes
1
answer
853
views
R limma alternatives in Python
The R package limma is ideal to perform differential expression analysis. Is there any limma alternative in Python?
I'm trying ...
- 51
0
votes
1
answer
137
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:
...
- 41
0
votes
1
answer
316
views
1
vote
1
answer
54
views
Help with DIA-Mass Spectrometry data analysis with several conditions (limma?)
I am trying to learn how to analyse normalised DIA-MS data and I am struggling with it ://
The original dataset I got is (6 conditions (2 samples each)) with 3 technical replicates (total: 36 sample ...
- 11
1
vote
1
answer
45
views
The biological meaning of the random variables and the responses in Seurat analysis
In the Seurat analysis, if we suppose that Xg and Xr denote the random variables that associate to the expression level of the gene g and of the gene r, respectively. Let Y and Yv represent the ...
- 153
1
vote
5
answers
258
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 ...
- 11
1
vote
1
answer
100
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 ...
1
vote
1
answer
38
views
Are there any databases for gene co-expression or expression pattern clustering?
I am currently working on gene clustering based on co-expression pattern in mouse brain. The problem is I do not have some solid way to test my result. Are there any suggestions for databases ...
- 13
0
votes
1
answer
116
views
Clustering RNAseq fold-change data
I have a dataset that looks something like this
Treatment1
Treatment 2
Control
Sample1
3.23.
0.87.
1
Sample2
1.71.
1.
1
Sample3
2.88.
5.65
1
Sample4
0.77.
1.34
1
The numbers describe a fold ...
1
vote
0
answers
32
views
How to call this : some dots in volcano plot aligned a line after add a very small value to padj
when I plotting the volcano plot, I got some dots very top of the Y-axis because their padj nearly reached 0（I checked the data, I think they are reliable）, so I followed the advice to add a very ...
- 13
0
votes
0
answers
56
views
Adding matrix as a confounder in EnrichmentBrowser
Brain data, particualrly human brain data have very heterogeneus cell types, so it's important to normalise by cell type/add proportoins of different cell types to the formula when performing ...
- 85
1
vote
1
answer
298
views
annotation using ChIPseeker package error
I have differential binding sites object obtain from diffBind (dba.report). I am using the ChIP Seeker package to annotate the peaks but keep getting the following error: Error in (function (classes, ...
- 95
0
votes
4
answers
96
views
Where can I find Single Cell Data with Location "Coordinates"?
Does single cell data typically have the following meta-data: the "coordinates" (e.g. on a tissue, adjacent tissues) saying where each cell in the sample was located relative to other cells? ...
0
votes
1
answer
89
views
Is there an established method for comparing expression of groups of genes (gene sets)?
I have a figure in a paper where I show the log2(fold-change) (obtained with DESeq2) between two groups (based on genotype) for genes within a specific hallmark gene set, and a reviewer is asking ...
- 367
1
vote
0
answers
19
views
Questions about EBSeq (in R)
I have made a UMAP of malignant cells and the result is split into 3 clusters, as seen below.
For the sake of example, let's say I want to find the differentially expressed genes in the uppermost ...
- 11
1
vote
0
answers
153
views
Differential Gene Expression with Replicates for some of the samples
[this question has also been posted on Biostars; some additional clarification from there has been copied into this question]
I've been asked to analyse a set of samples in which their control sample ...
- 169
0
votes
0
answers
38
views
a proper Design Matrix for several drug treatments with both control negative and control positive
I have a dataset of RNA-seq samples for testing different drugs on the presence of another drug. One of my samples is the normal cells with no drugs (control negative) and another is the cells with ...
- 169
2
votes
2
answers
180
views
Differentially expressed genes analysis in Seurat
For the differentially expressed genes analysis, is it possible to check for DEGs based on the levels already identified in the object? For example, my dataset contains cells from 11 subjects, which ...
0
votes
0
answers
44
views
Alternatives to imputation for single-cell RNA-Seq analysis to increase statistical power of DE analysis
I have a small scRNA-Seq dataset (n = 357, inhibitory neurons). This set of cells is split almost evenly between two conditions (Case and Control). I would like to test for differential expression ...
- 1,723
0
votes
2
answers
37
views
Separating peaks of chip-seq with specific length
My data contain several chip-seq results. I have the peaks called by MACS2.I wanted to only look at those peaks that their size is e.g 500bp to 1000bp. How can I separate those peaks efficiently?
I ...
- 95
0
votes
1
answer
35
views
GSEA enrichr with 10x genomics differential_expression ranks
I am attempting to use GSEA enrichr with 10x genomics differential_expression rankings. Reading what people seem to be doing with GSEA, there seems to be a pre-/post-singlecell gene expression ...
- 2,234
1
vote
2
answers
386
views
How to create a list of differentially expressed (DE) genes after normalization with RUVSeq?
I am using edgeR to perform differential expression (DE) analysis on a set of RNA-seq data samples (2 controls; 8 treatments). To correct for batch effects, I am using RUVSeq.
I am able to get a list ...
- 115
1
vote
3
answers
2k
views
Is there an "official" standard usage of Red and Green on differential gene expression heatmaps?
I have recently found myself making multiple heatmaps for visualization of differential gene expression results between replicates in two experimental conditions.
The default settings in the plotting ...
- 36
2
votes
1
answer
176
views
Comparing multiple treatments to multiple other treatments in edgeR for simple effects in a complex experimental design
I am working with a RNA-seq data set in maize that has a relatively complex design. There are two levels of treatment A (nitrogen fertilizer level in the field, high or low), two levels of treatment B ...
- 21
1
vote
1
answer
517
views
What unit (TPM, FPKM/RPKM, or other) to use when working across samples
I have raw gene read counts and would like to perform an analysis across multiple samples. I've found conflicting info online on how this should be done. One commonality however is that FPKM/RPKM ...
1
vote
1
answer
54
views
Is it conflicting: significant difference according to Wald test but not in PERMANOVA and ANOSIM to compare species abundances of metagenomic samples?
If there is no significant difference between two groups of metagenomic samples according to "Multivariate differential abundance tests" for example PERMANOVA and ANOSIM, does it makes sense to use "...
- 53
0
votes
0
answers
31
views
Do I need to study alternative splicing and isoformswitch separately?
I have done a study on isoform switch between different tissue type. Now, many people asking me to study alternative splicing also. But, dont they will give me same results? If there is a significant ...
- 171
2
votes
1
answer
104
views
DESeq2 compare within a condition
This might be a really stupid question, but I can't figure out how to do this. I've read through the DESeq2 vignette and manual pages but couldn't find an answer.
I have a bunch of samples split up ...
- 437
2
votes
1
answer
947
views
How to downsample some of the samples in RNA-seq data?
I have 40 samples and these are into two groups. I would like to perform a differential analysis between two groups. The library size of the samples is very low. But there are two samples (GroupA_12 ...
- 1,242
1
vote
1
answer
113
views
Differential Expression
I have genes differentially expressed between two groups (case and control). I would like to annotate them by classifying them according to their biological functions.
I work on parasites of ...
- 121