Seqanswers Leaderboard Ad

Collapse
X
 
  • Filter
  • Time
  • Show
Clear All
new posts
  • metheuse
    Member
    • Jan 2013
    • 84

    How to detect alternative splicing when there is no replicate

    I'm working on a pilot RNAseq project which compares paired-end reads between two samples, and there is only a single experiment in each sample.
    I've been searching for tools to handle alternative splicing and found things like DEXseq and SplicingCompass. But both of them seem to require replicates. I'm wondering how I can use them (or any other tools) on my data.
    I understand that without replicates, the results would lack statistic meaning, but since this is just a pilot project, the thing is to find an appropriate pipeline to analyze the data.
    By the way, could I randomly separate the reads into two sets and treat them as pseudo replicates? I know this is definitely against the definition of biological replicates, but is it a way that would worth a try?
    Thanks for any suggestions.
  • syfo
    Just a member
    • Nov 2012
    • 103

    #2
    Originally posted by metheuse View Post
    the thing is to find an appropriate pipeline to analyze the data.
    That's the point: an appropriate pipeline requires replicates.

    Originally posted by metheuse View Post
    could I randomly separate the reads into two sets and treat them as pseudo replicates? I know this is definitely against the definition of biological replicates, but is it a way that would worth a try?
    This is what I would do to set up a pilot pipeline.

    Comment

    • metheuse
      Member
      • Jan 2013
      • 84

      #3
      Thanks! By the way, do you know how to randomly separate the reads into two groups? What I think of is to shuffle the reads and then pick the first half and the last half. Should I separate them equally?
      Last edited by metheuse; 04-16-2013, 05:57 AM.

      Comment

      • Simon Anders
        Senior Member
        • Feb 2010
        • 995

        #4
        The wole point of a pilot experiment is to figure out whether it is worthwhile to do the real experiment, i.e., whether you have a chance to detect interesting biology once you scale things up.

        To address this question, you have to know effect size (how big are the typical differences between control and treatment) and noise strength (how big are meaningless differences that also appear between two control or two treatment samples). Without performing at least one condition in at least duplicates, there is no way of knowing how much noise you are facing and hence how strong effects have to be in order to be meaningful.

        So, I simply don't get it: What is the point of unreplicated pilot experiments, and why do people waste their time on them?

        If your aim is solely to tets your pipeline, simply take some published data. If you split your unreplicates data in two halves, you will only test the unrealistic situation of nearly every gene with more than a few dozen read counts being called significant.

        Comment

        • metheuse
          Member
          • Jan 2013
          • 84

          #5
          Thanks for your reply, Simon.
          I perfectly understand your point.
          This is not my own project. I'm just in a position to help some biologists to deal with bioinformatics stuff. The unreplicate data is what they have right now and they want me to explore the tools that can process the data. I told them there is no way to determine significance without replicates and they absolutely understand it. They will sequence more replicates as the next step.
          Simply, I just want to construct a pipeline that can take the data and output results. But I understand I can't well evaluate the pipeline since the result would lack statistic meaning.

          Comment

          • syfo
            Just a member
            • Nov 2012
            • 103

            #6
            At least you can set up a bioinformatics pipeline in terms of format, data processing etc. I would just start by splitting the fastq files in two. Shuffling, random subsampling etc sounds better but as Simon pointed out the biological variability cannot really be simulated anyway. I understand your approach though, it makes sense to ensure your workflow can actually process the data the way it comes from your lab. Published data from another source might differ in many ways (format, naming system, etc).

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