Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • De Novo Transcriptome Assembly from Kmers only

    Background:

    I have RNA-Seq data from 8 time points comparing study to control. Each sample has been processed one sample per lane. I do not trust the reference sequence. What I'm interested in is the most significantly changing transcripts.

    Hypothesis.

    The kmers that are associated with each transcript should change in a coherent manner as the transcript expression changes. Comparing unique kmers first and extracting the most significantly changing kmers should enrich the transcriptome for the genes that are changing the most.

    Problem:
    I have too much data to do a de-novo assembly. It's quite good quality and even eliminating low frequency reads (kmers) I still have too much data to feed to an assembler like Trinity.

    Question.
    Assuming I can select a much smaller set of kmers that are significantly changing. How would I feed the resulting set to an assembler to generate a transcriptome of enriched genes?

    Caveats.
    I don't mind if the contigs that are created from this process results in partial exons associated with the genes that are changing. I can identify them later.

    Soooo

    1) How would you process a set of kmers to feed to an assembler resulting in a fasta file of contigs

    2) If your answer suggests mapping the kmer back to the source read -- can you also suggest how you would do that efficiently (realistically in a decent time frame)

    All thoughts are welcome

    Joe Carl

  • #2
    Instead of trying to pre-select kmers you deem potentially interesting (which is just another way to say pre-biasing your result) you should look at using digital normalization of your data. Digital normalization effectively reduces the input data size and thus the required computational requirements by reducing the input of high abundance kmers. Logically, the 1000th copy of the same kmer does not add any new information to an de novo transcript assembly so you can safely remove copies of kmers above a certain threshold without adversely effecting your assembly.

    Trinity has a digital normalization module built in, described here.

    There is also digital normalization functionality in Titus Brown's khmer suite. Some links:


    I'm out at a Cloud Computing for the Human Microbiome Workshop and I've been trying to convince people of the importance of digital normalization....

    Comment


    • #3
      I have to second kmcarr’s suggest of Trinity’s own digital read normalization process. It took my 4 lanes worth of reads and cut it in about 1/10th the size. Now, there were some genes I picked up when assembling each sample individually that I didn’t find in the pooled and digitally down sampled assembly, so it has some caveats and I would suggest setting the max read depth as high as you machine will allow.

      However, here are some other ideas.

      1) Did you do any quality filtering/trimming? What about adapter and other contaminates clipping? That may help some too.

      2) How much RAM does your machine have? Might it be easier to get access to something with more (i.e. blacklight?).

      3) Have you thought about Trans-ABySS? After some quality filtering and dropping low occurring kmers, you should be able to assemble it on a reasonable number of nodes if you have access to a cluster.

      4) What is your experimental set up? If its a time course, maybe you can just assemble first, middle and end points and skip some samples in-between? There is no reason to have to include all your samples if you think you have enough data from a few to get a representative assembly. You could even do several of these partial assemblies, with different samples, and check to see if you’re finding the same genes. My guess is that in terms of the quality of your assembly you’ll hit a saturation point where more samples doesn’t give better results, just longer run times and higher RAM usage.

      Comment


      • #4
        It's too bad SGA is still in experimental phase, and it appears that it is making genome type assumptions about coverage, otherwise for read correcting and data reduction it is an awesome assembler. Maybe you could just use the read correcting and reduction feature.

        The problem with Kmer reductions, is precisely that if there are lots of data to begin with you get a situation where there are more errors than signal, since the signal is more likely to be redundant therefore discarded, especially in high coverage areas.

        Comment

        Latest Articles

        Collapse

        • seqadmin
          Investigating the Gut Microbiome Through Diet and Spatial Biology
          by seqadmin




          The human gut contains trillions of microorganisms that impact digestion, immune functions, and overall health1. Despite major breakthroughs, we’re only beginning to understand the full extent of the microbiome’s influence on health and disease. Advances in next-generation sequencing and spatial biology have opened new windows into this complex environment, yet many questions remain. This article highlights two recent studies exploring how diet influences microbial...
          02-24-2025, 06:31 AM
        • seqadmin
          Quality Control Essentials for Next-Generation Sequencing Workflows
          by seqadmin




          Like all molecular biology applications, next-generation sequencing (NGS) workflows require diligent quality control (QC) measures to ensure accurate and reproducible results. Proper QC begins at nucleic acid extraction and continues all the way through to data analysis. This article outlines the key QC steps in an NGS workflow, along with the commonly used tools and techniques.

          Nucleic Acid Quality Control
          Preparing for NGS starts with isolating the...
          02-10-2025, 01:58 PM

        ad_right_rmr

        Collapse

        News

        Collapse

        Topics Statistics Last Post
        Started by seqadmin, 03-03-2025, 01:15 PM
        0 responses
        32 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 02-28-2025, 12:58 PM
        0 responses
        138 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 02-24-2025, 02:48 PM
        0 responses
        498 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 02-21-2025, 02:46 PM
        0 responses
        245 views
        0 likes
        Last Post seqadmin  
        Working...
        X