Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Differentially expressed genes

    Hi,
    I am doing RNA-seq analysis using Tuxedo pipeline. I wonder when do you name it DIFFERENTIALY expressed genes? Does it have to be significant genes out of cuffdiff analysis or any sort of fold changes between two conditions.
    I am asking it because my two conditions show very few significant changes in gene expression. And that is between G2 and S phase cells I just get 16 significant differentialy expressed genes.
    Let me know if you need further information.

    thanks.

  • #2
    Significance in cuffdiff is when the q-value (FDR adjusted p-value) is less than the p-value itself and it is not dependent on the magnitude of the fold change.

    (Source, the cuffdiff manual : http://cufflinks.cbcb.umd.edu/manual.html#gene_exp_diff )

    It's impossible to tell why you have few significant genes without more information about your experimental set up. Too few reads per sample, or too few replicates would be possible explanations.
    Last edited by Bukowski; 02-24-2014, 02:04 PM. Reason: added link

    Comment


    • #3
      Yes, I understood how Cuffdiff defines significant, but my question is if differentially expressed genes are specifically those that are expressed significantly?

      Second, since we were establishing RNA-seq for the first time in our lab, we did only one replicate of each condition. We just wanted to go all the way to the end before doing actual experiment. I recently read that number of replicates have great impact on the accuracy of data analysis. My data has between 10 to 20 million reads per sample. I guess having one replicate is making the problem. What do you think?

      How many replicates do you suggest with how many reads in each?
      Please let me know what other information will make it more clear.

      Thanks.

      Comment


      • #4
        Great article about the replicates/read numbers
        The code used in this paper is provided on: http://home.uchicago.edu/∼jiezhou/replication/. The expression data is deposited in the Gene Expression Omnibus under the accession ID GSE51403.

        Comment


        • #5
          Most qPCR papers use a minimum fold change of 2 and statistical significance before calling anything differentially expressed. Since RNA-seq is basically high throughput qPCR it would make sense to use the same thresholds.

          Comment


          • #6
            How about using logFold change?

            Comment


            • #7
              Thanks people.

              Comment


              • #8
                When we have no replicates.. then i suggest u to try Gfold.
                Gfold > 0 || gfold < 0 means it Upregulated or down regulated.

                For more information check.
                Krishna

                Comment


                • #9
                  Thanks a lot Krish. Looks very interesting, how can I get the software? Apt-get I guess? I am new to Linux environment.

                  Comment


                  • #10
                    You can download at


                    For downlading to linux :

                    wget https://bitbucket.org/feeldead/gfold....V1.1.1.tar.gz

                    Check Readme and installation files in gfold.
                    Last edited by Krish_143; 02-27-2014, 11:35 PM.
                    Krishna

                    Comment


                    • #11
                      Thanks a million Krish. I will give it a try.

                      Comment


                      • #12
                        This Gfold was a game changing program. Thanks for suggesting it Krish.

                        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
                        10 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
                        67 views
                        0 likes
                        Last Post seqadmin  
                        Working...
                        X