Seqanswers Leaderboard Ad

Collapse

Announcement

Collapse
No announcement yet.
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • differentiated expression analysis on isoforms/transcripts level

    Hi there,

    I have a RNA-seq data (3 cases and 4 controls). I am looking for ways to identify differentiated isoforms/transcripts in these two conditions. Any suggestions besides t-test?

    Thanks.

  • #2
    DESeq and EdgeR have been frequently recommended in this forum

    Comment


    • #3
      As far as I understand, DESeq and EdgeR identify differentiated genes, but not transcripts/isoforms. I wish they have extended their software's ability to isoform level.

      Comment


      • #4
        Cufflinks is probably your best bet for isoform level analyses.

        Comment


        • #5
          ALEXA-SEQ could also be worth looking into.

          And DEXSeq looks at differential exon usage.

          Comment


          • #6
            There are also IsoEM and, a bit different, DEXSeq

            Comment


            • #7
              "As far as I understand, DESeq and EdgeR identify differentiated genes, but not transcripts/isoforms."

              This doesn't make sense to me. These programs analyze differential expression, whether that occurs at the gene level or at the transcript isoform level shouldn't matter. The mapping program that you use to identify splice junctions will help sort out the various isoforms and some mappers are better than others. Commonly used programs include ERANGE and SpliceMap among others. Bowtie/Tophat is not my favorite because it just doesn't perform quite as well as some other aligners and you end up locked into Cufflinks, which is not necessarily the best option. Whatever you do, don't just perform a t-test! Some useful recent reading includes Risso, et al., doi:10.1186/1471-2105-12-480, or Tarazona, et al., doi:10.1101/gr.124321.111 or several others on this topic over the last 6 months.

              Comment


              • #8
                Thanks for all responses.

                That's right. Cufflink-cuffcompare-cuffdiff gives results of paired samples, but not between two conditions.

                I am reading the paper and the manual of "edgeR", you might right. the authors said their method might apply to transcript level. Then I wonder how to get raw count from RNA-seq fastq files. Also the cpm function in edgeR is missing. Anybody know what it is?

                Thanks again.

                Comment


                • #9
                  I found cpm in their source package.

                  Comment


                  • #10
                    The problem is that it's hard to assign integer read counts to transcripts, and edgeR / DESeq work with integers. All the isoform/transcript quantification models that I know of - most are based on EM, like Cufflinks, although there are also linear equation system models etc. - probabilistically assign reads to multiple isoforms with different non-integer weights. One way to circumvent this is to only consider reads that map to a single isoform only - it may be that e g ALEXA-SEQ and NEUMA work in this way. That would achieve integer counts without double-counting reads.

                    Comment


                    • #11
                      Is it a reason that we can not fit a GLM model with FPKM value directly since FPKM is a probability?

                      Comment


                      • #12
                        found some helpful information here in case one is wondering this type of question like me:http://seqanswers.com/forums/showthread.php?t=4349

                        Comment


                        • #13
                          Originally posted by arrchi View Post
                          That's right. Cufflink-cuffcompare-cuffdiff gives results of paired samples, but not between two conditions.
                          That is wrong. cuffdiff can be used to compute the differential expression between two conditions. For two conditions, you need to give the sample names as comma separated values and differentiate the conditions with a space.
                          Code:
                          cuffdiff --labels treated,untreated cuffcompare_combined.gtf treated_1.sam,treated_2.sam,treated_3.sam untreated_1.sam,untreated_2.sam
                          The cuffcompare_combined.gtf is obtained by running cuffcompare with the above 5 samples together.

                          Comment

                          Latest Articles

                          Collapse

                          • seqadmin
                            Recent Developments in Metagenomics
                            by seqadmin





                            Metagenomics has improved the way researchers study microorganisms across diverse environments. Historically, studying microorganisms relied on culturing them in the lab, a method that limits the investigation of many species since most are unculturable1. Metagenomics overcomes these issues by allowing the study of microorganisms regardless of their ability to be cultured or the environments they inhabit. Over time, the field has evolved, especially with the advent...
                            09-23-2024, 06:35 AM
                          • seqadmin
                            Understanding Genetic Influence on Infectious Disease
                            by seqadmin




                            During the COVID-19 pandemic, scientists observed that while some individuals experienced severe illness when infected with SARS-CoV-2, others were barely affected. These disparities left researchers and clinicians wondering what causes the wide variations in response to viral infections and what role genetics plays.

                            Jean-Laurent Casanova, M.D., Ph.D., Professor at Rockefeller University, is a leading expert in this crossover between genetics and infectious...
                            09-09-2024, 10:59 AM

                          ad_right_rmr

                          Collapse

                          News

                          Collapse

                          Topics Statistics Last Post
                          Started by seqadmin, 10-02-2024, 04:51 AM
                          0 responses
                          8 views
                          0 likes
                          Last Post seqadmin  
                          Started by seqadmin, 10-01-2024, 07:10 AM
                          0 responses
                          13 views
                          0 likes
                          Last Post seqadmin  
                          Started by seqadmin, 09-30-2024, 08:33 AM
                          0 responses
                          18 views
                          0 likes
                          Last Post seqadmin  
                          Started by seqadmin, 09-26-2024, 12:57 PM
                          0 responses
                          16 views
                          0 likes
                          Last Post seqadmin  
                          Working...
                          X