Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • 48 SOLiD Small RNA Samples...help

    Dear all,

    As per the title, I have 48 small RNA samples that I've run as part of a barcoded experiment. I've downloaded the primary data which includes the cfasta and qual files for analysis.

    My interest is miRNAs and I understand this area is still very new and ever changing.

    Just wondering whether anyone has any recommendations for analysis software. I've run the samples through CLC Bio but having got raw miRNA counts, I am unsure how to procede with normalising these values (my various libraries have very different total counts in some cases as they come from different tissue pathologies).
    Moreover, CLC doesn't give me the aligned reads in a format I can easily manage downstream.

    I'm therefore looking for a complete solution from csfasta to normalised vaules...

    Any help would be much appreciated.

    D

  • #2
    Here's a rough workflow for RNA transcripts using the free/open source formal attire programs. I'm not quite sure how successful this would be for miRNAs, given that cufflinks is designed for transcripts and finding relative isoform abundancies:
    1. split reads based on index barcode somehow (I assume there's SOLiD software that can do this)
    2. create/retrieve genome or transcriptome index in colour-space (e.g. using bowtie-build)
    3. map reads in colour-space using bowtie (-C, --sam options) or tophat (-C option, default output is sam/bam). Do this once for each sample to produce 48 sam files.
    4. run cufflinks on all sam files together to generate a consensus gtf file
    5. run cuffdiff on all sam files together with the gtf file created in the previous step

    Cuffdiff will produce normalised read counts and probabilities, but because it's designed with the assumption that the input data is from RNA transcripts, it may not produce correct results for miRNAs.

    Comment


    • #3
      CLC can export to SAM format. This is necessary for downstream processing with other programs.

      You can obviously normalise for read count, i.e. to counts per million reads or similar. I'm not sure about any miRNA specific read normalisation.

      There is still some debate in the RNA-seq community about normalising for gene length, or in this case for miRNA length.

      Comment


      • #4
        I haven't tried it but I've run across it from time to time: miRanalyzer
        Can be downloaded or run from a webserver. I think it uses DESeq for differential expression.

        miRanalyzer: a microRNA detection and analysis tool for next-generation sequencing experiments. Nucl. Acids Res. (2009) 37 (suppl 2): W68-W76.

        From the site:
        miRanalyzer is a free web-server tool for processing small-RNA data obtained with next generation sequencing platforms such Illumina or SOLiD. The tool requires unique reads in read-count format (i.e., a list of sequences together with the number of times each has been sequenced in the experiment) which can be sequence space (Illumina) or color space (SOLiD). The main features of miRananlyzer are:

        * Mapping all reads against:
        o libraries of known mature microRNAs (including the mature-star libraries – the sequences which pair with the mature microRNAs in the secondary structure of the pre-microRNAs).
        o libraries of theoretically possible mature-star microRNAs which are currently not annotated in miRBase, i.e. have not been observed before.
        o other libraries of transcribed sequences, such as transcriptome and RFam to discard messenger and small non-coding RNAs.
        * Prediction of previously unknown microRNAs by data mining with models specifically trained for plant or animal microRNAs. Detection of differentially expressed microRNAs (including newly predicted ones)
        * Identification of putative target genes for the differentially expressed, new microRNAs.

        Comment


        • #5
          I will make a slightly different suggestion. Consult a service provider specializing in NGS analysis who can help you reduce the time taken to analyze this data and get the data normalized using the various options of RPKM, etc.
          This can be especially important with SOLiD where the protocols for analysis may not be as straightforward as Illumina.
          There is a service provider section on seqanswers. Otherwise some notable vendors who do alignment software and other solutions could help you out. Shop around and you will definitely find a good partner to work with.

          Comment


          • #6
            I suggest map to genome and convert them to fasta files, then use miRdeep.
            For normalization, it's different from the RNA-seq or transcriptome analysis, you can use Sequencing depth to normalize these libraries.

            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
            8 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, Yesterday, 06:07 PM
            0 responses
            8 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, 03-22-2024, 10:03 AM
            0 responses
            49 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, 03-21-2024, 07:32 AM
            0 responses
            66 views
            0 likes
            Last Post seqadmin  
            Working...
            X