Sim3C: Simulation Of HiC And Meta3C Proximity Ligation Sequencing Technologies

BackgroundChromosome conformation capture (3C) and HiC DNA sequencing methods have rapidly advanced our understanding of the spatial organization of genomes and metagenomes. Many variants of these protocols have been developed, each with their own strengths. Currently there is no systematic means for simulating sequence data from this family of sequencing protocols.\n\nFindingsWe describe a computational simulator that, given reference genome sequences and some basic parameters, will simulate HiC sequencing on those sequences. The simulator models the basic spatial structure in genomes that is commonly observed in HiC and 3C datasets, including the distance-decay relationship in proximity ligation, differences in the frequency of interaction within and across chromosomes, and the structure imposed by cells. A means to model the 3D structure of topologically associating domains (TADs) is provided. The simulator also models several sources of error common to 3C and HiC library preparation and sequencing methods, including spurious proximity ligation events and sequencing error.\n\nConclusionsWe have introduced the first comprehensive simulator for 3C and HiC sequencing protocols. We expect the simulator to have use in testing of HiC data analysis algorithms, as well as more general value for experimental design, where questions such as the required depth of sequencing, enzyme choice, and other decisions must be made in advance in order to ensure adequate statistical power to test the relevant hypotheses.