Selfish: discovery of differential chromatin interactions via a self-similarity measure.

MOTIVATION: High-throughput conformation capture experiments, such as Hi-C provide genome-wide maps of chromatin interactions, enabling life scientists to investigate the role of the three-dimensional structure of genomes in gene regulation and other essential cellular functions. A fundamental problem in the analysis of Hi-C data is how to compare two contact maps derived from Hi-C experiments. Detecting similarities and differences between contact maps are critical in evaluating the reproducibility of replicate experiments and for identifying differential genomic regions with biological significance. Due to the complexity of chromatin conformations and the presence of technology-driven and sequence-specific biases, the comparative analysis of Hi-C data is analytically and computationally challenging. RESULTS: We present a novel method called Selfish for the comparative analysis of Hi-C data that takes advantage of the structural self-similarity in contact maps. We define a novel self-similarity measure to design algorithms for (i) measuring reproducibility for Hi-C replicate experiments and (ii) finding differential chromatin interactions between two contact maps. Extensive experimental results on simulated and real data show that Selfish is more accurate and robust than state-of-the-art methods. AVAILABILITY AND IMPLEMENTATION: https://github.com/ucrbioinfo/Selfish.

Chromatin structure and var2csa - a tango in regulation of var gene expression in the human malaria parasite, Plasmodium falciparum?

Over the last few decades, novel methods have been developed to study how chromosome positioning within the nucleus may play a role in gene regulation. Adaptation of these methods in the human malaria parasite, Plasmodium falciparum, has recently led to the discovery that the three-dimensional structure of chromatin within the nucleus may be critical in controlling expression of virulence genes (var genes). Recent work has implicated an unusual, highly conserved var gene called var2csa in contributing to coordinated transcriptional switching, however how this gene functions in this capacity is unknown. To further understand how var2csa influences var gene switching, targeted DNA double-strand breaks (DSBs) within the sub-telomeric region of chromosome 12 were used to delete the gene and the surrounding chromosomal region. To characterize the changes in chromatin architecture stemming from this deletion and how these changes could affect var gene expression, we used a combination of RNA-seq, Chip-seq and Hi-C to pinpoint epigenetic and chromatin structural modifications in regions of differential gene expression. We observed a net gain of interactions in sub-telomeric regions and internal var gene regions following var2csa knockout, indicating an increase of tightly controlled heterochromatin structures. Our results suggest that disruption of var2csa results not only in changes in var gene transcriptional regulation but also a significant tightening of heterochromatin clusters thereby disrupting coordinated activation of var genes throughout the genome. Altogether our result confirms a strong link between the var2csa locus, chromatin structure and var gene expression.