BS-RNA: An efficient mapping and annotation tool for RNA bisulfite sequencing data

Cytosine methylation is one of the most important RNA epigenetic modifications. With the development of experimental technology, scientists attach more importance to RNA cytosine methylation and find bisulfite sequencing is an effective experimental method for RNA cytosine methylation study. However, there are only a few tools can directly deal with RNA bisulfite sequencing data efficiently. Herein, we developed a specialized tool BS-RNA, which can analyze cytosine methylation of RNA based on bisulfite sequencing data and support both paired-end and single-end sequencing reads from directional bisulfite libraries. For paired-end reads, simply removing the biased positions from the 5' end may result in "dovetailing" reads, where one or both reads seem to extend past the start of the mate read. BS-RNA could map "dovetailing" reads successfully. The annotation result of BS-RNA is exported in BED (.bed) format, including locations, sequence context types (CG/CHG/CHH, H=A,T, or C), reference sequencing depths, cytosine sequencing depths, and methylation levels of covered cytosine sites on both Watson and Crick strands. BS-RNA is an efficient, specialized and highly automated mapping and annotation tool for RNA bisulfite sequencing data. It performs better than the existing program in terms of accuracy and efficiency. BS-RNA is developed by Perl language and the source code of this tool is freely available from the website: http://bs-rna.big.ac.cn.

[1]  A. Noor,et al.  Mutation in NSUN2, which encodes an RNA methyltransferase, causes autosomal-recessive intellectual disability. , 2012, American journal of human genetics.

[2]  Zlatko Trajanoski,et al.  meRanTK: methylated RNA analysis ToolKit , 2016, Bioinform..

[3]  Steven L Salzberg,et al.  Fast gapped-read alignment with Bowtie 2 , 2012, Nature Methods.

[4]  Schraga Schwartz,et al.  Transcriptome-Wide Mapping of 5-methylcytidine RNA Modifications in Bacteria, Archaea, and Yeast Reveals m5C within Archaeal mRNAs , 2013, PLoS genetics.

[5]  Rakesh David,et al.  Conservation of tRNA and rRNA 5-methylcytosine in the kingdom Plantae , 2015, BMC Plant Biology.

[6]  Gonçalo R. Abecasis,et al.  The Sequence Alignment/Map format and SAMtools , 2009, Bioinform..

[7]  Cole Trapnell,et al.  TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions , 2013, Genome Biology.

[8]  Steven L Salzberg,et al.  HISAT: a fast spliced aligner with low memory requirements , 2015, Nature Methods.

[9]  Dagmar Wieczorek,et al.  Mutations in NSUN2 cause autosomal-recessive intellectual disability. , 2012, American journal of human genetics.

[10]  T. Preiss,et al.  Widespread occurrence of 5-methylcytosine in human coding and non-coding RNA , 2012, Nucleic acids research.

[11]  Jun Liu,et al.  Methylation modifications in eukaryotic messenger RNA. , 2014, Journal of genetics and genomics = Yi chuan xue bao.

[12]  Brian P. Brunk,et al.  Comparative analysis of RNA-Seq alignment algorithms and the RNA-Seq unified mapper (RUM) , 2011, Bioinform..

[13]  Serban Nacu,et al.  Fast and SNP-tolerant detection of complex variants and splicing in short reads , 2010, Bioinform..

[14]  Bradley R. Cairns,et al.  Identification of direct targets and modified bases of RNA cytosine methyltransferases , 2013, Nature Biotechnology.

[15]  Jizhen Li,et al.  5-hydroxymethylcytosine is detected in RNA from mouse brain tissues , 2016, Brain Research.

[16]  Thomas R. Gingeras,et al.  STAR: ultrafast universal RNA-seq aligner , 2013, Bioinform..

[17]  S. Gabriel,et al.  Whole exome sequencing identifies a splicing mutation in NSUN2 as a cause of a Dubowitz-like syndrome , 2012, Journal of Medical Genetics.

[18]  M. Gorospe,et al.  The tRNA methyltransferase NSun2 stabilizes p16INK4 mRNA by methylating the 3′-untranslated region of p16 , 2012, Nature Communications.