Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Finding diferentially expressed genes

    Hi All,

    I have the reads count for the following RNA-Seq samples:
    (x1.1, x1.2, x2.1, x2.2, x3.1, x3.2, y1.1, y1.2, y2.1, y2.2, z1.1, z1.2, z2.1, z2.2)

    What i need is to find differentially expressed genes in a group of the samples above, example: differentially expressed genes in all x samples but not in the y nor z.

    Not: x1.1 & x1.2 are biological replicates.

    Please let me know of any way/tool to do this job.

    Thank you!

  • #2
    DESeq
    EdgeR
    DEGseq
    BaySeq


    Theres a few for you to check out.

    Comment


    • #3
      Thanks,
      I tried the a couple of them before, but only to find differntiallly expressed genes between only two samples (x1.1, x1.2, y1.1, y1.2)
      Now the case is a bit more complicated, i want it as i described above to find those DEG is a group of samples that is not expressed in the remaining. If any of these packages provide this can you please let me know which one and if you have links to any examples that looks like mine I would really appreciate that.

      Comment


      • #4
        As far as I know, nobody has yet released (I know of groups who are working on such approaches) a package to deal with data in that fashion. Although if anyone here knows of one, I would also be interested in knowing.

        Otherwise, as far as I know, all approaches require a pairwise comparison between different treatments.

        So the approach I would use would be to do the pairwise comparisons and then filter accordingly.

        So find your differentially expressed genes between x-y, x-z, y-z then find genes that are differentially expressed in both x-y and x-z and filter out genes that are differentially expressed in y-z. You should be left with a genes that are only differentially expressed between x and both y and z conditions, but which are unchanged in between your y and z conditions.

        We actually used a similar approach to this when comparing transcriptomes between 3 different varieties of the same species, but looking for differences in only one of the conditions that could explain a phenotype.

        Comment


        • #5
          Thanks for the details

          When you say "find your differentially expressed genes between x-y", do you mean finding differentially expressed genes in all x conditions (x1.1, x1.2, x2.1, x2.2, x3.1, x3.2) vs all y conditions (y1.1, y1.2, y2.1, y2.2) ???

          Or do I have first to find out those differentially expressed between x conditions:
          (x1.1,x1.2) vs (x2.1 vs x2.2)
          (x1.1,x1.2) vs (x3.1 vs x3.2)
          (x2.1 vs x2.2) vs (x3.1 vs x3.2)
          and then the same for y conditions
          (y1.1, y1.2) vs (y2.1, y2.2).

          Up to my knowledge i can only use the packages above with only two conditions at a time
          (x1.1,x1.2) vs (x2.1 vs x2.2)

          I am not sure if there is a way to do two samples with more than one condition at the same time, as follow,
          (x1.1, x1.2, x2.1, x2.2, x3.1, x3.2) vs (y1.1, y1.2, y2.1, y2.2) ???

          Or that is considered same as my first question!!

          Thanks

          Comment


          • #6
            Originally posted by Amative View Post
            Thanks for the details

            When you say "find your differentially expressed genes between x-y", do you mean finding differentially expressed genes in all x conditions (x1.1, x1.2, x2.1, x2.2, x3.1, x3.2) vs all y conditions (y1.1, y1.2, y2.1, y2.2) ???

            Or do I have first to find out those differentially expressed between x conditions:
            (x1.1,x1.2) vs (x2.1 vs x2.2)
            (x1.1,x1.2) vs (x3.1 vs x3.2)
            (x2.1 vs x2.2) vs (x3.1 vs x3.2)
            and then the same for y conditions
            (y1.1, y1.2) vs (y2.1, y2.2).

            Up to my knowledge i can only use the packages above with only two conditions at a time
            (x1.1,x1.2) vs (x2.1 vs x2.2)

            I am not sure if there is a way to do two samples with more than one condition at the same time, as follow,
            (x1.1, x1.2, x2.1, x2.2, x3.1, x3.2) vs (y1.1, y1.2, y2.1, y2.2) ???

            Or that is considered same as my first question!!

            Thanks
            According to your opening post, all your x samples are biological reps, all your y samples are biological reps, and dido for all your z samples.

            You NEED biological reps to do differential expression. So what I mean by doing pairwise comparisons of the conditions is all biological reps of x versus all biological reps of y and so on. There is no point in separating these replicates as they are what allows you to estimate the biological variance between samples.

            Comment


            • #7
              I am sorry, I was not clear enough before,

              All x samples are from the same tissue but different stages
              (x1.1,x1.2) are biological replicates of tissue x stage 1
              (x2.1,x2.2) are biological replicates of tissue x stage 2
              .
              .

              (y1.1,y1.2) are biological replicates of tissue y stage 1
              (y2.1,y2.2) are biological replicates of tissue y stage 2
              .
              .
              .

              I hope it is more clear now!

              Comment


              • #8
                Ah, this vastly complicates things. Can you clarify then exactly what your aim is. Since X1 is different from X2 are you interested only in a specific stage or a specific condition regardless of stage?

                Comment


                • #9
                  Yea, I know it is complicated !!

                  Actually, I need both,
                  Find those differentially expressed between any two conditions
                  (x1.1,x1.2) vs (y1.1,y1.2)
                  which can be easily done with the packages listed above.

                  The one I am still confused how to do it, is finding those genes differentially expressed genes in all x conditions but not in any of y conditions.

                  1) I'm thinking of finding the DEG between each possible pair of conditions from sample x
                  then
                  2) I'm thinking of finding the DEG between each possible pair of conditions from sample y

                  then
                  3) Union all results from 1
                  4) Union all results from 2

                  find out those genes in 3 and not 4


                  what do you think?!! make sense?

                  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
                  12 views
                  0 likes
                  Last Post seqadmin  
                  Started by seqadmin, Yesterday, 06:07 PM
                  0 responses
                  10 views
                  0 likes
                  Last Post seqadmin  
                  Started by seqadmin, 03-22-2024, 10:03 AM
                  0 responses
                  51 views
                  0 likes
                  Last Post seqadmin  
                  Started by seqadmin, 03-21-2024, 07:32 AM
                  0 responses
                  68 views
                  0 likes
                  Last Post seqadmin  
                  Working...
                  X