Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • RNA-seq analysis- finding isomiRs in BAM files?

    I am dealing with seq data for the first time and currently I am analyzing some Illumina RNA-seq libraries for identifying isoMiRs.
    I have processed the raw data using different tools (FASTX,FASTQC,GALAXY) and I am stuck at the mapping step.
    I have aligned my reads to the pri-miRNA human sequences (from MirBase) using Bowtie and I have converted the obtained sam file to a bam file which I've planned to analyze using RSamtools...However I am not sure if that is the best option...


    What is the "standard" procedure for analyzing the alignment files for dentifying isoMiRs?

    Thanks.

  • #2
    Hi,

    I am not sure about the "standard" procedure bur there are few steps you should take before any analysis of micorRNA
    1. You mentioned that you are using RNA-seq library is it small RNA-seq or the total transcriptome? In the second case you can't analyze microRNAs from this data.
    2. How long are your reads? The standard read length of Illumina is 35-38nt or 50nt. As you might notice none of these lengths resemble average microRNA length. In this case you have to trim the adapters that were ligated to the mature microRNAs and sequenced with them as well.

    Regarding the mapping : If you are using Bowtie you have to carefully examine the parameters that you use, take a look at the bowtie manual and decide for yourself how many mismatches you want to allow (since you are looking for Isomirs you might want to be somewhat flexible with this parameter) other parameters that you might want to carefully use and adjust -m, --best,--strata.

    If I remember correctly there are no pri-mir sequences in mirBase but pre-mirs (hairpin loop).

    I am not sure what kind of analysis you want to run with Samtools on the data when looking for isomirs, but I do strongly suggest that you use the IGV viewer to look at the *.bam files before running any analyses so that you will get an idea on how your data looks like...

    Good Luck!

    Comment


    • #3
      Hi moriah and thanks for the reply.

      1. It is only small RNAs.

      2.I have trimmed the adapters and discarded the the reads lower than 18 nts and the ones which had low quality scores. Also I have trimmed the reads and now I am using reads with an average length of 22 nts.

      Regarding the mirBase sequences, I am using the pre-mirs from mirBase.

      I will try to have a look at the data using IGV viewer...Otherwise I am not really sure how I can proceed after generating the (assumingly right) alignment file . Do you have any suggestions?

      Comment


      • #4
        I am no expert with isomirs identification but there are some tools the you might want to get familiar with, like BedTools (coverageBed might be useful). Since the isomirs are defined as variants of the mature microRNA there also align differently to the reference with small shifts from the mature miR coordinates of they fail to align to few nucleotides on the reference seq downstream to the mature seq usually with A's or U's (in rare cases upstream to mature mir as well) , coverageBed tool intersects between your reads and a bed file which contains the coordinates that you wish to intersect with (you can use the mature mirna coordinates). For the reads that will deviate from matures you can check for their abundance to claim for an Isomir pattern. Regarding the Isomirs that went through ADAR modifications you can try to use the "pileup" Samtools function which supposed to find SNPs but might be useful in this case although I am not sure if it is the best way.

        There might be a more accurate and perhaps easy pipelines but for the sake of brainstorming I tried to suggest possible directions...

        Good Luck

        Comment


        • #5
          There might be, but I think that at this point this is very valuable information..I haven't found anywhere documented the whole procedure..The rna-seq analysis descriptions usually jump from the alignment straight to the findings... Thank you.

          Comment

          Latest Articles

          Collapse

          • seqadmin
            Latest Developments in Precision Medicine
            by seqadmin



            Technological advances have led to drastic improvements in the field of precision medicine, enabling more personalized approaches to treatment. This article explores four leading groups that are overcoming many of the challenges of genomic profiling and precision medicine through their innovative platforms and technologies.

            Somatic Genomics
            “We have such a tremendous amount of genetic diversity that exists within each of us, and not just between us as individuals,”...
            05-24-2024, 01:16 PM
          • seqadmin
            Recent Advances in Sequencing Analysis Tools
            by seqadmin


            The sequencing world is rapidly changing due to declining costs, enhanced accuracies, and the advent of newer, cutting-edge instruments. Equally important to these developments are improvements in sequencing analysis, a process that converts vast amounts of raw data into a comprehensible and meaningful form. This complex task requires expertise and the right analysis tools. In this article, we highlight the progress and innovation in sequencing analysis by reviewing several of the...
            05-06-2024, 07:48 AM

          ad_right_rmr

          Collapse

          News

          Collapse

          Topics Statistics Last Post
          Started by seqadmin, 05-24-2024, 07:15 AM
          0 responses
          142 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 05-23-2024, 10:28 AM
          0 responses
          155 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 05-23-2024, 07:35 AM
          0 responses
          154 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 05-22-2024, 02:06 PM
          0 responses
          11 views
          0 likes
          Last Post seqadmin  
          Working...
          X