Question: Comparing gene expression for individual gene between samples in R
gravatar for mrz132435
8 months ago by
mrz13243520 wrote:

Hi. I have whole exome raw count data for a bunch of samples, and I'm looking to compare gene expression for a single gene between samples, e.g., between samples that have some type of mutation in said gene and samples that don't. I'd like to do this analysis over a set of genes, but each one individually (that is, I only care about how a mutation in gene X affects expression in gene X, but not how it might affect expression elsewhere). The issue is, I'm not sure how to do this in DESeq2, which seems to be the most popular gene expression package for R. The main DEseq2 function requires a design in which each sample is assigned a condition (e.g., treated vs. untreated) and then applies its negative binomial model to the unnormalized counts, but it would be prohibitively time-consuming to rerun DESeq2 for each gene, using each gene's mutation status as the condition for each run.

So, what would be the best way to do just normalize my count data so that I can compare individual genes (either with DESeq2 or another package)? Also, any comments on what the optimal statistical approach would be for comparing expression between two groups of samples for an individual gene are also appreciated. Thanks.

rna-seq R • 197 views
ADD COMMENTlink modified 8 months ago by swbarnes28.1k • written 8 months ago by mrz13243520
gravatar for swbarnes2
8 months ago by
United States
swbarnes28.1k wrote:

The size normalization step of DESeq can be done with a design of ~1. It doesn't take the design into account at all.

Unless you have an enormous number of samples, it doesn't take all that long to run DESeq multiple times. Or run it once, and ask for the result of multiple contrasts.

ADD COMMENTlink written 8 months ago by swbarnes28.1k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1355 users visited in the last hour