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
                            Strategies for Sequencing Challenging Samples
                            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...
                            03-22-2024, 06:39 AM
                          • seqadmin
                            Techniques and Challenges in Conservation Genomics
                            by seqadmin



                            The field of conservation genomics centers on applying genomics technologies in support of conservation efforts and the preservation of biodiversity. This article features interviews with two researchers who showcase their innovative work and highlight the current state and future of conservation genomics.

                            Avian Conservation
                            Matthew DeSaix, a recent doctoral graduate from Kristen Ruegg’s lab at The University of Colorado, shared that most of his research...
                            03-08-2024, 10:41 AM

                          ad_right_rmr

                          Collapse

                          News

                          Collapse

                          Topics Statistics Last Post
                          Started by seqadmin, Yesterday, 06:37 PM
                          0 responses
                          7 views
                          0 likes
                          Last Post seqadmin  
                          Started by seqadmin, Yesterday, 06:07 PM
                          0 responses
                          7 views
                          0 likes
                          Last Post seqadmin  
                          Started by seqadmin, 03-22-2024, 10:03 AM
                          0 responses
                          49 views
                          0 likes
                          Last Post seqadmin  
                          Started by seqadmin, 03-21-2024, 07:32 AM
                          0 responses
                          66 views
                          0 likes
                          Last Post seqadmin  
                          Working...
                          X