Question

Using contrast argument of DESeq2 to estimate linear combination coefficients

0

Entering edit mode

sebastiano.curreli • 0

@sebastianocurreli-9937

Last seen 7.6 years ago

Dear Community,

I'm characterizing transcriptional responses to transcription factors (TFs) overexpression using RNA-seq.

In brief, I have 5 groups of samples (x, y, z, xyz, Ctrl), where I provide either x, y or z TFs individually, or simultaneously using a multicistronic vector (xyz); Ctrl represents a control vector infection. Each group has 5 replicates. Data has been filtered for 0 counts removal.

So far, I've used the contrast argument as follows:

numerator<-c("conditionxyz")
denominator<-c("conditionx","conditiony","conditionz")

LC<-list(numerator,denominator)

DEcomposite<-results(dds, contrast=LC, listValues = c(1,-1/3), alpha=0.05, lfcThreshold = 1, altHypothesis = "greaterAbs", independentFiltering = T)

This allowed me to identify the DE genes between the average of the transcriptional effects of x, y, z and their simultaneous effect (xyz).

Currently, I'm trying to estimate the "dose" of x,y and z that better resemble the transcriptional profile of simultaneous xyz overexpression. Which basically, translates to identify the coefficients for the results function argument listValues that minimize the differentially expressed genes.

Does anybody know how to approach this problem?

Thanks in advance!

deseq2 contrast • 1.5k views

ADD COMMENT • link updated 7.9 years ago by Michael Love 43k • written 7.9 years ago by sebastiano.curreli • 0

score 1 · Answer 1 · 2017-02-14

1

Entering edit mode

Michael Love 43k

@mikelove

Last seen 2 days ago

United States

I don't see how to do this within the DESeq2 context. Let me see if I can turn it into a mathematic formulation:

Suppose you collapse the 5 replicates into an average vector X, Y, Z, and XYZ (these are vectors over genes).

You want to find alpha, beta, gamma > 0 to minimize:

d(alpha X + beta Y + gamma Z, XYZ)

where d(., .) is some distance.

ADD COMMENT • link 7.9 years ago Michael Love 43k

0

Entering edit mode

Thank you Michael.

Your formulation describes somehow what I'm looking for, the problem is that the "distance" I'm willing to use is the number of DE transcripts identified with DESeq2.

I thought to write a function which iterates n times the results extraction while modifying the contrast definition, such as numerator remains constant, while the members of denominator are multiplied with set of alpha, beta, gamma coefficients randomly chosen between 0 and 1. This to populate a list of n elements containing the number of DE transcripts, and the coefficients used at each iteration step. However, I've realized that this isn't easily implemented, because we can just provide one single coefficient respectively for the numerator and the denominator arguments through listValues.

Do you envision any strategy to get around this? Or do you think that using DESeq2 in this way is just hopeless?

ADD REPLY • link 7.9 years ago sebastiano.curreli • 0

1

Entering edit mode

I personally would prefer to use a distance rather than a testing framework for this. You might find the DESeq2 transformations useful though.

ADD REPLY • link 7.9 years ago Michael Love 43k