Question: Advice regarding small RNA annotation to a newly sequenced genome
0
gravatar for Konstantinos Yeles
3 months ago by
Italy
Konstantinos Yeles110 wrote:

Hello Biostars!
Background: I'm working on finding small RNA (piRNA more precisely) sequences on a recently sequenced genome. The genome is in the form of scaffolds currently. I have 3 biological replicates from small RNA seq, and I would like to predict piRNA sequences in these samples. I know that there are some prediction tools for that (either for individual sequences or piRNA clusters) but I would like to see if there is another more "robust" way (Bioconductor packages, or a workflow) or some filtering steps that I haven't used yet.

Questions:

  1. Is it better to collapse all the samples together or use them individually and then search for identical sequences between them?

  2. To predict piRNAs, you need to predict sequences that are tRNAs, rRNAs, miRNAs (this information has to be generated) first and then look on the other sequences for prediction?

My workflow until now: Aligned reads are collapsed and filtered for counts > 9, length of 26 - 34 with Uridine on the 1st base or Adenine on the 10th base (piRNA characteristic). With regard to this article piRNAs are not conserved among species so it would be probably useless to "blast" the sequences to other known piRNAs.

I'm not sure how to proceed.
I used the function plyranges::join_overlap_self_directed(minoverlap = 20) in order to see if there are reads with more or less the same sequence. From that, I got ~ 1.5 million overlaps (starting Genomic ranges(GR) are ~300k) making it more confusing than before. I also tried to "reduce" the GRs (using the function plyranges::reduce_ranges_directed()) but it gave me some ranges of 2kbs long...

As I don't have any experience with annotation, am looking for advice, workflows, or suggestions on the current work.

Thank you for your time,
Konstantinos

ADD COMMENTlink modified 12 weeks ago • written 3 months ago by Konstantinos Yeles110
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1674 users visited in the last hour
_