FlaHMM: unistrandflamenco-like piRNA cluster prediction inDrosophilaspecies using hidden Markov models

Abstract

AbstractPIWI-interacting RNAs are a class of small non-coding RNAs that are essential for transposon control in animal gonads. InDrosophilaovarian somatic cells, piRNAs are transcribed from large genomic regions called piRNA clusters, which are enriched for transposon fragments and acts as a memory of past invasions. Despite being widely present acrossDrosophilaspecies, somatic piRNA clusters are notoriously difficult to identify and study due to their lack of sequence conservation and limited synteny. Current identification methods rely either on extensive manual curation or availability of high-throughput small RNA-seq data, limiting large-scale comparative studies. We now present FlaHMM, a hidden Markov model developed to automate genomic annotation offlamenco-like unistrand piRNA clusters inDrosophilaspecies without the need of experimental data beyond a genome assembly. FlaHMM uses transposable element content across 5 or 10 kb bins calculated from genome sequence alone and is thus able to detect candidate piRNA clusters without the need to obtain flies and experimentally perform small RNA sequencing. We show that FlaHMM performs on par with piRNA-guided or manual methods, and thus provides a scalable and efficient approach to piRNA cluster annotation in new genome assemblies. FlaHMM is freely available athttps://github.com/Hannon-lab/FlaHMMunder an MIT licence.