52: M6A Peak Calling Accounting for Sequencing Bias Across Regions and Samples

Music City Center 

N6-Methyladenosine (m6A) is the most abundant type of mRNA methylation and is most widely measured by methylated RNA immunoprecipitation sequencing (MeRIP-seq). In MeRIP-seq, an immunoprecipitation (IP) sample and a pairing control (input) sample are sequenced for each biological sample. Methylated regions are identified as peaks showing increased counts in the IP sample versus the input. We report that technical bias in sequencing can vary substantially in the IP and input samples depending on the local sequence context. Current sequencing depth-based normalization does not appropriately account for the varying technical bias along the transcriptome and leads to inaccurate identification of m6A regions. We describe a method to estimate a local size factor that reflects the RNA sequence context and show that peak calling using these region-specific size factors identifies more accurate peak regions.

transcription

RNA methylation

m6A

MeRIP-seq 

Section on Statistics in Genomics and Genetics