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.
79
questions
0
votes
0answers
13 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
42 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
55 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
15 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 ...
0
votes
0answers
52 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
34 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
48 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
23 views
particular heatmap combinning to barplot for up and down-regulated genes
I have a list of genes (regulated up and down) by Deseq2 . I'd like to generate a heatmap of a particular combination for the barplot like in this graph :
which I read in this paper: https://www....
0
votes
0answers
39 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
27 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
19 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
99 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
155 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
29 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
56 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
31 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
26 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 ...
0
votes
0answers
10 views
Cuffdiff multiple comparison
I have four groups G1,G2,G3 and G4 all in triplicates, I want to make the comparison of differential expression between the four groups together used cuffdiff following:
...
1
vote
1answer
28 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
183 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
93 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 ...
-1
votes
1answer
21 views
Filtering genes from cuffdiff results
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
47 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
938 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
147 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
341 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
0answers
30 views
Filtering criteria for selecting DEGs from Monocle3 regression analysis with categorical terms
I'm testing the package Monocle3 (version 0.2.0) with the following code. The main aim of this analysis is to identify differentially expressed genes for every time ...
0
votes
1answer
52 views
How I deal with this expression set?
I have a gene expression raw counts like below for 16 patients
...
0
votes
1answer
87 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
38 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
46 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
0answers
55 views
How could I find Differential Genes in Expression Matrix File by R?
I tried to get gset of my dataset by getGEO function, but unfortunately it didn't work( I installed GEOquery package) and i get errors like this:
Found 1 file(s)
GSE69606_series_matrix.txt.gz
trying ...
0
votes
1answer
54 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
283 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
278 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 ...
1
vote
0answers
80 views
Gene ratio as imput in limma
I have a data frame with gene-expression ratio. Is it possible to input this into limma/voom to find signinficannt gene-ratios between groups of samples?
my data:
...
3
votes
1answer
77 views
Merging transcriptomes coming from different experiments
I'm planning to build a transcriptome by pooling all existing transcriptomes in SRA for a non-model species (which has no reference genome) to study differentially expressed genes and the like. It ...
3
votes
1answer
197 views
load the phenotype data for ballgown
I am trying to reproduce the work of this paper
[1], and I have run StringTie successfully, but after that I have to run Ballgown but could not understand this command:
...
0
votes
0answers
24 views
How to retrieve the best-scoring Trinity isoform from blastx results
I'm working on a somewhat unusual transcriptome focused on a killer X chromosome. This chromosome has some new genes, with expression levels quite different from their autosomal paralogs.
I've found ...
0
votes
1answer
64 views
I need some tips and suggestions for further analysis of NGS expression data (log2cpm)
I am a PhD student who inherited some log2cpm data of expression data from bulk kidney tissue from a UUO(unilateral urethral obstruction) experiment that tests a new drug. The sample material consists ...
0
votes
0answers
53 views
Common Mouse and human DEG analysis
I have DEGs from human and mouse (equivalent of a human disease model), and aim to generate a logCPM correlation plot between the overlapping mouse and human gene (275 genes). Prior to joining, mouse ...
5
votes
1answer
570 views
How to normalise scRNASeq data for differential expression analysis
I wish to perform differential expression analysis for cluster-specific gene expression in single-cell data (with a tool such as MAST or SCDE).
I have data for 3 biological replicates. I performed ...
5
votes
2answers
4k views
What does an FDR value of 1 in RNA-seq mean?
I am looking at the supplemental data from the paper "An allelic series of miR-17 ā¼ 92-mutant mice uncovers functional specialization and cooperation among members of a microRNA polycistron" which ...
1
vote
1answer
304 views
Compare clusters from different datasets
I have two scRNAseq datasets, and I select a cluster from one of them, and another cluster from the other. I want to find differentially expressed genes between the ...
1
vote
1answer
595 views
DESeq2 complicated design - effect of replicated samples
I have RNAseq data from a relatively complicated experimental design with variables = genotype, treatment, time, and batch. I have 2 biological replicates for each genotype/condition, however ...
3
votes
1answer
136 views
Change sample names while using cummeRbund
I want to analyze RNA seq data using R package cummeRbund. Using this package I have a problem changing sample names. When I run this command in cummeRbund
...
2
votes
3answers
206 views
Selection of differential expressed genes
I'm working with RNA-seq data. I have 40 tumor samples and 5 Normal samples. Differential analysis with Deseq2 based on Fold change > 1.2 and alpha < 0.05 gave ...
4
votes
1answer
114 views
Is removing samples based on clustering for downstream analysis a right choice?
I'm using TCGA Lung cancer data. I'm interested in doing differential analysis between Lung vs Normal. Before DEA, to check the distance between each pair of samples I plotted an MDS plot:
In this, I ...
3
votes
2answers
227 views
Detecting differentially expressed genes with foldchange >= 2 and FDR < 0.05
I'm using edgeR for differential analysis. Using glmTreat function I'm detecting differentially expressed genes between Tumor and Normal. I have set the arguments like below:
...
2
votes
2answers
41 views
Can I ask DESeq2 which variable alone explains which gene's behavior the best?
Imagine I have a dataset of highly periodic gene expression, taken at densely spaced timepoints. I know that the expression is periodic, but different genes have different offsets.
It is trivial to ...
