methylation detection assays used methylation modification sensitive restriction endonucleases. genomic dna was digested with both methylation-sensitive and insensitive restriction enzymes recognizing the same restriction site. the idea being that whenever the site was methylated, only the methylation insensitive enzyme could cleave at that position. by comparing restriction fragment sizes generated from the methylation-sensitive enzyme to those of the methylation-insensitive enzyme, it was possible to determine the methylation pattern of the region. this analysis step was done by amplifying the restriction fragments via pcr, separating them through gel electrophoresis and analyzing them via southern blot with probes for the region of interest. this technique was used to compare the dna methylation modification patterns in the human adult and hemoglobin gene loci. different regions of the gene (gamma delta beta globin) were known to be expressed at different stages of development. consistent with a role of dna methylation in gene repression, regions that were associated with high levels of dna methylation were not actively expressed. this method was limited not suitable for studies on the global methylation pattern, or ‘methylome’. even within specific loci it was not fully representative of the true methylation pattern as only those restriction sites with corresponding methylation sensitive and insensitive restriction assays could provide useful information. further complications could arise when incomplete digestion of dna by restriction enzymes generated false negative results. genome wide approaches dna methylation profiling on a large scale was first made possible through the restriction landmark genome scanning (rlgs) technique. like the locus-specific dna methylation assay, the technique identified methylated dna via its digestion methylation sensitive enzymes. however it was the use of two-dimensional gel electrophoresis that allowed be characterized on a broader scale. however it was not until the advent of microarray and next generation sequencing technology when truly high resolution and genome-wide dna methylation became possible. as with rlgs, the endonuclease component is retained in the method but it is coupled to new technologies. one such approach is the differential methylation hybridization (dmh), in which one set of genomic dna is digested with methylation-sensitive restriction enzymes and a parallel set of dna is not digested. both sets of dna are subsequently amplified and each labelled with fluorescent dyes and used in two-colour array hybridization. the level of dna methylation at a given loci is determined by the relative intensity ratios of the two dyes. adaptation of next generation sequencing to dna methylation assay provides several advantages over array hybridization. sequence-based technology provides higher resolution to allele specific dna methylation, can be performed on larger genomes, and does not require creation of dna microarrays which require adjustments based on cpg density to properly function. bisulfite sequencing bisulfite sequencing relies on chemical conversion of unmethylated cytosines exclusively, such that they can be identified through standard dna sequencing techniques. sodium bisulfate and alkaline treatment does this by converting unmethylated cytosine residues into uracil while leaving methylated cytosine unaltered. subsequent amplification and sequencing of untreated dna and sodium bisulphite treated dna allows for methylated sites to be identified. bisulfite sequencing, like the traditional restriction based methods, was historically limited to methylation patterns of specific gene loci, until whole genome sequencing technologies became available. however, unlike traditional restriction based methods, bisulfite sequencing provided resolution on a nucleotide level. limitations of the bisulfite technique include the incomplete conversion of cytosine to uracil, which is a source of false positives. further, bisulfite treatment also causes dna degradation and requires an additional purification step to remove the sodium bisulfite. next-generation sequencing is well suited in complementing bisulfite sequencing in genome-wide methylation analysis. while this now allows for methylation pattern to be determined on the highest resolution possible, on the single nucleotide level, challenges still remain in the assembly step because of reduced sequence complexity in bisulphite treated dna. increases in read length seek to address this challenge, allowing for whole genome shotgun bisulphite sequencing (wgbs) to be performed. the wgbs approach using an illumina genome analyzer platform and has already been implemented in arabidopsis thaliana. reduced representation genomic methods based on bisulfite sequencing exist as well, and they are particularly suitable for species with large genome sizes. chromatin accessibility assays chromatin accessibility is the measure of how "accessible" or "open" a region of genome is to transcription or binding of transcription factors. the regions which are inaccessible (i.e. because they're bound by nucleosomes) are not actively transcribed by the cell while open and accessible regions are actively transcribed. changes in chromatin accessibility are important epigenetic regulatory processes that govern cell- or context-specific expression of genes. assays such as mnase-seq, dnase-seq, atac-seq or faire-seq are routinely used to understand the accessible chromatin landscape of cells. the main feature of all these methods is that they're able to selectively isolate either the dna sequences that are bounded to the histones, or those that are not. these sequences are then compared to a reference genome that allows to identify their relative position. mnase-seq and dnase-seq both follow the same principles, as they employ lytic enzymes that target nucleic acids to cut the dna strands unbounded by nucleosomes or other proteic factors, while the bounded pieces are sheltered, and can be retrieved and analysed. since active, unbound regions are destroyed, their detection can only be indirect, by sequencing with a next generation sequencing technique and comparison with a reference. mnase-seq utilises a micrococcal nuclease that produces a single strand cleavage on the opposite strand of the target sequence. dnase-seq employs dnase i, a non-specific double strand-cleaving endonuclease. this technique has been used to such an extent that nucleosome-free regions have been labelled as dhss, dnase i hypersensitive sites, and has been encode consortium's election method for genome wide chromatin accessibility analyses. the main issue of this technique is that the cleavage distribution can be biased, lowering the quality of the results. faire-seq (formaldehyde-assisted isolation of regulatory elements) requires as its first step crosslinking of the dna with nucleosomes, then dna shearing by sonication. the free and linked fragments are separated with a traditional phenol-chloroform extraction, since the proteic fraction is stuck in the interphase while the unlinked dna shifts to the aqueous phase and can be analysed with various methods. sonication produces random breaks, and therefore is not subject to any kind of bias, and is also the bigger length of the fragments (200-700 nt) makes this technique suitable for wider regions, while it's unable to resolve the single nucleosome. unlike the nuclease-based methods, faire-seq allows the direct identification of the transcriptionally active sites, and a less laborious sample preparation. atac-seq is based on the activity of tn5 transposase. the transposase is used to insert tags in the genome, with higher frequency on regions not covered by proteic factors. the tags are then used as adapters for prc or other analytical tools. direct detection polymerase sensitivity in single-molecule real-time sequencing made it possible for scientists to directly detect epigenetic marks such as methylation as the polymerase moves along the dna molecule being sequenced. several projects have demonstrated the ability to collect genome-wide epigenetic data in bacteria. nanopore sequencing is based on changes of electrolytic current signals according to base modifications (e.g. methylation). a polymerase mediates the entrance of ssdna in the pore: the ion-current variation is modulated by a section of the pore and the consequently generated difference is recorded revealing the position of cpg. discrimination between hydroxymethylation and methylation is possible thanks to solid-state nanopores even if the current while passing through the high-field region of the pore may be slightly influenced in it. as a reference amplified dna is used which will not present copied methylationed sites after the pcr process. the oxford nanopore technologies minion sequencer is a technology where, according to a hidden markov model, it is possible to distinguish unmethylated cytosine from the methylated one even without chemical treatment that acts to enhance the signal of that modification. the data are registered commonly in picoamperes during established time. other devices are the nanopolish and the signaalign: the former expresses the frequency of a methylation in a read while the latter gives a probability of it derived from the sum of all the reads. single-molecule real-time sequencing (smrt) is a single-molecule dna sequencing method. single-molecule real-time sequencing utilizes a zero-mode waveguide (zmw). a single dna polymerase enzyme is bound to the bottom of a zmw with a single molecule of dna as a template. each of the four dna bases is attached to one of four different fluorescent dyes. when a nucleotide is incorporated by the dna polymerase, the fluorescent tag is cleaved off and the detector detects the fluorescent signal of the nucleotide incorporation. as the sequencing occurs, the polymerase enzyme kinetics shift when it encounters a region of methylation or any other base modification. when the enzyme encounters chemically modified bases, it will slow down or speed up in a uniquely identifiable way. fluorescence pulses in smrt sequencing are characterized not only by their emission spectra but also by their duration and by the interval between successive pulses. these metrics, defined as pulse width and interpulse duration (ipd), add valuable information about dna polymerase kinetics. pulse width is a function of all kinetic steps after nucleotide binding and up to fluorophore release, and ipd is determined by the kinetics of nucleotide binding and polymerase translocation. in 2010 a team of scientists demonstrated the use of single-molecule real-time sequencing for direct detection of modified nucleotide in the dna template including n6-methyladenosine, 5-methylcytosine and 5-hydroxylcytosine. these various modifications affect polymerase kinetics differently, allowing discrimination between them. in 2017, another team proposed a combined bisulfite conversion with third-generation single-molecule real-time sequencing, it is called single-molecule real-time bisulfite sequencing (smrt-bs), which is an accurate targeted cpg methylation analysis method capable of a high degree of multiplying and long read lengths (1.5 kb) without the need for pcr amplicon sub-cloning. theoretical modeling approaches first mathematical models for different nucleosome states affecting gene expression were introduced in 1980s . later, this idea was almost forgotten, until the experimental evidence has indicated a possible role of covalent histone modifications as an epigenetic code. in the next several years, high-throughput data have indeed uncovered the abundance of epigenetic modifications and their relation to chromatin functioning which has motivated new theoretical models for the appearance, maintaining and changing these patterns,. these models are usually formulated in the frame of one-dimensional lattice approaches.