Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • RNA-editing pipeline for miRNAs?

    Hi,

    This is my very first post and I'm also a newbie in performing bioinformatic analysis, so forgive me if something sounds off. I'm planning on editing this post with additional data in case I'm missing something relevant here.

    I'm analyzing SE miRNA-seq data and for this, I used miRDeep2 software and followed by performing DESeq analysis of both known miRNAs and predicted novel miRNAs.

    I'm now interested in identifying if some of these miRNAs could have been edited post transcriptionally (most likely by substitutional A-to-I RNA editing), and distinguishing them from SNPs.

    I found a pipeline (http://genome.cshlp.org/content/22/8/1533.long) which seems to have been used and cited by quite a few publications. The protocol used is explained here step by step: http://www.tau.ac.il/~elieis/miR_editing/

    I'd like to ask if anyone here has ever used it or know a better/up-to-date protocol. The trimming step and the mapping step seem to work fine but then I proceed to step 3:

    Mapping the mismatches against the pre-miRNA sequences:
    The purpose of this step is to move from reads aligned against the genome (the end-point of the previous step) to counts of each of the four possible nucleotides at each position along the pre-miRNA sequence, for all the pre-miRNAs. Performing this transformation will allow us to focus our analysis only on bona fide miRNA and to use, in the following step, binomial statistics in order to detect significant modifications inside them...

    ... The main output of this script is a text file (termed 'main_output.txt') containing the counts of each of the four possible nucleotides at each position along all the pre-miRNA sequences...

    So after this step, I get an output file which I don't really understand but I don't think it looks right (basically, all 0 at the last section of the file and 0s alongside the file for each of the positions of each miRNA). It looks like this:

    Click image for larger version

Name:	Screenshot_2017-09-18_19-59-38.png
Views:	1
Size:	38.3 KB
ID:	309307

    Click image for larger version

Name:	Screenshot_2017-09-18_20-04-47.png
Views:	1
Size:	40.6 KB
ID:	309308

    Click image for larger version

Name:	Screenshot_2017-09-18_17-10-22.png
Views:	1
Size:	61.2 KB
ID:	309306

    I thought maybe the problem had something to do with the formatting of my input data (.fq files) so I even tried with the data they used for their analysis and still get the same results. If someone has ever followed this pipeline, I would gladly accept some guidance, I'm a bit lost here.

    Thank you in advance

Latest Articles

Collapse

  • 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
  • seqadmin
    The Impact of AI in Genomic Medicine
    by seqadmin



    Artificial intelligence (AI) has evolved from a futuristic vision to a mainstream technology, highlighted by the introduction of tools like OpenAI's ChatGPT and Google's Gemini. In recent years, AI has become increasingly integrated into the field of genomics. This integration has enabled new scientific discoveries while simultaneously raising important ethical questions1. Interviews with two researchers at the center of this intersection provide insightful perspectives into...
    02-26-2024, 02:07 PM

ad_right_rmr

Collapse

News

Collapse

Topics Statistics Last Post
Started by seqadmin, 03-14-2024, 06:13 AM
0 responses
32 views
0 likes
Last Post seqadmin  
Started by seqadmin, 03-08-2024, 08:03 AM
0 responses
71 views
0 likes
Last Post seqadmin  
Started by seqadmin, 03-07-2024, 08:13 AM
0 responses
80 views
0 likes
Last Post seqadmin  
Started by seqadmin, 03-06-2024, 09:51 AM
0 responses
68 views
0 likes
Last Post seqadmin  
Working...
X