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.

Filter by
Sorted by
Tagged with
1
vote
1answer
25 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
1answer
52 views

Different contrasts in DESeq2

I have a treatment and control in two time points like this ...
2
votes
1answer
30 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 ...
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 ...
0
votes
0answers
23 views

Different results in Differential Expression analysis of Microarray Data

I'm new to R, so be patient with me. 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 ...
0
votes
0answers
18 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 ...
4
votes
1answer
219 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 ...
0
votes
1answer
59 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: ...
-1
votes
1answer
81 views

Seurat heatmap for two conditions

I have a Seurat object of four cancers and four controls ...
1
vote
1answer
27 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 ...
1
vote
1answer
41 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 ...
1
vote
5answers
138 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
58 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
1answer
37 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 ...
0
votes
0answers
8 views

How to perform KOG classification?

I have a list of DEGs and I want to perform KOG functional analysis, I try this method but WebMGA is not updated and seems to be deprecated. Are there any other online or standalone tools available ...
0
votes
1answer
61 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 ...
0
votes
0answers
22 views

How to merge several microarray-based gene expression datasets while not excluding genes that are not present in all datasets?

In integrative gene expression analysis when two expression matrices from different datasets are merged together, they are usually merged by common genes in both expression matrices while all other ...
1
vote
0answers
24 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 ...
0
votes
0answers
53 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 ...
1
vote
1answer
192 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, ...
0
votes
3answers
64 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
1answer
69 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 ...
1
vote
0answers
17 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 ...
1
vote
0answers
110 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 ...
0
votes
0answers
36 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 ...
1
vote
1answer
112 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
0answers
43 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 ...
0
votes
2answers
31 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 ...
0
votes
1answer
29 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 ...
1
vote
2answers
289 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 ...
1
vote
3answers
986 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 ...
2
votes
1answer
78 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 ...
1
vote
1answer
369 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
1answer
36 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 "...
0
votes
0answers
30 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 ...
2
votes
1answer
54 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 ...
2
votes
1answer
603 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
vote
1answer
104 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 ...
0
votes
1answer
35 views

Filtering genes from cuffdiff results

This question has also been asked on Biostars I have run cuffdiff (with statistics turned ON) to compare two groups of samples: Control group and Late AD group. This is the command I ran to be precise:...
-2
votes
1answer
54 views

How are these definitions related to differential expression?

I have two groups of patients; for each patient I have an output file (RNA-seq) contains this information ...
2
votes
2answers
2k views

STAR quantMode vs RSEM vs Kallisto

I recently discovered this Snakemake pipeline for RNASeq that uses STAR's quantMode to quantify gene expression for DESeq2 differential expression analysis. In the past I've always seen the workflow ...
3
votes
2answers
224 views

Differential gene expression bias due to effect of an individual sample

I am analysing a human single cell RNA seq experiment, where we have 4 groups, four samples each. Data has been analysed using Seurat, with the canonical workflow. I have tried DE using various ...
0
votes
1answer
903 views

Seurat define manually a cluster and find markers

I am aware of this question Manually define clusters in Seurat and determine marker genes that is similar but I couldn't make tit work for my use case. So I have a single cell experiments and the ...
0
votes
1answer
53 views

How I deal with this expression set?

I have a gene expression raw counts like below for 16 patients ...
0
votes
1answer
89 views

How I deal with this kind of gene expression comparision

I have RNA-seq from two sequencing batches; Lab technician says that he has run the RNA expression quantification two times in bathes 1 and 2 for example tumor 1 in batch 1 and tumor 1 in batch 2 , ...
0
votes
1answer
41 views

How to perform DE analysis for each sample

I am new to R and biocondunctor. I have the normalized expression values for 20 samples for a disease and for 10 controls. I wanted to get the differential expressed values for each sample with all 10 ...
1
vote
0answers
51 views

Identification of differential genes across 8 groups [closed]

I need to identify differential genes across 8 groups. I know this can be done using DESeq2 or EdgeR as these are better suited and equipped. However, I was thinking if something like this piece of ...
0
votes
1answer
84 views

Single sample in group: normal pipeline or Kal's Z test

As stated, which one is better for differential expression analysis? When I say normal pipeline I mean limma-voom, edgeR and DESeq2 pipeline for standard analysis. Kal's z test is mentioned in this ...
1
vote
2answers
496 views

Strange p-value histogram for differential gene expression analysis

I'm trying to perform pretty standard differential expression analysis using RNA-seq data. I've used Kallisto to perform RNA quantification and am using Sleuth to perform the differential expression ...
0
votes
1answer
677 views

Why does DESeq2 convert numeric columns to factor during differential expression analysis?

I'm attempting to perform differential expression analysis using DESeq2 on a dataset of ~90 samples and ~20000 transcripts. I have a number of variables I'd like to test against, some of which are ...