I am still working on preparing the data for running DEXSeq. I have a question. Will DEXSeq work if there are more than two conditions, e.g. cond1,cond2,cond3.. Thanks!
Seqanswers Leaderboard Ad
Collapse
Announcement
Collapse
No announcement yet.
X
-
The original paper has both single-end and paired-end data:
RNA-seq is a powerful tool for the study of alternative splicing and other forms of alternative isoform expression. Understanding the regulation of these processes requires sensitive and specific detection of differential isoform abundance in comparisons between conditions, cell types, or tissues. W …
Comment
-
DEXSeq multiple tissues
Originally posted by areyes View PostYes it will. If you run an analysis with all your conditions, it will show those exons in which at least in one condition is different from the other conditions.
Hello,
I have 37 RNA-seq datasets for 18 human tissues - some have 3 replicates, some 2 and some do not have any (sorry to say that).
I'm thinking of a way to do differential exon usage analysis - using DEXseq.
I could not see any example in the documentation where multiple conditions are compared (eg., tissues) other than pairwise comparisons. In edgeR, I was able to use the ANOVA-like test for getting differentially expressed genes (across an average of all 18 tissues). I was wondering if there is such a test I could do with DESeq/DEXseq.
Thanks,
Manoj.
Comment
-
Originally posted by areyes View PostYes it will. If you run an analysis with all your conditions, it will show those exons in which at least in one condition is different from the other conditions.
I have RNA-Seq paired-ends data for 3 conditions: Control (4 biological replicates), Leukemic cells with mutation in gene X (6 biological replicates) and Leukemic cells without mutation in gene X (4 biological replicates).
One of the questions I try to answer is the differential exon usage induced by this mutation. I am using DEXSeq_1.8.0.
Basically, I need to compare exon level in leukemic cells with and without the mutation.
Since I have additional information with the control, I would like to output the genes that show a differential expression in at least one exon between both conditions leukemic cells with and without the mutation, with in addition the level for the controls.
To be clear, I don't want to see the genes that show a differential expression in at least one exon between conditions leukemic cells with or without the mutation and the control.
Is it possible with DEXSeq?
If not, I have 2 choices: either I "loose" the information from the control condition or I see all the changes including those that do not come from the mutation (that's what I did for now).
Thank you in advance,
Jane
Comment
-
Hi @Jane M,
You could try to change the GLM formulas to include all the variables of your analysis. First, I would recommend updating your R and DEXSeq version. if your design data frame looks like this:
Code:phenotype mutatedGene control no control no control no control no leukemia yes leukemia yes leukemia yes leukemia yes leukemia yes leukemia yes leukemia no leukemia no leukemia no leukemia no
Code:formulaFullModel <- ~ sample + exon + phenotype:exon + mutatedGene:exon formulaReducedModel <- ~ sample + exon + phenotype:exon
Comment
-
Originally posted by areyes View Postyou could try to include the formulas:
Code:formulaFullModel <- ~ sample + exon + phenotype:exon + mutatedGene:exon formulaReducedModel <- ~ sample + exon + phenotype:exon
Comment
-
Thinking about it again, I am not sure what will be the output of the testing. To add the "phenotype" as a blocking factor we would also need samples with phenotype "control" and mutatedGene "yes". Of course this is not possible in this experiment! How do the "plotDEXSeq" plots look like?
Maybe would be good to try to use the "formulaFullModel" to estimate dispersions, and the testing doing ignoring the "phenotype" variable...
Let me know what comes out!
Comment
-
Hi Jose,
The default formula would test for differences in exon usage in any of the conditions.
If you are interested in the specific changes between H24-CTRL and between H48-CTRL it is better to subset the data accordingly and run the analysis separately.
Alejandro
Comment
-
Thanks Alejandro,
I thought something similar to to what Mike uses in DESeq2 could be used in DEXSeq but with exons behind is harder.
I'll subset the two and do it separately.
Jose
Originally posted by areyes View PostHi Jose,
The default formula would test for differences in exon usage in any of the conditions.
If you are interested in the specific changes between H24-CTRL and between H48-CTRL it is better to subset the data accordingly and run the analysis separately.
Alejandro
Comment
Latest Articles
Collapse
-
by seqadmin
Proteins are often described as the workhorses of the cell, and identifying their sequences is key to understanding their role in biological processes and disease. Currently, the most common technique used to determine protein sequences is mass spectrometry. While still a valuable tool, mass spectrometry faces several limitations and requires a highly experienced scientist familiar with the equipment to operate it. Additionally, other proteomic methods, like affinity assays, are constrained...-
Channel: Articles
04-04-2024, 04:25 PM -
-
by seqadmin
Despite advancements in sequencing platforms and related sample preparation technologies, certain sample types continue to present significant challenges that can compromise sequencing results. Pedro Echave, Senior Manager of the Global Business Segment at Revvity, explained that the success of a sequencing experiment ultimately depends on the amount and integrity of the nucleic acid template (RNA or DNA) obtained from a sample. “The better the quality of the nucleic acid isolated...-
Channel: Articles
03-22-2024, 06:39 AM -
ad_right_rmr
Collapse
News
Collapse
Topics | Statistics | Last Post | ||
---|---|---|---|---|
Started by seqadmin, 04-11-2024, 12:08 PM
|
0 responses
31 views
0 likes
|
Last Post
by seqadmin
04-11-2024, 12:08 PM
|
||
Started by seqadmin, 04-10-2024, 10:19 PM
|
0 responses
32 views
0 likes
|
Last Post
by seqadmin
04-10-2024, 10:19 PM
|
||
Started by seqadmin, 04-10-2024, 09:21 AM
|
0 responses
28 views
0 likes
|
Last Post
by seqadmin
04-10-2024, 09:21 AM
|
||
Started by seqadmin, 04-04-2024, 09:00 AM
|
0 responses
53 views
0 likes
|
Last Post
by seqadmin
04-04-2024, 09:00 AM
|
Comment