The MACPET package can be used for general analysis of paired-end (PET) data like ChIA-PET.
MACPET currently implements the following five stages: Linker filtering (stage 0), mapping to
the reference genome (stage 1), PET classification (stage 2), peak-calling (stage 3) and interaction analysis (stage 4).
The interaction analysis stage (4) currently supports only intra-chromosomal data.
All of the MACPET stages can be run at once, or separately.
In stage 0, MACPET identifies the linkers in the fastq files and classifies the
reads as usable, chimeric or ambiguous. Usable reads are considered in the
subsequent stages. In stage 1, MACPET maps the usable reads to the reference
bowtie and produces a paired-end BAM file. This BAM file is further used in
stage 2 to classify the PETs as self-ligated/intra- or inter-chromosomal.
Self-ligated PETs are used in stage 3 for the identification of significant peaks.
In stage 3, MACPET segments the genome into regions and
applies 2D mixture models for identifying candidate peaks using
skewed generalized students-t distributions (SGT). It then uses a local poisson
model for finding significant binding sites. Finally, in stage 4, MACPET uses the peaks from stage 3 as well as the
Intra-chromosomal PETs from stage 2 to running an additive interaction-analysis model.
MACPET is mainly written in C++,
and it supports the BiocParallel package.
MACPETUlt runs the whole analysis at once.
S4 class for Self-ligated PETs.
S4 class for Self-ligated PETs after peak-calling.
S4 class for Inter-chromosomal PETs.
S4 class for Intra-chromosomal PETs.
S4 class for the interactions found in the data.
Method for plotting different objects.
Method for summarizing different objects.
Method for converting Tags to
Method for converting peaks to
Method for exporting peaks in cvs file format.
Method for converting a
GInteractions class of
Self-ligated PETs to object of
Method for subseting the significant interactions of a
Method for finding the shortest path between two peaks given a set of interactions.
Method for converting peaks to narrowPeak (BED) format for use in
interaction analysis using the MANGO/MICC algorithm. Alternatively the user can run stage 4 at
Function for converting two separate BAM files into one paired-end BAM file.
Prints summary of multiple objects.
Sample ChIA-PET data.
First reads from the sample ChIA-PET data.
Second reads from the sample ChIA-PET data.
Ioannis Vardaxis, email@example.com
Vardaxis I, Drabløs F, Rye M and Lindqvist BH (2018). MACPET: Model-based Analysis for ChIA-PET. To be published.
Vardaxis I, Drabløs F, Rye M and Lindqvist BH (2018). MACPET: Complete pipeline for ChIA-PET. To be published.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.