Changes in deepTools2.0¶
The major changes encompass features for increased efficiency, new sequencing data types, and additional plots, particularly for QC.
Moreover, deepTools modules can now be used by other python programs. The deepTools API example is now part of the documentation.
- correlation and comparisons can now be calculated for bigWig files (in addition to BAM files) using
- RNA-seq: split-reads are now natively supported
- MNase-seq: using the new option
bamCoverage, one can now compute read coverage only taking the 2 central base pairs of each mapped fragment into account.
- All modules have comprehensive and automatic tests that evaluate proper functioning after any modification of the code.
- Virtualization for stability: we now provide a
dockerimage and enable the easy deployment of deepTools via the Galaxy
- Our documentation is now version-aware thanks to readthedocs and
- The API is public and documented.
- We dramatically improved the speed of bigwig related tools (multiBigwigSummary and
computeMatrix) by using the new pyBigWig module.
- It is now possible to generate one composite heatmap and/or meta-gene image based on multiple bigwig files in one go (see computeMatrix, plotHeatmap, and plotProfile for examples)
computeMatrixnow also accepts multiple input BED files. Each is treated as a group within a sample and is plotted independently.
- We added additional filtering options for handling BAM files, decreasing the need for prior filtering using tools other than deepTools: The
--samFlagExcludeparameters can, for example, be used to only include (or exclude) forward reads in an analysis.
- We separated the generation of read count tables from the calculation of pairwise correlations that was previously handled by
bamCorrelate. Now, read counts are calculated first using
multiBigWigCoverageand the resulting output file can be used for calculating and plotting pairwise correlations using
plotCorrelationor for doing a principal component analysis using
Correlation analyses are no longer limited to BAM files – bigwig files are possible, too! (see multiBigwigSummary)
Correlation coefficients can now be computed even if the data contains NaNs.
Added the possibility for hierarchical clustering, besides k-means to
computeMatrixcan now read files with DOS newline characters.
--missingDataAsZerowas renamed to
--skipNonCoveredRegionsfor clarity in
- Read extension was made optional and we removed the need to specify a default fragment length for most of the tools:
--fragmentLengthwas thus replaced by the new optional parameter
- Added option
multiBigwigSummary, which can be used to, for example, skip all ‘random’ chromosomes.
- Added the option for adding titles to QC plots.
- Resolved an error introduced by
numpy version 1.10in
- Improved plotting features for
plotProfilewhen using as plot type: ‘overlapped_lines’ and ‘heatmap’
- Fixed problem with BED intervals in
multiBamSummarythat returned wrongly labeled raw counts.
multiBigwigSummarynow also considers chromosomes as identical when the names between samples differ by ‘chr’ prefix, e.g. chr1 vs. 1.
- Fixed problem with wrongly labeled proper read pairs in a BAM file. We now have additional checks to determine if a read pair is a proper pair: the reads must face each other and are not allowed to be farther apart than 4x the mean fragment length.
bamCompare, the behavior of
scaleFactorwas updated such that now, if given in combination with the normalization options (
--normalizeUsingRPKM), the given scaling factor will be multiplied with the factor computed by the respective normalization method.