a statistical framework to detect differences in chromatin contact density at genomic regions of interest
Deposited scripts and files for calculation of melting scores of long genes in oligodendroglia (OLG) from the somatosensory cortex, dopaminergic neurons (DNs) from the midbrain VTA and pyramidal glutamatergic neurons (PGNs) from the hippocampus CA1 relative to embryonic stem cells (ESCs).
library(tidyverse)
library(ggpubr)- code/MELTRON.R:
Compares insulation score (IS) distributions over long genes in OLGs, DNs, PGNs to ESCs and calculates melting score per gene. - code/plot_ecdf.R:
Plots empirical cumulative density functions (ECDF) for IS values of individual genes. - code/plot_rna_atac_over_meltingScore.R:
Plots expression and chromatin accessibility as a function of the melting score per cell-type. Density plots with median lines indicate population trends. - code/plot_domain_melting_gene_characteristics.R Plots heatmaps and summary violin plots for each of the clusters.
- command_line_apps/matrix_wide_to_long.R: Converts a square matrix into a long matrix for IS calculation. Accepts wildcards for processing of multiple chromosomes.
- command_line_apps/long_matrix_to_IS.R: Calculates insulation scores at multiple distances (default 100kb - 1Mb, steps of 100kb).
- command_line_apps/Meltron.R: Calculation of melting scores over genomic regions of interest
type
Rscript command_line_apps/matrix_wide_to_long.R --help for explanations
Developed and tested with R version 3.6.0 Planting of a Tree.
Developed and maintained by Dominik Szabó 
with intellectual input from Christoph Thieme .
Please get in touch for questions and issues: dominik.szabo at mdc-berlin.de
Please check our manuscript:
Winick-Ng, W., Kukalev, A., Harabula, I., Zea Redondo, L., Szabó, D. et al.:
Cell-type specialization is encoded by specific chromatin topologies