Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Help!!!!!! RNA-Seq

    Hi all,

    We did RNA-Seq using RNA isolated from six E11.5 mouse embryonic hearts. And we put all six samples in one lane for Illunima Hiseq 2000 RNA-Seq. Each sample contains around 800 ng total RNA. rRNA was depleted before RNA sequencing.

    As far as I know, one lane can produce around 150 million reads. My aim is to see differentially expressed isoforms.

    I am worrying if the depth coverage for each of my samples is enough.
    Did I make a wrong decision to put all six samples in one lane? The reason for this decision is that we thought we already deleted rRNA which normally eat up most of the reads. The mRNA left in the sample won't need so much reads. This is why we put all six samples in one lane.

    Could you please let me know your opinion? Do we still have enough depth coverage for isoforms?

    Thanks

  • #2
    At 25M per sample, my guess is that you're on the low side of the number of reads you want. It also depends on how you 'depleted' your sample of rRNA. The most 'aggressive' way of doing this is with some sort of poly dT pulldown to isolate the polyA fraction. With that method, 25M would be ok, but you'll likely only see the medium to highly expressed isoforms. If you used RiboMinus, you may find yourself very short of the ideal # of reads, both because some of your reads will be dedicated to ncRNAs and also because RiboMinus tends not to deplete all of the rRNA. (I've heard that Ribo-Zero is better, but I don't know how much better). My guess is that you might want 2-3X more reads per sample (or more if you want the rare isoforms). Curious to hear what others think.

    Comment


    • #3
      Originally posted by guzhi100 View Post
      Hi all,

      We did RNA-Seq using RNA isolated from six E11.5 mouse embryonic hearts. And we put all six samples in one lane for Illunima Hiseq 2000 RNA-Seq. Each sample contains around 800 ng total RNA. rRNA was depleted before RNA sequencing.

      As far as I know, one lane can produce around 150 million reads. My aim is to see differentially expressed isoforms.

      I am worrying if the depth coverage for each of my samples is enough.
      Did I make a wrong decision to put all six samples in one lane? The reason for this decision is that we thought we already deleted rRNA which normally eat up most of the reads. The mRNA left in the sample won't need so much reads. This is why we put all six samples in one lane.

      Could you please let me know your opinion? Do we still have enough depth coverage for isoforms?

      Thanks
      Multiplexing all 6 samples is the right approach. It is the optimal design to run all your replicates in the same lane to minimize effects across lanes.

      ~25 million reads per sample, you will be able to detect differential expression in your more highly expressed genes. But for more lowly expressed genes, your depth will likely be too low.

      However, since you multiplexed your samples (the right thing to do), if you find you do not have not sequenced enough, then the solution is simple, just run all six samples on another lane to get the extra depth you need.

      Comment


      • #4
        Hi, thanks so much for your answers, which are really helpful.

        Comment


        • #5
          I have almost the exact same concern. Thanks very much, the comments have been very useful. I analyzed my 3 samples results using cuffdiff, and what I was wondering was, what is the minimum number of mapped fragments for a gene to still make statements about its expression?

          Comment


          • #6
            bump?

            Sorry to be pest, if this isn't the most interesting question, but I really would love some more info on it. I had a run with my 3 samples on one lane of HiSeq1000. I'm about to run my duplicates, and my lab can only afford one more lane. What I'm wondering is, should I repeat the first run and put three samples in one lane? Or, should I put 6 samples (2 sets of the 3 conditions) on one lane?

            I'm trying to balance the loss of coverage depth versus duplicates. I'm using human epithelial cells.

            To make this determination, I thought I would need to look at my data. I have about 400 "significantly" differentially expressed genes. I have about 100 M mapped hits per condition. So for a FPKM of 1, I guess I would have 100 counts for that gene. If I ran twice as many samples, I guess my FPKM would be 0.5, so with a total count of 50, could I still make determination on gene expression?

            Comment


            • #7
              If all you are looking to do is get differentially expressed genes, you are better off sequencing more samples to help estimate biological variability.

              Comment


              • #8
                Thanks very much!

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