Seqanswers Leaderboard Ad

Collapse
X
 
  • Filter
  • Time
  • Show
Clear All
new posts
  • zinky
    Member
    • Dec 2011
    • 48

    How to decide the two sample are the same in trancription level?

    Hi, guys.
    Is there any criterion to help us to draw a conclusion that two samples are the same or unchanged in trancription level based on RNA-seq data? Correlation test ? or No differential gene indentified ?
  • fanli
    Senior Member
    • Jul 2014
    • 197

    #2
    Do you mean absolute or relative transcription levels? If absolute, then you'd need an external spike-in like the ERCCs. Relative, I'd say the lack of differentially expressed genes as even in a normal two-condition comparison the correlation values tend to be very high (>0.95).

    Comment

    • zinky
      Member
      • Dec 2011
      • 48

      #3
      Originally posted by fanli View Post
      Do you mean absolute or relative transcription levels? If absolute, then you'd need an external spike-in like the ERCCs. Relative, I'd say the lack of differentially expressed genes as even in a normal two-condition comparison the correlation values tend to be very high (>0.95).
      thanks for your suggestion. yes i mean relative transcripts level, because the library size may result in differ counts number( PS: there are alseo some other variables). After normalization, my data show correlation test << 0.05, and the correlation values is 0.97, but there are still some(300-500) DE genes based on Foldchange (abs(threshold) >2). How to explan this ?
      thanks!

      Comment

      • amitm
        Member
        • Feb 2011
        • 52

        #4
        hi zinky,
        You say based on fold change. Did you have replicates? If yes then maybe you can apply a standard statistical testing like a t-test.
        If you are not expecting DE and seeing some (by fold change) then those might be due to sample to sample variation.
        A t-test would consider standard dev. and then give p-values

        Comment

        • zinky
          Member
          • Dec 2011
          • 48

          #5
          hi amitm, i have no replicates but only two samples here. I know the way to detect DE features while my purpose is to demonstrate that there is no difference between them. You advice is useful for understanding but may not work for making a conclusion.
          thanks

          Comment

          • fanli
            Senior Member
            • Jul 2014
            • 197

            #6
            Are you getting 300-500 DE genes from a pooled variance model (e.g. DESeq or edgeR) with no replicates? Or is that based solely on the fold change?

            I'd be really surprised if you managed to get 300 hits with no replicates even between two very different cell types.

            Comment

            • SNPsaurus
              Registered Vendor
              • May 2013
              • 525

              #7
              Originally posted by amitm View Post
              hi zinky,
              You say based on fold change. Did you have replicates? If yes then maybe you can apply a standard statistical testing like a t-test.
              If you are not expecting DE and seeing some (by fold change) then those might be due to sample to sample variation.
              A t-test would consider standard dev. and then give p-values
              Care should be taken not to interpret a t-test p-value that is above some threshold like 0.05 to mean that the samples are the same. The samples could, for instance, be very different but the experiment lacks the statistical power to detect the differences, or the variance is too high (as you say). But the opposite of differently expressed with statistical significance is not similarly expressed with statistical significance, it is "not differently expressed with statistical significance". There are tests that try to measure statistically if two distributions are the same, and this was the basis of a flurry of reports on problems with how statistics are used in many clinical studies, but the exact tests escape me!
              Providing nextRAD genotyping and PacBio sequencing services. http://snpsaurus.com

              Comment

              • zinky
                Member
                • Dec 2011
                • 48

                #8
                Originally posted by fanli View Post
                Are you getting 300-500 DE genes from a pooled variance model (e.g. DESeq or edgeR) with no replicates? Or is that based solely on the fold change?

                I'd be really surprised if you managed to get 300 hits with no replicates even between two very different cell types.
                yes, FoldChange is taken in consideration mainly. I also performed a test which was proposed by Audic S et al in 1997(ref:Audic S, Claverie JM. The significance of digital gene expression profiles. Genome Res. 1997;7:986–995. ). This statistic
                model combine the bayes and poisson distribution, it finally gives a p value to help screening DE features(fdr was also calculated). So based on this method, i got 300-500 DE genes while the correlation of two samples transcript profile is bigger than 0.97.

                Comment

                Latest Articles

                Collapse

                • seqadmin
                  Pathogen Surveillance with Advanced Genomic Tools
                  by seqadmin




                  The COVID-19 pandemic highlighted the need for proactive pathogen surveillance systems. As ongoing threats like avian influenza and newly emerging infections continue to pose risks, researchers are working to improve how quickly and accurately pathogens can be identified and tracked. In a recent SEQanswers webinar, two experts discussed how next-generation sequencing (NGS) and machine learning are shaping efforts to monitor viral variation and trace the origins of infectious...
                  03-24-2025, 11:48 AM
                • seqadmin
                  New Genomics Tools and Methods Shared at AGBT 2025
                  by seqadmin


                  This year’s Advances in Genome Biology and Technology (AGBT) General Meeting commemorated the 25th anniversary of the event at its original venue on Marco Island, Florida. While this year’s event didn’t include high-profile musical performances, the industry announcements and cutting-edge research still drew the attention of leading scientists.

                  The Headliner
                  The biggest announcement was Roche stepping back into the sequencing platform market. In the years since...
                  03-03-2025, 01:39 PM

                ad_right_rmr

                Collapse

                News

                Collapse

                Topics Statistics Last Post
                Started by seqadmin, 03-20-2025, 05:03 AM
                0 responses
                49 views
                0 reactions
                Last Post seqadmin  
                Started by seqadmin, 03-19-2025, 07:27 AM
                0 responses
                57 views
                0 reactions
                Last Post seqadmin  
                Started by seqadmin, 03-18-2025, 12:50 PM
                0 responses
                50 views
                0 reactions
                Last Post seqadmin  
                Started by seqadmin, 03-03-2025, 01:15 PM
                0 responses
                201 views
                0 reactions
                Last Post seqadmin  
                Working...