Integrative Analysis of Histone ChIP-seq and RNA-seq Data

The R package epigenomix has been designed to detect differentially transcribed gene isoforms that, in addition, exhibit altered histone modifications at their respective genomic loci. The package provides methods to map histone ChIP-seq profiles to isoforms and estimate their transcript abundances from RNA-seq data. Based on the differences observed between case and control samples in the RNA-seq and ChIP-seq data, a correlation measure is calculated for each isoform. The distribution of this correlation measure is further investigated by a Bayesian mixture model to (i) reveal the relationship between the studied histone modification and transcriptional activity, and (ii) detect specific isoforms with differences in both transcription values and histone modifications. The method is designed for experiments with a few or no replicates, and is superior to separate analyses of both data types in that setting. This unit illustrates the integrative analysis of ChIP-seq and RNA-seq data with epigenomix. © 2016 by John Wiley & Sons, Inc.