Differential expression without replicates
1
0
Entering edit mode
@srikanth-manda-srinivas-4590
Last seen 10.2 years ago
HI all, I have paired-end RNAdata from five different tissues of the same organism, I calculated the read counts using HTSeq tools and now I want to see which are the transcripts which are differentially expressed (if any) among the five tissues (I do not have any replicates). I have tried using edgeR and DESeq packages but not getting even a single value to be significant. Can anyone suggest the possibilities?? Regards, Srikanth -- Srinivas M. Srikanth Ph.D. Student Institute of Bioinformatics Discoverer, 7th Floor, International Technology Park, Bangalore, India Mob:+919019114878 [[alternative HTML version deleted]]
• 1.9k views
ADD COMMENT
0
Entering edit mode
Mark Robinson ▴ 880
@mark-robinson-4908
Last seen 6.1 years ago
Hi Srikanth, No replicates is not an ideal situation. In the user manual, edgeR has a section "What to do if you have no replicates": http://www.bioconductor.org/packages/release/bioc/vignettes/edgeR/inst /doc/edgeRUsersGuide.pdf Similarly, DESeq has a section "Working without any replicates": http://www.bioconductor.org/packages/release/bioc/vignettes/DESeq/inst /doc/DESeq.pdf Did you try these? Best regards, Mark On 28.08.2012, at 08:51, Srinivas M. Srikanth wrote: > HI all, > > I have paired-end RNAdata from five different tissues of the same organism, > I calculated the read counts using HTSeq tools and now I want to see which > are the transcripts which are differentially expressed (if any) among the > five tissues (I do not have any replicates). I have tried using edgeR and > DESeq packages but not getting even a single value to be significant. Can > anyone suggest the possibilities?? > > > Regards, > Srikanth > > -- > Srinivas M. Srikanth > Ph.D. Student > Institute of Bioinformatics > Discoverer, 7th Floor, > International Technology Park, > Bangalore, India > Mob:+919019114878 > > [[alternative HTML version deleted]] > > _______________________________________________ > Bioconductor mailing list > Bioconductor at r-project.org > https://stat.ethz.ch/mailman/listinfo/bioconductor > Search the archives: http://news.gmane.org/gmane.science.biology.informatics.conductor ---------- Prof. Dr. Mark Robinson Bioinformatics Institute of Molecular Life Sciences University of Zurich Winterthurerstrasse 190 8057 Zurich Switzerland v: +41 44 635 4848 f: +41 44 635 6898 e: mark.robinson at imls.uzh.ch o: Y11-J-16 w: http://tiny.cc/mrobin ---------- http://www.fgcz.ch/Bioconductor2012 http://www.eccb12.org/t5
ADD COMMENT
0
Entering edit mode
Sriniva you said "I have tried ... but not getting even a single value to be significant." Significance relates to the rejection of a (null) hypothesis in a test. What null hypothesis did you test? - was it appropriate for the biological question you are after? - if so, may be there is just not enough signal? Best wishes Wolfgang Aug/28/12 9:07 AM, Mark Robinson scripsit:: > Hi Srikanth, > > No replicates is not an ideal situation. > > In the user manual, edgeR has a section "What to do if you have no replicates": > http://www.bioconductor.org/packages/release/bioc/vignettes/edgeR/in st/doc/edgeRUsersGuide.pdf > > Similarly, DESeq has a section "Working without any replicates": > http://www.bioconductor.org/packages/release/bioc/vignettes/DESeq/in st/doc/DESeq.pdf > > Did you try these? > > Best regards, > Mark > > > On 28.08.2012, at 08:51, Srinivas M. Srikanth wrote: > >> HI all, >> >> I have paired-end RNAdata from five different tissues of the same organism, >> I calculated the read counts using HTSeq tools and now I want to see which >> are the transcripts which are differentially expressed (if any) among the >> five tissues (I do not have any replicates). I have tried using edgeR and >> DESeq packages but not getting even a single value to be significant. Can >> anyone suggest the possibilities?? >> >> >> Regards, >> Srikanth >> >> -- >> Srinivas M. Srikanth >> Ph.D. Student >> Institute of Bioinformatics >> Discoverer, 7th Floor, >> International Technology Park, >> Bangalore, India >> Mob:+919019114878 >> >> [[alternative HTML version deleted]] >> >> _______________________________________________ >> Bioconductor mailing list >> Bioconductor at r-project.org >> https://stat.ethz.ch/mailman/listinfo/bioconductor >> Search the archives: http://news.gmane.org/gmane.science.biology.informatics.conductor > > ---------- > Prof. Dr. Mark Robinson > Bioinformatics > Institute of Molecular Life Sciences > University of Zurich > Winterthurerstrasse 190 > 8057 Zurich > Switzerland > > v: +41 44 635 4848 > f: +41 44 635 6898 > e: mark.robinson at imls.uzh.ch > o: Y11-J-16 > w: http://tiny.cc/mrobin > -- Best wishes Wolfgang Wolfgang Huber EMBL http://www.embl.de/research/units/genome_biology/huber
ADD REPLY
0
Entering edit mode
Dear Dr. Huber, I want to see that if there is any difference in the expression levels of isoforms in different tissues (five in this case). So, the null here is there is no significant difference in the expression levels across tissues. I am trying to get the top differentials and than try to correlate the data with proteogenomics. Regards, Srikanth On Tue, Aug 28, 2012 at 1:04 PM, Wolfgang Huber <whuber@embl.de> wrote: > Sriniva > > you said "I have tried ... but not getting even a single value to be > significant." Significance relates to the rejection of a (null) hypothesis > in a test. What null hypothesis did you test? > > - was it appropriate for the biological question you are after? > - if so, may be there is just not enough signal? > > Best wishes > Wolfgang > > > Aug/28/12 9:07 AM, Mark Robinson scripsit:: > > Hi Srikanth, >> >> No replicates is not an ideal situation. >> >> In the user manual, edgeR has a section "What to do if you have no >> replicates": >> http://www.bioconductor.org/**packages/release/bioc/** >> vignettes/edgeR/inst/doc/**edgeRUsersGuide.pdf<http: www.bioconduc="" tor.org="" packages="" release="" bioc="" vignettes="" edger="" inst="" doc="" edgerusersguide="" .pdf=""> >> >> Similarly, DESeq has a section "Working without any replicates": >> http://www.bioconductor.org/**packages/release/bioc/** >> vignettes/DESeq/inst/doc/**DESeq.pdf<http: www.bioconductor.org="" pa="" ckages="" release="" bioc="" vignettes="" deseq="" inst="" doc="" deseq.pdf=""> >> >> Did you try these? >> >> Best regards, >> Mark >> >> >> On 28.08.2012, at 08:51, Srinivas M. Srikanth wrote: >> >> HI all, >>> >>> I have paired-end RNAdata from five different tissues of the same >>> organism, >>> I calculated the read counts using HTSeq tools and now I want to see >>> which >>> are the transcripts which are differentially expressed (if any) among the >>> five tissues (I do not have any replicates). I have tried using edgeR and >>> DESeq packages but not getting even a single value to be significant. Can >>> anyone suggest the possibilities?? >>> >>> >>> Regards, >>> Srikanth >>> >>> -- >>> Srinivas M. Srikanth >>> Ph.D. Student >>> Institute of Bioinformatics >>> Discoverer, 7th Floor, >>> International Technology Park, >>> Bangalore, India >>> Mob:+919019114878 >>> >>> [[alternative HTML version deleted]] >>> >>> ______________________________**_________________ >>> Bioconductor mailing list >>> Bioconductor@r-project.org >>> https://stat.ethz.ch/mailman/**listinfo/bioconductor<https: stat.="" ethz.ch="" mailman="" listinfo="" bioconductor=""> >>> Search the archives: http://news.gmane.org/gmane.** >>> science.biology.informatics.**conductor<http: news.gmane.org="" gman="" e.science.biology.informatics.conductor=""> >>> >> >> ---------- >> Prof. Dr. Mark Robinson >> Bioinformatics >> Institute of Molecular Life Sciences >> University of Zurich >> Winterthurerstrasse 190 >> 8057 Zurich >> Switzerland >> >> v: +41 44 635 4848 >> f: +41 44 635 6898 >> e: mark.robinson@imls.uzh.ch >> o: Y11-J-16 >> w: http://tiny.cc/mrobin >> >> > > -- > Best wishes > Wolfgang > > Wolfgang Huber > EMBL > http://www.embl.de/research/**units/genome_biology/huber<http: www.="" embl.de="" research="" units="" genome_biology="" huber=""> > > > ______________________________**_________________ > Bioconductor mailing list > Bioconductor@r-project.org > https://stat.ethz.ch/mailman/**listinfo/bioconductor<https: stat.et="" hz.ch="" mailman="" listinfo="" bioconductor=""> > Search the archives: http://news.gmane.org/gmane.** > science.biology.informatics.**conductor<http: news.gmane.org="" gmane.="" science.biology.informatics.conductor=""> > -- Srinivas M. Srikanth Ph.D. Student Institute of Bioinformatics Discoverer, 7th Floor, International Technology Park, Bangalore, India Mob:+919019114878 [[alternative HTML version deleted]]
ADD REPLY

Login before adding your answer.

Traffic: 857 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6