PEALS 1.2.3.1

Peak-based Enhancement Analysis PipeLine for MeRIP-seq

PEALS: Peak-based Enhancement Analysis PipeLine for MeRIP-Seq

Latest Release:

• Github:
• PyPI:

Introduction

MeRIP-seq stands for methylated RNA immunoprecipitation sequencing, which is the most popular method for detection of post-transcriptional RNA modifications. Here, we presented the P-based Enhancement Analysis PipeLine for MeRIP-Seq (PEALS), for identifying enriched methylation regions. PEALS is designed for capturing the enriched methylated regions from MeRIP-seq, where IP signals are significantly over the input signals. Unlike other tools, PEALS not only can detect the peaks located across the exons, but also can detect peaks located in introns. It will be very usefull to detect the methylated peaks in caRNAs, primary transcripts, etc. After peak detections, PEALS wrap DESeq2 framework to estimate the enrichment and their significance. When detecting the differentially methylated (DM) peaks, PEALS use complex strategies to normalized the reads count under different conditions to avoid over-normalization.

Required third party software

PEALS relies on the following third party software. Before installing PEALS, please to ensure the following software have been properly installed.

• bedTools (>=2.30.0)
• samtools (>=1.14)
• featureCounts (>=2.0.2)

The following R pacakges should be also installed.

• getopt (>=1.20.3)
• DESeq2 (>=1.32.0)
• ggplot2 (>=3.4.1)
• ashr (>=2.2.54)
• apeglm (>=1.14.0)
• glmGamPoi (>=1.4.0)

Install

The common way to install PEALS is through PYPI).

Please check INSTALL document for detail.

Usage

Example for regular peak calling on MeRIP-seq.

For sample matrix, please chek the sample matrix here as reference.

peals callpeak -i <bam directory> -m <sample matrix> -P <prefix> -o <output directory>

Example for differentially methylated peak calling on MeRIP-seq.

For sample matrix, please chek the sample matrix here as reference.

peals diffpeak -i <bam directory> -m <sample matrix> -P <prefix> -o <output directory>

Subcommand	Description
callpeak	Main function to call peaks from alignment results.
diffpeak	Main function to call differentially methylated peaks from alignment results.

Contribute

If you have any questions, suggestion/ideas, or just want to have conversions with developers and other users in the community, we recommand you using discussion section or posting to our Issues page.

Ackowledgement

PEALS project also relies on the following python packages.

• CSAPS
• findpeaks