Cobind: quantitative analysis of the genomic overlaps

Analyzing the overlap between two sets of genomic intervals is a frequent task in the field of bioinformatics. Typically, this is accomplished by counting the number (or proportion) of overlapped regions, which applies an arbitrary threshold to determine if two genomic intervals are overlapped. By making binary calls but disregarding the magnitude of the overlap, such an approach often leads to biased, non-reproducible, and incomparable results.

We developed the cobind package, which incorporates six statistical measures: the Jaccard coefficient, Sørensen-Dice coefficient, Szymkiewicz-Simpson coefficient, collocation coefficient, pointwise mutual information (PMI), and normalized PMI. These measures allow for a quantitative assessment of the collocation strength between two sets of genomic intervals. To demonstrate the effectiveness of these methods, we applied them to analyze CTCF's binding sites identified from ChIP-seq, cancer-specific open-chromatin regions (OCRs) identified from ATAC-seq of 17 cancer types, and oligodendrocytes-specific OCRs identified from scATAC-seq. Our results indicated that these new approaches effectively re-discover CTCF's cofactors, as well as cancer-specific and oligodendrocytes-specific master regulators implicated in disease and cell type development.

The cobind package is implemented in Python and freely available at https://cobind.readthedocs.io/en/latest/.