Seqanswers Leaderboard Ad

Collapse
X
 
  • Filter
  • Time
  • Show
Clear All
new posts
  • TimK
    Member
    • Jul 2014
    • 13

    Qualitative comparison of different RNA-Seq experiments

    Hi,

    I'm trying to do a qualitative comparison of publicly available RNA-Seq datasets from different bacteria. Given that the growth conditions are comparable is it possible, e.g. based on RPKM, to make a qualitative comparison of expression levels of conserved genes?

    Cheers,
    Tim
  • Richard Finney
    Senior Member
    • Feb 2009
    • 701

    #2
    What do you think the problems might be?

    Comment

    • TimK
      Member
      • Jul 2014
      • 13

      #3
      I'm not sure how to compare them. Usually I'd go with something like edgeR or similar for DE analysis. But in this case I'm not sure if that's suitable because it's different bacteria, different experimental setups etc. Therefore, I guess it's more a "gene A of orgA relatively highly expressed, in organismB it's not" ... but what normalization/approach to us? Does that make my problem more clear?

      Tim

      Comment

      • Richard Finney
        Senior Member
        • Feb 2009
        • 701

        #4
        You are, indeed, going to have the problems you mentioned.
        The data is from various sources, prepared differently, using different methods, on different machines analyzing expression for sequences that , though conserved, are different sequences.

        Comparing them with deep certainty about the results is going to be challenging.

        Back in the day, I remember some guys using all the CEL files for Affymetrix expression chips and comparing expression between the various expressions from various sources in NIH's GEO (gene expression omnibus) database. A lot easier with a single chip type and single manufacturer's recommended procedure.

        This EdegeR or DESEQ type approach are good for getting a measure of expression differences, that's not your problem.

        The noise and biases of the various methods used to create the reads is your problem.

        Perhaps limiting data to the best/most used protocols and sequencing machine might a good approach.

        Comment

        • Asaf
          Member
          • Jul 2014
          • 20

          #5
          Maybe you should compare pairs

          You're going to stumble into many normalization issues here. I'm having trouble comparing experiments from the same lab in two different experimental settings, let alone different labs and bacteria.
          I would try another approach using simple statistics. I assume you have enough replicates so you can pair the experiments, taking the two conditions in each bacteria, growth phase, lab etc. and in each pair determine, for each gene, if it was up or down regulated (you'll have to normalize them, using DESeq normalization or other method). Then, if you have about 15 (or even less, depends on the number of genes you're testing) such pairs you can run the sign test on each gene. You should get p-values with sufficient power even after FDR correction.

          Comment

          • mbblack
            Senior Member
            • Aug 2009
            • 245

            #6
            If it were me, I would not even try to merge the two datasets. Analyze each independently just as you would any other cross-species comparison. Pull out your two independent gene lists and look for the overlap between them. Then you can make some comparisons of the commonly detected differential expression in each. First thing I'd want to see is do the two gene lists even rank order correlate at all. You can then do some stats if you wished to compare the independently derived gene lists.
            Michael Black, Ph.D.
            ScitoVation LLC. RTP, N.C.

            Comment

            • TimK
              Member
              • Jul 2014
              • 13

              #7
              Thanks everybody.
              The "simple stats" solution is what I was thinking of, too. I will have a look into sign test and rank order conservation. Thanks for the help so far.

              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...