Tag: deseq2

InteractiveComplexHeatmap on DESeq2 object with more than 2 groups

InteractiveComplexHeatmap on DESeq2 object with more than 2 groups 1 Hello all, I’m writing with the hope someone can clarify a doubt I have about the differential heatmap generated by the package InteractiveComplexHeatmap via the simple command interactivate(dds). I read the package documentation at bioconductor.org/packages/release/bioc/html/InteractiveComplexHeatmap.html, but I couldn’t find the…

Continue Reading InteractiveComplexHeatmap on DESeq2 object with more than 2 groups

tReasure: R-based GUI package analyzing tRNA expression profiles from small RNA sequencing data | BMC Bioinformatics

tReasure (tRNA Expression Analysis Software Utilizing R for Easy use) is a graphical user interface (GUI) tool for the analysis of tRNA expression profiles from deep-sequencing data of small RNAs (small RNA-seq) using R packages. The whole analysis workflow, including the uploading of FASTQ files of small RNA-seq, quantification of…

Continue Reading tReasure: R-based GUI package analyzing tRNA expression profiles from small RNA sequencing data | BMC Bioinformatics

DESeq2 how to specify contrast to test difference of differences

I am trying to take the “difference of differences” in contrasts from two factors (sex and group). We have male and female animals (sex factor) that were untrained or trained for 1, 2, 4, or 8 weeks (group factor, i.e., “control”, “1w”, “2w”, “4w”, “8w”). I want to know the…

Continue Reading DESeq2 how to specify contrast to test difference of differences

Error in SummarizedExperiment

I have installed DESeq2 version 1.36.0 samples <- colnames(txi$counts) group <- as.factor(c(“control”,”control”,”control”,”control”,”control”,”diet”,”diet”,”diet”,”diet”,”diet”, “control”,”control”,”control”,”control”,”control”,”diet”,”diet”,”diet”,”diet”,”diet”,”diet”)) coldata <- data.frame(samples, group, stringsAsFactors = F) coldata <- coldata[,c(“samples”,”group”)] coldata$samples <- factor(coldata$samples) coldata$group <- factor(coldata$group) rownames(coldata) <- sub(“fb”, “”, rownames(coldata)) all(rownames(coldata$samples) %in% colnames(txi)) all(rownames(coldata) == colnames(txi)) TRUE library(DESeq2) ddsTxi <- DESeqDataSetFromTximport(txi, colData = coldata, design =…

Continue Reading Error in SummarizedExperiment

Genes looking differential abundant are not accoring to DESeq2

Genes looking differential abundant are not accoring to DESeq2 0 I have a metagenomic dataset crossing three time points from which I have mined CAZymes and am using DESeq2 to identify differentially abundant CAZymes from using trimmed mean depth generated my CoverM (very similar to Q2Q3 contig coverage). From this…

Continue Reading Genes looking differential abundant are not accoring to DESeq2

Extremely different results for both EdgeR and DESeq2 analysis

Extremely different results for both EdgeR and DESeq2 analysis 1 @373f98d7 Last seen 23 hours ago Singapore Dear all, Upon comparing my results for the analysis between DESeq2 and EdgeR, I have realized that the 2 results obtained after DEG analysis are extremely different from each other. The thresholds I…

Continue Reading Extremely different results for both EdgeR and DESeq2 analysis

RNA interference triggers that target SARS-CoV-2 genome

Coronavirus disease 2019 (COVID-19) vaccines have played a critical role in reducing transmission of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) transmission. However, with emerging reports of waning vaccine efficacy, there remains an urgent need to develop prophylactic measures against COVID-19. In a recent study published on the bioRxiv*…

Continue Reading RNA interference triggers that target SARS-CoV-2 genome

Can Differential Isoform expression analysis can be performed using DESeq2 package

Can Differential Isoform expression analysis can be performed using DESeq2 package 0 @03ddb485 Last seen 9 hours ago India Hello, I am want to perform differential isoform expression (DIE) analysis for RNAseq data from human. Can I use DESeq2 for this by inputting the transcript level abundance and getting differential…

Continue Reading Can Differential Isoform expression analysis can be performed using DESeq2 package

Counts from recount3 for DESeq2 analysis

Hi, As the topic title suggests, I want to use RSE objects from recount3 for differential expression analysis with DESeq2. This is an ongoing project that I took over. Looking at what has already been done, there are two points I wonder about : 1 – I found that the…

Continue Reading Counts from recount3 for DESeq2 analysis

deseq2 problem

deseq2 problem 0 Hi I am trying to draw a PCA plot with DESeq2 but somehow I cannot use DESeq2 functions. It is a really simple code i wil be pasting below. > transform <- DESeq2::rlog(eliminated_data, blind = TRUE) Error in (function (classes, fdef, mtable) : unable to find an…

Continue Reading deseq2 problem

scrnaseq – Normalization methods to combine scRNA-seq experiments with different sequencing depths

I don’t think you need to complicate the idea of normalisation by introducing machine learning classifiers as a necessary component. Normalisation is common when comparing different datasets for all differential analysis. If you have single cell data, have a look at integration techniques in the Seurat workflows: satijalab.org/seurat/articles/integration_rpca.html If you’ve…

Continue Reading scrnaseq – Normalization methods to combine scRNA-seq experiments with different sequencing depths

DESEQ2 Results log2fold and pvalue changes after removing genes

DESEQ2 Results log2fold and pvalue changes after removing genes 1 @bine-23912 Last seen 9 hours ago UK Dear all, we have been discussing internally but couldnt find an answer, maybe you could help us. The following happened: I run my DESEQ2 analysis with the full dataset and got genes with…

Continue Reading DESEQ2 Results log2fold and pvalue changes after removing genes

different results between DESeq model with multiple groups or with specific groups

different results between DESeq model with multiple groups or with specific groups 1 @b295d7f1 Last seen 3 hours ago Italy Hi! I am working with DESeq2 to perform a differential expression analysis between different treatments.I have 4 conditions and 4 four biological replicates for each conditions. I’ve performed differential expression…

Continue Reading different results between DESeq model with multiple groups or with specific groups

Tximport in usegalaxy

Tximport in usegalaxy 0 Devon Ryan: Please help in resolving this issue. How to use tximport in usegalaxy to convert transcript ID(DESEQ2-SALMON) to gene ID. I want to get gene ids from the results of deseq2(salmon) . Which GTF should be used for tximport. Iam getting the following error in…

Continue Reading Tximport in usegalaxy

About outliers and non -separated samples in PCA

About outliers and non -separated samples in PCA 0 Hi all, I have plotted PCA for my samples(Tumor and Normal) in some cancer types. I have used the HTSeq-counts data from TCGA. Then I’ve normalized them by DESeq2 and the total normalized counts are in cnt dataframe. Head of cnt:…

Continue Reading About outliers and non -separated samples in PCA

How to merge miRDeep2 counts for DeSeq2 analysis?

How to merge miRDeep2 counts for DeSeq2 analysis? 0 Hi, I am trying to identify DE miRNAs using small RNA-Seq data. I ran miRDeep2 tools mapper.pl an also miRDeep2.pl scripts and got counts for my cases. How can I merge count.cvs files of multiple cases for DE miRNA detection using…

Continue Reading How to merge miRDeep2 counts for DeSeq2 analysis?

What is the relationship between the value of the bigwig file generated by MACS2 and the read count?

What is the relationship between the value of the bigwig file generated by MACS2 and the read count? 0 What is the relationship between the value of the bigwig file generated by MACS2 and the read count? When using Deseq2 or edgeR to identify the differential peak, because the value…

Continue Reading What is the relationship between the value of the bigwig file generated by MACS2 and the read count?

Three Differential Expression Analysis Methods for RNA Sequencing: limma, EdgeR, DESeq2

A detailed protocol of differential expression analysis methods for RNA sequencing was provided: limma, EdgeR, DESeq2. Three differential expression analysis methods for RNA sequencing:limma, EdgeR, and DESeq2. Open the RStudio program and then load R file, DEGs. The file can be acquired from supplementary files.One. Downloading and pre-processing of data.1.1….

Continue Reading Three Differential Expression Analysis Methods for RNA Sequencing: limma, EdgeR, DESeq2

Error in DESeq2

Error in DESeq2 1 @4bede5c6 Last seen 11 hours ago India Hi, All I am new to Rstudio , i Have to study expression data for the following But my DESeq2 data is having an error of Error in h(simpleError(msg, call)) : error in evaluating the argument ‘x’ in selecting…

Continue Reading Error in DESeq2

How to identify differential enhancers?

How to identify differential enhancers? 0 Deseq2 or edgeR only allow the input of reads count to identify differential regions, but I get less differential enhancer with reads count, can I use the signal from the bigwig file as input to deseq2 or edgeR to identify the differential enhancer? EdgeR…

Continue Reading How to identify differential enhancers?

DESeq2 (or EdgeR) Exploratory Analysis with no Replicates

DESeq2 (or EdgeR) Exploratory Analysis with no Replicates 1 My pipeline so far is hisat2->featureCounts->DESeq2. I have generated heatmaps after rlog and log2 transformation of the genes with the most variance, which is somewhat meaningful. What I really want to do is compare everything to the control sample and take…

Continue Reading DESeq2 (or EdgeR) Exploratory Analysis with no Replicates

Advice on gene expression method through heatmap

Advice on gene expression method through heatmap 0 Hello! I want to kindly get advice on a method I am using for showing gene expression on a heatmap. I have three groups: Negative Control (4 samples) Control with injury (4 samples) Treatment with injury (4 samples) I took the intersecting…

Continue Reading Advice on gene expression method through heatmap

deseq2 datasetfrommatrix

deseq2 datasetfrommatrix 1 @235c1561 Last seen 1 day ago United Kingdom countsTable <-read.csv(“counts.csv”, header = T, stringsAsFactors = FALSE) head(countsTable) sample_info <-read.table(“sample.txt”, header = T, stringsAsFactors = FALSE) #rownames(countsTable) <-countsTable[,27902] #countsTable <-countsTable[,1] class(sample_info) sample_info$Type <- as.factor(as.character(sample_info$Type)) library(DESeq2) rownames(countsTable) <-countsTable[,1] countsTable <-countsTable[,-1] (dds<-DESeqDataSetFromMatrix(countsTable[,c(1:128)],colData = sample_info,design = ~Type)) i get this error…

Continue Reading deseq2 datasetfrommatrix

DESeq2: Error in checkFullRank(modelMatrix) | Tech Programing

I want to build a DEseq2 model matrix using this coldata but it shows this error: Error in checkFullRank(modelMatrix) :the model matrix is not full rank, so the model cannot be fit as specified.Levels or combinations of levels without any samples have resulted incolumn(s) of zeros in the model matrix….

Continue Reading DESeq2: Error in checkFullRank(modelMatrix) | Tech Programing

TCGA transcriptome data to R (DESeq2)

This seems to be frequently asked question, so here is a robust method to fully recapitulate the counts given by TCGA and port it to DESeq2. Why the long way? Tanya and I noticed via TCGA-Biolinks and Firehose did not generate the full count matrix. ~5-10% of genes were missing…

Continue Reading TCGA transcriptome data to R (DESeq2)

Isolate contrasts using DESeq2 when use a combination of categorical and continuous variables

Dear community members, I am currently using DESeq2 to analyze RNASeq data from plant roots that were colonized by a fungus. But I am facing some problems regarding contrast with a continuous variable. If you have any suggestions or comments it would be great if you can let me know….

Continue Reading Isolate contrasts using DESeq2 when use a combination of categorical and continuous variables

need some help on how to use DESeq2 for TCGA data

need some help on how to use DESeq2 for TCGA data 0 Hello, I am sorry for this newbie question, but I spent all morning trying to find it out but can’t find a clear answer anywhere. I want to normalise RNA seq data from TCGA using DESeq2. I use…

Continue Reading need some help on how to use DESeq2 for TCGA data

Error with hicdcdiff function in HiCDCPlus package

I am trying to find differential interactions from my allvalidpairs data. I was able to create an index file with the union of significant interactions but when proceeding to the next step to perform differential analysis with hicdcdiff function I am getting the following error. Error in gi_list_read(path.expand(prefix)) : Some…

Continue Reading Error with hicdcdiff function in HiCDCPlus package

Paired sample design in DESeq2

I am analyzing a dataset that is looking into the effects of age on cancer using patient samples within one cancer type. In short, we are interested in finding genes involved in tumorigenesis that are altered by age. Each patient is classified as young or old, and had tumor tissue…

Continue Reading Paired sample design in DESeq2

Map mouse RNA-seq Ensembl gene ids to “org.Mm.eg.db”, but get some “UNKNOWN” gene symbol, like “NA..25729”, is it normal?

Map mouse RNA-seq Ensembl gene ids to “org.Mm.eg.db”, but get some “UNKNOWN” gene symbol, like “NA..25729”, is it normal? 0 Hi I was trying to map mouse RNA-seq Ensembl gene ids to “org.Mm.eg.db” database, but get some “UNKNOWN” gene symbol, like “NA..25729”, is it normal? How to explain it? >…

Continue Reading Map mouse RNA-seq Ensembl gene ids to “org.Mm.eg.db”, but get some “UNKNOWN” gene symbol, like “NA..25729”, is it normal?

DEseq2 replicate normalization

There are maybe two questions going on here. Normalizing by “replicate” Unless there is something common among the different replicates in your experiment, I’m not sure that including the replicate term in your model as you are doing and trying to normalize is possible/reasonable. Allow me to outline a few…

Continue Reading DEseq2 replicate normalization

Deseq2 dataframematrix error

Error in DESeqDataSet(se, design = design, ignoreRank) : some values in assay are negative I have raw counts of rna seq data with approxx 68 samples, im getting error in deseq dataset can you guys help me countsTable <-read.csv(“counts.csv”, header = T, stringsAsFactors = FALSE) head(countsTable) sample_info <-read.table(“sample.txt”, header =…

Continue Reading Deseq2 dataframematrix error

Installing packages by bioconductor version 3.14

Installing packages by bioconductor version 3.14 1 @dfefb504 Last seen 5 hours ago Vienna I installed the latest R version 4.1.2, but I got problems with installing the latest bioconductor version 3.14. I can not install any package by bioconductor. Can anyone tell me the reason? THANKS! > if (!requireNamespace(“BiocManager”,…

Continue Reading Installing packages by bioconductor version 3.14

values (NA) in p value Deseq2 (reopen)

Hello how are you? I reopen this question because the following has happened: I am doing a differential expression exercise using the hisat2, stringie & DESeq2 workflow. Finally I use the python prepDE.py script recommended in the StringTie manual to extract the counts. So far so good, I have rows…

Continue Reading values (NA) in p value Deseq2 (reopen)

Help with replicating results of differential expression analysis in DESeq2 : bioinformatics

**Also posted this on Biostars** Hi everyone! I’m fairly new to bioinformatics (self-taught for my MSc) and I’m trying to replicate a study using publicly available transcriptomic data (GSE107934). I’m struggling to get the same results as the authors. I’m following the study’s methods and using DESeq2 to conduct the…

Continue Reading Help with replicating results of differential expression analysis in DESeq2 : bioinformatics

Differential expression analysis – issue with replicating results

Hi Biostars! I’m fairly new to bioinformatics (self-taught for my MSc) and I’m trying to replicate a study using publicly available transcriptomic data (GSE107934). I’m struggling to get the same results as the authors. I’m following the study’s methods and using DESeq2 to conduct the differential expression analysis. Most of…

Continue Reading Differential expression analysis – issue with replicating results

Errors when ajusting pvalue on R

Errors when ajusting pvalue on R 2 Hello. I’m trying to ajust the pvalues from the file I’ve got from Deseq2. But I keep getting error mesages: p.adjust(dds, method = “bonferroni”) Error in as.vector(x, mode = “numeric”) : no method to coerce this S4 class to vector Then I tried…

Continue Reading Errors when ajusting pvalue on R

Alignment of de novo assembled transcripts to reference transcriptome?

Alignment of de novo assembled transcripts to reference transcriptome? 0 Hi all, I started with 40 samples of raw (but trimmed) RNASeq samples. I used these as inputs for SPAdes-rna and Trinity, to assemble them without a reference. I now have 80 assembled transcriptomes (?), and I tried to follow…

Continue Reading Alignment of de novo assembled transcripts to reference transcriptome?

Normalization by variance stabilizing transformation VST

Normalization by variance stabilizing transformation VST 1 @6c372dab Last seen 8 hours ago Sweden Hello! I am a bit confused about the normalization performed by the DESeq2 varianceStabilizingTransformation() and vst() functions in addition to the actual variance stabilization. My understanding is that the normalization by division by size factors (which…

Continue Reading Normalization by variance stabilizing transformation VST

Time series RNA seq design query

Hello, I am working with Deseq2 on RNA seq time-series data. I have RNA seq data at baseline before any intervention (timepoint0), RNA seq data with two different interventions (Placebo, and Condition) at TP1 and TP2. I went through the time-series tutorial for Deseq2 and don’t think the same example…

Continue Reading Time series RNA seq design query

Error when running ‘DESeqDataSetFromMatrix’ from Deseq2

Error when running ‘DESeqDataSetFromMatrix’ from Deseq2 1 Hello. I’m having problem at the very beggining of the DESeq2 pipeline. When I try to run the command dds <- DESeqDataSetFromMatrix(countData=countData, colData=metaData, design=~dex, tidy = TRUE) I keep getting the error mesage: Error in `.rowNamesDF<-`(x, value = value) : duplicate ‘row.names’ are…

Continue Reading Error when running ‘DESeqDataSetFromMatrix’ from Deseq2

Biostar DESEQ2

Showing : DESEQ2 • reset 2 days ago • updated 1 day ago adR &utrif; 10 10 days ago • updated 8 days ago adeler001 • 0 15 days ago • updated 9 days ago Cheryl • 0 updated 17 days ago by swbarnes2 &utrif; 920 • written 18 days…

Continue Reading Biostar DESEQ2

Treat time variable as continuous or discrete?

Hello, I have two questions regarding an RNA-seq experiment with 10 time points and 20 replicates per time point. So far I’ve run an LRT using the following models: design = ~ timepoint reduced = ~ 1 Throughout the experiment, RNA was extracted at days 3, 4, 15, 16, 44,…

Continue Reading Treat time variable as continuous or discrete?

DESeq2 input file from featureCount

DESeq2 input file from featureCount 0 Hello. I am analysing RNA-Seq data for 8 samples (2 treatments, 4 repetition each). I just finished the gene count using featureCount (Subread). How do I generate a matrix to use on DeSeq2 from the files (.RDS) genereted from featureCount? Do I just create…

Continue Reading DESeq2 input file from featureCount

Is it not possible to use LFCshrink on the “intersect” coefficient?

Is it not possible to use LFCshrink on the “intersect” coefficient? 0 I have been trying to use the LFCshrink function on the intersect coefficient. I thought this would be possible on any coefficents from resultsNames() but it gives an error: Error in apeglm::apeglm(Y = Y, x = design, log.lik…

Continue Reading Is it not possible to use LFCshrink on the “intersect” coefficient?

normalizing EXOM-seq read using DESeq2

normalizing EXOM-seq read using DESeq2 0 @do-it-23093 Last seen 9 hours ago Germany, München Hi, I have exom read aligned to hg38. I used featureCounts to quantify the reads. Therefore, I would like to normalize the reads in the different sample using DESeq2. May I know if there especial consideration…

Continue Reading normalizing EXOM-seq read using DESeq2

just trying to launching my DE-seq2 by csv but isn’t work…

hello sir, im korean and currently work in dpt, bioenginearing… first of all, my major is basically experiment and im suddenly transferred into Industry-University Cooperation Foundation. its my first time coding r script. so i have to run DE-seq2 so that make my MA-plot(and have no idea any of r…

Continue Reading just trying to launching my DE-seq2 by csv but isn’t work…

no DE genes found in subset

Troubleshooting DESeq2 Analysis: no DE genes found in subset 1 @8e2958c3 Last seen 5 hours ago United States Hi, I’m quite new to R and DESeq2 and need some help troubleshooting! I’m working on analyzing a dataset comparing 3 groups: patients before and after a treatment against a healthy control…

Continue Reading no DE genes found in subset

Read Counts from miRDeep2 for posterior DE analysis

Hello everyone, I would appreciate some help with my data: I’ve been analysing some NGS miRNA data in order to do Differential Expression profiling. Using the miRDeep2 program and algorithm, I’m in trouble when dealing with results to make a proper table in order to use it as input for…

Continue Reading Read Counts from miRDeep2 for posterior DE analysis

GOStats “genes being tested do not have corresponding GO terms”

Hi all, Has anyone ever had problems with the org.Hs.eg.db database not picking up any GO classifications for a list of genes? For example, I have a list of 21 genes, and when I try referencing the org.Hs.eg.db database, I get “genes being tested do not have corresponding GO terms”…

Continue Reading GOStats “genes being tested do not have corresponding GO terms”

Normalizing exonic reads and then log transforming or variance stabilising

Normalizing exonic reads and then log transforming or variance stabilising 0 @a-14337 Last seen 38 minutes ago United Kingdom Hi, I have raw counts per exon matrix and want to normalize and do a VST or log2 transform without running any differential expression (for downstream usage separate to DE), so…

Continue Reading Normalizing exonic reads and then log transforming or variance stabilising

Genome sequencing of the multicellular alga Astrephomene provides insights into convergent evolution of germ-soma differentiation

1. Smith, J. M. & Száthmary, E. The Major Transitions in Evolution (Freeman, 1995). Google Scholar  2. Queller, D. C. Relatedness and the fraternal major transitions. Philos. Trans. R. Soc. Lond. B Biol. Sci. 355, 1647–1655 (2000). CAS  PubMed  PubMed Central  Google Scholar  3. Grosberg, R. K. & Strathmann, R….

Continue Reading Genome sequencing of the multicellular alga Astrephomene provides insights into convergent evolution of germ-soma differentiation

Question: DESeq2 multiple factor design

I have an experiment where I got 2 factor, the first one is the cell line, and the second one is the treatment. In this, one cell line is derived from the other, the derived cell line received a inducible gene. I then treated both cell lines with the inducer…

Continue Reading Question: DESeq2 multiple factor design

Using normalised counts for pairwise gene correlation (pearson) across samples

Using normalised counts for pairwise gene correlation (pearson) across samples 0 Hi, I have carried out gene pairwise correlations between every pair of genes in my RNAseq data. I have done this with non-log transformed normalised counts generated from estimatesizefactors in DESeq2, as well as with VST counts and log2…

Continue Reading Using normalised counts for pairwise gene correlation (pearson) across samples

Filtering Gene by mean expression

Filtering Gene by mean expression 0 Dear All, Sorry for another post. I have RNAseq data and used DESeq2 pipeline for differntially expressed genes. Now I am doing WGCNA and the author suggest to filter genes based on mean expression rather than variance. I am wondering whether anyone has ever…

Continue Reading Filtering Gene by mean expression

Elastic Net for RNAseq data

Elastic Net for RNAseq data 1 @1c386ab3 Last seen 1 hour ago Canada I have RNAseq results from NovaSeq and after doing some differential gene expression with DESeq2, I will like to explore the Elastic Net model to select features that differentiate my two conditions. However, I haven’t been able…

Continue Reading Elastic Net for RNAseq data

Deseq2 model matrix error

Deseq2 model matrix error 1 I have trying to do Deseq2 analysis for mainly mutants vs controls. I have two covariates influencing the dataset – Gender and Genotype. This is how my sample annotation file looks – sample condition genotype gender SK01 mutant hr F SK02 mutant hr F SK03…

Continue Reading Deseq2 model matrix error

DESeq2 – VST data

DESeq2 – VST data 1 @3a675fdc Last seen 4 hours ago United States I have been using DESeq2 package for RNA-sq data analysis and really like the VST data in log2 units. But was unsure about usage of VST data for certain analyses. Specifically, can the VST data be used…

Continue Reading DESeq2 – VST data

How to define the gene length for RPKM calculation

How to define the gene length for RPKM calculation 4 Hi guys, I would like to calculate the RPKM of my RNA seq experiment. To do this, as from the formula, I need to know the gene length. My starting point are the row reads (single end) counts resulting from:…

Continue Reading How to define the gene length for RPKM calculation

Difference between two groups of conditions

Deleted:DESeq2: Difference between two groups of conditions 0 Hello. I intend to perform three contrasts including: i) mock.24h vs mock.4h; ii) S2.24h vs S2.4h and iii) the difference between these two comparisons (S2.24h vs S2.4h) – (mock.24h vs mock.4h). The design of my colData is as follows condition time group…

Continue Reading Difference between two groups of conditions

Calculating gene expression dominance coefficient via DESeq2

I am interested in calculating dominance coefficients of gene expression. For this we have RNAseq data consisting of two homozygous lines (i.e. Genotype AA and BB, respectively) and of the heterozygous crosses between them (i.e. Genotype AB and BA, respectively). We have data for both sexes in each group and…

Continue Reading Calculating gene expression dominance coefficient via DESeq2

Retrieving and interpreting specific contrast results from an interaction experiment in DESeq2

I have an experiment with 2 different groups, and several different time points (for simplicity, let’s say it is 2 time points), and I’m trying to find the best way to model both group and time effects. Note that time is a categorical variable. # in dds <- makeExampleDESeqDataSet(n=100,m=32) dds$time…

Continue Reading Retrieving and interpreting specific contrast results from an interaction experiment in DESeq2

S4Vectors, strange metadata behavior

S4Vectors, strange metadata behavior 0 @hauken_heyken-13992 Last seen 12 hours ago Bergen Is this internal logic intentional? makes no sense to me. # Goal: insert b or bmeta in index 1 of a # Will fail if metadata is of type GRangesList and length > 1. library(GenomicRanges) b <- bmeta…

Continue Reading S4Vectors, strange metadata behavior

betaPrior and lfcShrink

Hi, I am using DESe2 (v 1.33.4) to compare three shade avoidance phenotype in Arabidopsis. I observed that the betaPrior in current version of DESeq2 was set to FALSE. If this is the case, we can compare C vs A and B vs A but we can’t compare C vs…

Continue Reading betaPrior and lfcShrink

Annotating DESeq2 Differential Expression Results

Annotating DESeq2 Differential Expression Results 1 @adeler001-21743 Last seen 11 minutes ago Canada Hello I am trying to annotate a DESeq2 significantly differentially expressed genes results table with the gene start (bp) , gene end (bp) and SNP ID info using the biomaRt package on RStudio, when I look at…

Continue Reading Annotating DESeq2 Differential Expression Results

filtering most variable genes

filtering most variable genes 1 Dear All, I am trying to filter most variable genes for my specific analysis. I have normalized count from DEseq2 attached here with row for genes and column for sample ID. I have found code chunk in Biostar with the following **data$variance= apply(data, 1, var)…

Continue Reading filtering most variable genes

How to get DEGs compare one sample to another?

How to get DEGs compare one sample to another? 0 I have a count matrix data (23 samples with different condition), and coldata is shown above. I want to get DEGS like this: I tried this code: # Read data data_matrix <- read.csv(“matrix_data”, row.names=1) coldata <- read.csv(“data/coldata.csv”, row.names=1) # Creating…

Continue Reading How to get DEGs compare one sample to another?

When I change samples grouping for normalization by DESeq2 I get different results, Why?

When I change samples grouping for normalization by DESeq2 I get different results, Why? 0 Hi all I have downloaded the raw counts of RNA-Seq from TCGA for DE analysis. the number of primary tumor samples was 223 and the normal adjacent tissue was 41. I have performed the normalization…

Continue Reading When I change samples grouping for normalization by DESeq2 I get different results, Why?

Which interactions terms to choose?

Which interactions terms to choose? 0 I would appreciate your help with choosing the right structure of result names to answer my research question. I defined 2 conditions: Strain (WT is the reference level) cell-Fraction (FractionA is the reference level) I want to calculate the difference between FractionA and FractionB…

Continue Reading Which interactions terms to choose?

Which rlog should I use for DESeq2 analysis

Which rlog should I use for DESeq2 analysis 0 @rafaelsolersanblas-22935 Last seen 1 minute ago Alicante I have a question regarding the rlog normalization. I have many samples to compare, with only one factor. A vs Treat, B vs Treat, C vs Treat … So, should I put everything in…

Continue Reading Which rlog should I use for DESeq2 analysis

DEseq2 contrast in two time points

I need help with some clarifications please I have control versus treatment in two time points like > y X condition time 1 CTRL_24_hrs_replicate1 control 24 2 CTRL_24_hrs_replicate2 control 24 3 CTRL_24_hrs_replicate3 control 24 4 treatment_24_hrs_replicate1 t 24 5 treatment_24_hrs_replicate2 t 24 6 treatment_24_hrs_replicate3 t 24 7 CTRL_48_hrs_replicate1 control 48…

Continue Reading DEseq2 contrast in two time points

Principal Component Analysis Rna Seq

Listing Results Principal component analysis rna seq Genomatix Principal Component Analysis For RNASeq Data Preview 2 hours agoThis is explained in detail on “RNA–Seq workflow: gene-level exploratory analysis and differential expression”. The matrix of raw counts is input to the DESeq2 rlog function and the resulting transformed matrix is used…

Continue Reading Principal Component Analysis Rna Seq

95% confidence intervals for PCA in DESEQ2

95% confidence intervals for PCA in DESEQ2 2 @add481ab Last seen 1 hour ago United States Dear Help, We have used your package DESEQ2 (including the vst transform) on some RNA-seq data, in order to perform PCA analysis. We were hoping to add Monte-Carlo noise to the data in order…

Continue Reading 95% confidence intervals for PCA in DESEQ2

Multiple stages of evolutionary change in anthrax toxin receptor expression in humans

Human research participants We have complied with all relevant ethical regulations and informed consent was obtained from all participants. This work was approved by the Cornell University IRB under protocol 1506005662. Animal research This work was approved by the Cornell University IACUC under protocol 2009-0044. Welfare and handling of all…

Continue Reading Multiple stages of evolutionary change in anthrax toxin receptor expression in humans

Bioconductor – IHWpaper

DOI: 10.18129/B9.bioc.IHWpaper     This package is for version 3.11 of Bioconductor; for the stable, up-to-date release version, see IHWpaper. Reproduce figures in IHW paper Bioconductor version: 3.11 This package conveniently wraps all functions needed to reproduce the figures in the IHW paper (www.nature.com/articles/nmeth.3885) and the latest arXiv preprint available…

Continue Reading Bioconductor – IHWpaper

Do I have to separate my interest genes from my count matrix and then perform differential expression analysis for them?

Do I have to separate my interest genes from my count matrix and then perform differential expression analysis for them? 0 Hi all, I am trying to study the differential expression of my interest genes in colon cancer. First, I’ve downloaded the RNA-Seq raw counts from TCGA and have built…

Continue Reading Do I have to separate my interest genes from my count matrix and then perform differential expression analysis for them?

how to extract normalized expression from DESeq2 object (dds)?

how to extract normalized expression from DESeq2 object (dds)? 1 Hi I have multiple groups to compare using DESeq2 and used contrast argument in results function to find DE. But now i would like to extract the normalized expression based on compared group for the heatmap and the further merging…

Continue Reading how to extract normalized expression from DESeq2 object (dds)?

Gene list for GO enrichment after DESeq2. Should I use a gene list of differentially expressed genes (DEG) before or after shrinkage?

Gene list for GO enrichment after DESeq2. Should I use a gene list of differentially expressed genes (DEG) before or after shrinkage? 0 Hi, I have been starting by first steps into RNA-Seq data and I have been performing my first analysis since two months ago. So I apologise if…

Continue Reading Gene list for GO enrichment after DESeq2. Should I use a gene list of differentially expressed genes (DEG) before or after shrinkage?

DESeq2 log2Foldchange issue

Hi, when I got my DESeq2 analysis result, I found the log2Foldchange values are between -1 to 1. But when I checked the normalized read counts, I found they are hugely different between my two types of tissue (3 replicates for each type). Is that normal? My code is like:…

Continue Reading DESeq2 log2Foldchange issue

Complex Contrast in DESeq2

Hi, I am using DESeq2 to look at the effect of different diet over time. I have 3 different diet (D1, D2, D3), and 3 time points (T1, T2, T3). I am trying to create best model for this, therefore I tried dds <- DESeqDataSetFromMatrix(countData=rawCounts_totat, colData=mapping_file, design=~Diet+Time+Diet:Time) dds.1 <- DESeq(dds,…

Continue Reading Complex Contrast in DESeq2

how to combine counts and DESeq data into one CSV?

Noob question: how to combine counts and DESeq data into one CSV? 2 I am trying to replicate the format of the below file, which seems to include both the raw counts AND the DESeq data into one csv WTS1 WTS2 KOS1 KOS2 baseMean log2FoldChange lfcSE stat pvalue padj Gm6166…

Continue Reading how to combine counts and DESeq data into one CSV?

Getting samples.txt to use in the salmon to DESeq2 pipeline

Getting samples.txt to use in the salmon to DESeq2 pipeline 0 Hi all, I’d appreciate some help on the salmon-to-DESeq2 pipeline outlined here. It seems to require a samples.txt file, which I don’t have anywhere. I used Salmon to quantify my assembled transcripts (40 x 3 = 120 quant.sf files),…

Continue Reading Getting samples.txt to use in the salmon to DESeq2 pipeline

Assessment of CircRNA Expression Profiles and Potential Functions in Brown Adipogenesis

doi: 10.3389/fgene.2021.769690. eCollection 2021. Affiliations Expand Affiliations 1 Department of Biotechnology, College of Life Sciences, Xinyang Normal University, Xinyang, China. 2 Institute of Animal Science and Veterinary Medicine, Hainan Academy of Agricultural Sciences, Haikou, China. 3 Institute for Conservation and Utilization of Agro-Bioresources in Dabie Mountains, Xinyang Normal University, Xinyang,…

Continue Reading Assessment of CircRNA Expression Profiles and Potential Functions in Brown Adipogenesis

collapseReplicates in DESeq2 when the number of technical replicates per biosample is different

collapseReplicates in DESeq2 when the number of technical replicates per biosample is different 0 Good afternoon, I have a question about using collapseReplicates in DESeq2. As far as I understand, this function adds up the counts belonging to one biosample. I understand the meaning of this if the number of…

Continue Reading collapseReplicates in DESeq2 when the number of technical replicates per biosample is different

Bioconductor – zinbwave

DOI: 10.18129/B9.bioc.zinbwave     This package is for version 3.12 of Bioconductor; for the stable, up-to-date release version, see zinbwave. Zero-Inflated Negative Binomial Model for RNA-Seq Data Bioconductor version: 3.12 Implements a general and flexible zero-inflated negative binomial model that can be used to provide a low-dimensional representations of single-cell…

Continue Reading Bioconductor – zinbwave

Differentail expression analysis between 3 groups

Differentail expression analysis between 3 groups 0 Hi, I have a question on differential expression analysis. I have 3 groups(A, B, C) and each group has 20 samples. I want to do differential expression analysis to get genes differently expressed in A, not in B and C. I would appreciate…

Continue Reading Differentail expression analysis between 3 groups

DESeq2 interaction terms in 2x2x2 factorial design

Hello, I’m trying to analyze RNASeq data from an multifactorial experimental setup with DESeq2. In brief, the experiment is about multiple stressor effects in stream organisms, whereby different stressors (here: added sediment, increased salinity or reduced flow velocity) are applied to the study organisms, either as single stressors or in…

Continue Reading DESeq2 interaction terms in 2x2x2 factorial design

Functional profiling visualization enrichplot error: NAs introduced by coercion

Hello, I’m having trouble generating a dotplot of functionally enriched pathways across different groups. Most of them are mapped as NA via Warning message: In order(as.numeric(unique(result$Cluster))) : NAs introduced by coercion I’m trying to figure out why this happens. I’ve included the code to the original code I am adapting…

Continue Reading Functional profiling visualization enrichplot error: NAs introduced by coercion

Comparing the whole genelist by GSEA GO BP terms to adjusted pvalue genelist by DAVID GO BP terms

Comparing the whole genelist by GSEA GO BP terms to adjusted pvalue genelist by DAVID GO BP terms 0 Say I have done DESeq2 on my RNA-Seq dataset: experimental vs. control DESeq2 has a column for BH – adjusted Pvalues and I plan to take the genes with less than…

Continue Reading Comparing the whole genelist by GSEA GO BP terms to adjusted pvalue genelist by DAVID GO BP terms

Chapter 9 Differential abundance analysis

Here, we analyse abundances with three different methods: Wilcoxon test (CLR), DESeq2, and ANCOM-BC. All of these test statistical differences between groups. We will analyse Genus level abundances. We might want to first perform prevalence filtering to reduce the amount of multiple tests. In this particular dataset, all genera pass…

Continue Reading Chapter 9 Differential abundance analysis

How we can create a shiny app to run Deseq2 ? – shiny

Hello,how we can create a shiny app to run Deseq2 ?suppose i have count table and meta table as below : count_table GENE sample1 sample2 sample3 sample4 sample5 A1BG 8 8 8 7 7 A1BG-AS1 7 6 5 4 4 A1CF 6 6 6 7 6 A2M 6 6 2…

Continue Reading How we can create a shiny app to run Deseq2 ? – shiny

Bioinformatics Scientist – Job at DAWSON in Bethesda, MD

Bioinformatics Scientist Full Time Prof-Entry Bethesda, MD, US DAWSON is a Native Hawaiian Organization 8(a) small business that brings the Spirit of Aloha to our employees. As part of the DAWSON “Ohana”, you will be provided a best-in-class benefits program that strives to ensure our great people have peace of…

Continue Reading Bioinformatics Scientist – Job at DAWSON in Bethesda, MD

How to create shiny app for deseq2 object creation

How to create shiny app for deseq2 object creation 0 Hello, how we can create a shiny app to run Deseq2 ? suppose i have count table and meta table as below : count_table GENE sample1 sample2 sample3 sample4 sample5 A1BG 8 8 8 7 7 A1BG-AS1 7 6 5…

Continue Reading How to create shiny app for deseq2 object creation

Account for parent cell line when comparing isogenic mutants

Account for parent cell line when comparing isogenic mutants 1 @oliverziff-22402 Last seen 1 day ago United Kingdom Hi there, I am using DESeq2 to analyse RNAseq from mutants and control samples. I have 3 controls and 2 mutant samples from distinct individuals but in addition I have a sample…

Continue Reading Account for parent cell line when comparing isogenic mutants

Different output when generating step by step airway dataset

rnaseqGene: Different output when generating step by step airway dataset 1 @dec8f401 Last seen 13 hours ago Spain I have been following the workflow ‘RNA-seq workflow: gene-level exploratory analysis and differential expression’ and I have a problem: I tried to generate step by step the airway dataset, downloading the fastq…

Continue Reading Different output when generating step by step airway dataset

Bug#994443: [R-pkg-team] Bug#994443: Bug#994443: marked as done (r-bioc-deseq2: autopkgtest regression)

Control: reopen -1 Hi Nilesh On Fri, 29 Oct 2021 at 21:15, Nilesh Patra <nil…@nileshpatra.info> wrote: > I think this has got nothing to do with our changes in unstable — I’d be > very surprised if so, > since essentially atleast stable is disconnected, unless some other package >…

Continue Reading Bug#994443: [R-pkg-team] Bug#994443: Bug#994443: marked as done (r-bioc-deseq2: autopkgtest regression)

Help for extraction of fasta sequences

Hello everyone, I hope you are well. I am writing this post because I have a question or rather I have a problem with my workflow. Perform a workflow for RNA-seq processing as follows: quality control – Hisat2 – Stringtie – Deseq2 A simple, normal workflow that threw me important…

Continue Reading Help for extraction of fasta sequences

Bioinformatics analysis to screen DNA methylation-driven genes for prognosis of patients with bladder cancer

This article was originally published here Transl Androl Urol. 2021 Sep;10(9):3604-3619. doi: 10.21037/tau-21-326. ABSTRACT BACKGROUND: Bladder cancer (BLCA) is the most prevalent tumor affecting the urinary system, and has contributed to a rise in morbidity and mortality rates. Herein, we sought to identify the methylation-driven genes (MDGs)of BLCA in an…

Continue Reading Bioinformatics analysis to screen DNA methylation-driven genes for prognosis of patients with bladder cancer

Intersection of multiple methods for RNA-Seq differential expression: conservative or crazy?

Intersection of multiple methods for RNA-Seq differential expression: conservative or crazy? 1 Hi, For a while now, I’ve been looking at the intersect of the significant results from DESeq2 and EdgeR as my standard for determining differentially expressed genes, treating EdgeR as a filter over the DESeq2 results. I usually…

Continue Reading Intersection of multiple methods for RNA-Seq differential expression: conservative or crazy?

Analyzing counts over 3 different timepoints and 2 conditions with DESeq2

I am trying to analyze data that has three different timepoints (0, 6, and 12 hr) and two conditions (treated and control), but every attempt I make to carry this out with DESeq2 is met with error. The latest is below, which is an attempt I made using this brief…

Continue Reading Analyzing counts over 3 different timepoints and 2 conditions with DESeq2