Unconfigured Ad

Collapse
X
 
  • Time
  • Show
Clear All
new posts
  • arvid
    Senior Member
    • Jul 2011
    • 156

    #16
    Originally posted by pardonliang View Post
    ... I have used susceptible、resistance 、susceptible and resistance assebled transcriptomes as reference to compare the number of paired-maping reads、singled-mapping reads、up expressed transcripts.The result of stored in Attach Files.
    According to your question(the software and statistics used, the amount and type of replicates, the biological system you are working on, and the way the samples were collected).I haven't used statistics、replicateds yet...
    Alright, so don't expect much from your experiment, since you have different strains and no biological replicates... You specify "RPKM>2" in your table, was this log-transformed RPKM values or just raw RPKM? In any case, I really wouldn't recommended using an arbitrary expression-level focused cutoff, and I guess you can't evaluate those thousands of genes anyway. You might want to feed your counts through edgeR or DESeq in R/Bioconductor instead and use the options for experiments without biological replicates (there are good instructions in their respective manuals, or just search SEQanswers). You will not see a lot of differentially expressed transcripts then, but hopefully you can reduce the number of falsely called differentially expressed transcripts. You should afterwards check that the observed differential expression is not due to lower mapping due to polymorphisms in one of the samples, however.
    You could try IsoEM or RSEM as well, though I'm not sure how sensitive their respective EM algorithms are to polymorphisms in the underlying mappings. We used RSEM in a similar case to yours (I wasn't involved until a late stage), and were able to identify some qualitatively DE genes (confirmed by RT-qPCR). I can't recommend it again though - replicates are really necessary to get information of high quality out of such an experiment.

    Comment

    • shoegame2001
      Member
      • Dec 2010
      • 21

      #17
      I apologize if this question has been addressed elsewhere--I did a little bit of looking but could not find a consensus:

      When combining assemblies of different kmer lengths (in my case using Trinity), what is the best method for eliminating redundancy in the resulting combined assembly to make the most comprehensive reference possible? I have tried using cdhit-est, TGICL, and CAP3, but am still observing redundancy when I blast my assembly to itself. For those who have had success with these tools, what options were you using?

      Thanks!

      Comment

      Latest Articles

      Collapse

      • SEQadmin2
        Multiomic Liquid Biopsy: What DNA/RNA Panels Reveal That DNA Alone Can’t
        by SEQadmin2



        Tumors contain complex mixtures of cell populations that evolve in response to treatment and host-specific factors. Tissue biopsies provide only static snapshots of a tumor’s state. But what if we could follow that evolution over time using a minimally invasive technique?

        Liquid biopsies can detect tumor-derived material circulating in the bloodstream, providing a way to study treatment response, detect emerging resistance, and discover new biomarkers. Circulating
        ...
        10-14-2026, 09:55 AM
      • SEQadmin2
        New Genomics Technologies Take Aim at Long-Standing Limits
        by SEQadmin2


        Researchers using sequencing and genomics tools often have to make trade-offs. They can choose between speed or scale, short reads or long-range information, or targeted panels or a view of the whole transcriptome. New technologies that have been released this year are built to address those tough choices.

        We asked six companies the same four questions to learn about their latest products. The new technologies bring a lot to the table, including rethinking sequencing
        ...
        09-28-2026, 10:25 AM

      ad_right_rmr

      Collapse

      News

      Collapse

      Topics Statistics Last Post
      Started by SEQadmin2, Yesterday, 11:13 AM
      0 responses
      16 views
      0 reactions
      Last Post SEQadmin2  
      Started by SEQadmin2, Yesterday, 10:48 AM
      0 responses
      15 views
      0 reactions
      Last Post SEQadmin2  
      Started by SEQadmin2, 10-15-2026, 10:52 AM
      0 responses
      35 views
      0 reactions
      Last Post SEQadmin2  
      Started by SEQadmin2, 10-14-2026, 05:39 AM
      0 responses
      42 views
      0 reactions
      Last Post SEQadmin2  
      Working...