Effects of non-CpG site methylation on DNA thermal stability: a fluorescence study

High-Resolution Melting (HRM) analysis of DNA methylation using semiconductor chip-based digital PCR

BACKGROUND

Digital PCR (dPCR) technology allows absolute quantification and detection of disease-associated rare variants, and thus the use of dPCR technology has been increasing in clinical research and diagnostics. The high-resolution melting curve analysis (HRM) of qPCR is widely used to distinguish true positives from false positives and detect rare variants. In particular, qPCR-HRM is commonly used for methylation assessment in research and diagnostics due to its simplicity and high reproducibility. Most dPCR instruments have limited fluorescence channels available and separate heating and imaging systems. Therefore, it is difficult to perform HRM analysis using dPCR instruments.

OBJECTIVE

A new digital real-time PCR instrument (LOAA) has been recently developed to integrate partitioning, thermocycling, and imaging in a single dPCR instrument. In addition, a new technique to perform HRM analysis is utilized in LOAA. The aim of the present study is to evaluate the efficiency and accuracy of LOAA dPCR on HRM analysis for the detection of methylation.

METHODS

In this study, comprehensive comparison with Bio-Rad qRT-PCR and droplet-based dPCR equipment was performed to verify the HRM analysis-based methylation detection efficiency of the LOAA digital PCR equipment. Here, sodium bisulfite modification method was applied to detect methylated DNA sequences by each PCR method.

RESULTS

Melting curve analysis detected four different Tm values using LOAA and qPCR, and found that LOAA, unlike qPCR, successfully distinguished between different Tm values when the Tm values were very similar. In addition, melting temperatures increased by each methylation were about 0.5℃ for qPCR and about 0.2 ~ 0.6℃ for LOAA. The melting temperature analyses of methylated and unmethylated DNA samples were conducted using LOAA dPCR with TaqMan probes and EvaGreen, and the result found that Tm values of methylated DNA samples are higher than those of unmethylated DNA samples.

CONCLUSION

The present study shows that LOAA dPCR could detect different melting temperatures according to methylation status of target sequences, indicating that LOAA dPCR would be useful for diagnostic applications that require the accurate quantification and assessment of DNA methylation.

Skeletal muscle feature of different populations in large yellow croaker (Larimichthys crocea): from an epigenetic point of view

Fish skeletal muscle is composed of well-defined fiber types. In order to identify potential candidate genes affecting muscle growth and development under epigenetic regulation. Bisulfite sequencing was utilized to analyze and compare the muscle DNA methylation profiles of Larimichthys crocea inhabiting different environments. The results revealed that DNA methylation in L. crocea was predominantly CG methylation, with 2,396 differentially methylated regions (DMRs) identified through comparisons among different populations. The largest difference in methylation was observed between the ZhouShan and JinMen wild populations, suggesting that L. crocea may have undergone selection and domestication. Additionally, GO and KEGG enrichment analysis of differentially methylated genes (DMGs) revealed 626 enriched GO functional categories, including various muscle-related genes such as myh10, myf5, myf6, ndufv1, klhl31, map3k4, syn2b, sostdc1a, bag4, and hsp90ab. However, significant enrichment in KEGG pathways was observed only in the JinMen and XiangShan populations of L. crocea. Therefore, this study provides a theoretical foundation for a better understanding of the epigenetic regulation of skeletal muscle growth and development in L. crocea under different environmental conditions.

Kinetic Referencing Allows Identification of Epigenetic Cytosine Modifications by Single-Molecule Hybridization Kinetics and Superresolution DNA-PAINT Microscopy

We develop a DNA origami-based internal kinetic referencing system with a colocalized reference and target molecule to provide increased sensitivity and robustness for transient binding kinetics. To showcase this, we investigate the subtle changes in binding strength of DNA oligonucleotide hybrids induced by cytosine modifications. These cytosine modifications, especially 5-methylcytosine but also its oxidized derivatives, have been increasingly studied in the context of epigenetics. Recently revealed correlations of epigenetic modifications and disease also render them interesting biomarkers for early diagnosis. Internal kinetic referencing allows us to probe and compare the influence of the different epigenetic cytosine modifications on the strengths of 7-nucleotide long DNA hybrids with one or two modified nucleotides by single-molecule imaging of their transient binding, revealing subtle differences in binding times. Interestingly, the influence of epigenetic modifications depends on their position in the DNA strand, and in the case of two modifications, effects are additive. The sensitivity of the assay indicates its potential for the direct detection of epigenetic disease markers.

Z-Form Adoption of Nucleic Acid is a Multi-Step Process Which Proceeds through a Melted Intermediate

The left-handed Z-conformation of nucleic acids can be adopted by both DNA and RNA when bound by Zα domains found within a variety of innate immune response proteins. Zα domains stabilize this higher-energy conformation by making specific interactions with the unique geometry of Z-DNA/Z-RNA. However, the mechanism by which a right-handed helix contorts to become left-handed in the presence of proteins, including the intermediate steps involved, is poorly understood. Through a combination of nuclear magnetic resonance (NMR) and other biophysical measurements, we have determined that in the absence of Zα, under low salt conditions at room temperature, d(CpG) and r(CpG) constructs show no observable evidence of transient Z-conformations greater than 0.5% on either the intermediate or slow NMR time scales. At higher temperatures, we observed a transient unfolded intermediate. The ease of melting a nucleic acid duplex correlates with Z-form adoption rates in the presence of Zα. The largest contributing factor to the activation energies of Z-form adoption as calculated by Arrhenius plots is the ease of flipping the sugar pucker, as required for Z-DNA and Z-RNA. Together, these data validate the previously proposed "zipper model" for Z-form adoption in the presence of Zα. Overall, Z-conformations are more likely to be adopted by double-stranded DNA and RNA regions flanked by less stable regions and by RNAs experiencing torsional/mechanical stress.

A qPCR technology for direct quantification of methylation in untreated DNA

DNA methylation is important for gene expression and alterations in DNA methylation are involved in the development and progression of cancer and other major diseases. Analysis of DNA methylation patterns has until now been dependent on either a chemical or an enzymatic pre-treatment, which are both time consuming procedures and potentially biased due to incomplete treatment. We present a qPCR technology, EpiDirect®, that allows for direct PCR quantification of DNA methylations using untreated DNA. EpiDirect® is based on the ability of Intercalating Nucleic Acids (INA®) to differentiate between methylated and unmethylated cytosines in a special primer design. With this technology, we develop an assay to analyze the methylation status of a region of the MGMT promoter used in treatment selection and prognosis of glioblastoma patients. We compare the assay to two bisulfite-relying, methyl-specific PCR assays in a study involving 42 brain tumor FFPE samples, revealing high sensitivity, specificity, and the clinical utility of the method.

On the thermodynamics of DNA methylation process

DNA methylation is an epigenetic mechanism that plays important roles in various biological processes including transcriptional and post-transcriptional regulation, genomic imprinting, aging, and stress response to environmental changes and disease. Consistent with thermodynamic principles acting within living systems and the application of maximum entropy principle, we propose a theoretical framework to understand and decode the DNA methylation process. A central tenet of this argument is that the probability density function of DNA methylation information-divergence summarizes the statistical biophysics underlying spontaneous methylation background and implicitly bears on the channel capacity of molecular machines conforming to Shannon's capacity theorem. On this theoretical basis, contributions from the molecular machine (enzyme) logical operations to Gibb entropy (S) and Helmholtz free energy (F) are intrinsic. Application to the estimations of S on datasets from Arabidopsis thaliana suggests that, as a thermodynamic state variable, individual methylome entropy is completely determined by the current state of the system, which in biological terms translates to a correspondence between estimated entropy values and observable phenotypic state. In patients with different types of cancer, results suggest that a significant information loss occurs in the transition from differentiated (healthy) tissues to cancer cells. This type of analysis may have important implications for early-stage diagnostics. The analysis of entropy fluctuations on experimental datasets revealed existence of restrictions on the magnitude of genome-wide methylation changes originating by organismal response to environmental changes. Only dysfunctional stages observed in the Arabidopsis mutant met1 and in cancer cells do not conform to these rules.

5-Methyl-cytosine stabilizes DNA but hinders DNA hybridization revealed by magnetic tweezers and simulations

5-Methyl-cytosine (5mC) is one of the most important DNA modifications and plays versatile biological roles. It is well known that 5mC stabilizes DNA duplexes. However, it remains unclear how 5mC affects the kinetics of DNA melting and hybridization. Here, we studied the kinetics of unzipping and rezipping using a 502-bp DNA hairpin by single-molecule magnetic tweezers. Under constant loading rates, 5mC increases the unzipping force but counterintuitively decreases the rezipping force at various salt and temperature conditions. Under constant forces, the non-methylated DNA hops between metastable states during unzipping and rezipping, which implies low energy barriers. Surprisingly, the 5mC DNA can't rezip after fully unzipping unless much lower forces are applied, where it rezips stochastically in a one-step manner, which implies 5mC kinetically hinders DNA hybridization and high energy barriers in DNA hybridization. All-atom molecular dynamics simulations reveal that the 5mC kinetically hinders DNA hybridization due to steric effects rather than electrostatic effects caused by the additional methyl groups of cytosines. Considering the possible high speed of DNA unzipping and zipping during replication and transcription, our findings provide new insights into the biological roles of 5mC.

DNA methylation: Precise modulation of chromatin structure and dynamics

DNA methylation plays a vital role in epigenetic regulation in both plants and animals, and typically occurs at the 5-carbon position of the cytosine pyrimidine ring within the CpG dinucleotide steps. Cytosine methylation can alter DNA's geometry, mechanical and physico-chemical properties - thus influencing the molecular signaling events vital for transcription, replication and chromatin remodeling. Despite the profound effect cytosine methylation can have on DNA, the underlying atomistic mechanisms remain enigmatic. Many studies so far have produced controversial findings on how cytosine methylation dictates DNA flexibility and accessibility, nucleosome stability and dynamics. Here, we review the most recent experimental and computational studies that provide precise characterization of structure and function of cytosine methylation and its versatile roles in modulating DNA mechanics, nucleosome and chromatin structure, stability and dynamics. Moreover, the review briefly discusses the relationship between DNA methylation and nucleosome positioning, and the crosstalk between DNA methylation and histone tail modifications.

Cytosine methylation regulates DNA bendability depending on the curvature

Cytosine methylation plays an essential role in many biological processes, such as nucleosome inactivation and regulation of gene expression. The modulation of DNA mechanics may be one of the regulatory mechanisms influenced by cytosine methylation. However, it remains unclear how methylation influences DNA mechanics. Here, we show that methylation has contrasting effects on the bending property of dsDNA depending on DNA curvature. We directly applied bending force on 30 base pairs of dsDNA using a D-shaped DNA nanostructure and measured the degree of bending using single-molecule fluorescence resonance energy transfer without surface immobilization. When dsDNA is weakly bent, methylation increases the stiffness of dsDNA. The stiffness of dsDNA increased by approximately 8% with a single methylation site for 30 bp dsDNA. When dsDNA is highly bent by a strong force, it forms a kink, i.e., a sharp bending of dsDNA. Under strong bending, methylation destabilizes the non-kink form compared with the kink form, which makes dsDNA near the kink region apparently more bendable. However, if the kink region is methylated, the kink form is destabilized, and dsDNA becomes stiffer. As a result, methylation increases the stiffness of weakly bent dsDNA and concurrently can promote kink formation, which may stabilize the nucleosome structure. Our results provide new insight into the effect of methylation, showing that cytosine methylation has opposite effects on DNA mechanics depending on its curvature and methylation location.

Effect of DNA modifications on the transition between canonical and non-canonical DNA structures in CpG islands during senescence

Patterns and levels of DNA modifications play important roles in senescence. Two major epigenetic modifications of DNA, 5-methylcytosine (5mC) and 5-hydroxymethylcytosine (5hmC), target CpG sites. Importantly, CpG concentrated regions, known as CpG islands, contain GC-rich sequences, which have the potential to fold into non-canonical DNA structures such as i-motifs and G-quadruplexes. In this study, we investigated the effect of 5mC and 5hmC modifications on the transition between a duplex, and i-motif and G-quadruplexes. To examine the transition, we firstly investigated the stability and structure of the i-motif and G-quadruplexes, considering the molecular environment in senescent cells. Analyses of their stability showed that the modifications did not drastically affect the stability. However, noteworthily, the modification can weaken the (de)stabilisation effect on G-quadruplexes caused by cosolute(s) and cations. Circular dichroism analyses indicated that the surrounding environments, including the molecular crowding and the type of cations such as K⁺ and Na⁺, regulate the topology of G-quadruplexes, while neither 5mC nor 5hmC had a drastic effect. On the other hand, the modifications changed the transition between duplexes and quadruplexes. Unmodified DNA preferred to fold into quadruplexes, whereas DNA with 5mC and 5hmC preferred to fold into duplexes in the absence of PEG200; on the other hand, DNA with or without modifications tended to fold into i-motifs under crowded conditions. Furthermore, an investigation of quadruplexes forming sequences in CpG islands, which are hyper- or hypomethylated during senescence, followed by gene ontology enrichment analysis for each gene group classified by the presence of quadruplexes, showed a difference in function between genes with and without quadruplexes in the CpG region. These results indicate that it is important to consider the effects of patterns and levels of DNA modifications on the transition between canonical and non-canonical DNA structures to understand gene regulation by epigenetic modification during senescence.

The modification of DNA can regulate the transition between a duplex and quadruplexes during senescence responding to surrounding environments.

Double-stranded flanking ends affect the folding kinetics and conformational equilibrium of G-quadruplexes forming sequences within the promoter of KIT oncogene

G-quadruplexes embedded within promoters play a crucial role in regulating the gene expression. KIT is a widely studied oncogene, whose promoter contains three G-quadruplex forming sequences, c-kit1, c-kit2 and c-kit*. For these sequences available studies cover ensemble and single-molecule analyses, although for kit* the latter were limited to a study on a promoter domain comprising all of them. Recently, c-kit2 has been reported to fold according to a multi-step process involving folding intermediates. Here, by exploiting fluorescence resonance energy transfer, both in ensemble and at the single molecule level, we investigated the folding of expressly designed constructs in which, alike in the physiological context, either c-kit2 or c-kit* are flanked by double stranded DNA segments. To assess whether the presence of flanking ends at the borders of the G-quadruplex affects the folding, we studied under the same protocols oligonucleotides corresponding to the minimal G-quadruplex forming sequences. Data suggest that addition of flanking ends results in biasing both the final equilibrium state and the folding kinetics. A previously unconsidered aspect is thereby unravelled, which ought to be taken into account to achieve a deeper insight of the complex relationships underlying the fine tuning of the gene-regulatory properties of these fascinating DNA structures.

Genome-wide profiling reveals functional interplay of DNA sequence composition, transcriptional activity, and nucleosome positioning in driving DNA supercoiling and helix destabilization in C. elegans

DNA topology and alternative DNA structures are implicated in regulating diverse biological processes. Although biomechanical properties of these structures have been studied extensively in vitro, characterization in vivo, particularly in multicellular organisms, is limited. We devised new methods to map DNA supercoiling and single-stranded DNA in Caenorhabditis elegans embryos and diapause larvae. To map supercoiling, we quantified the incorporation of biotinylated psoralen into DNA using high-throughput sequencing. To map single-stranded DNA, we combined permanganate treatment with genome-wide sequencing of induced double-stranded breaks. We found high levels of negative supercoiling at transcription start sites (TSSs) in embryos. GC-rich regions flanked by a sharp GC-to-AT transition delineate boundaries of supercoil propagation. In contrast to TSSs in embryos, TSSs in diapause larvae showed dramatic reductions in negative supercoiling without concomitant attenuation of transcription, suggesting developmental-stage-specific regulation. To assess whether alternative DNA structures control chromosome architecture and gene expression, we examined DNA supercoiling in the context of X-Chromosome dosage compensation. We showed that the condensin dosage compensation complex creates negative supercoils locally at its highest-occupancy binding sites but found no evidence for large-scale supercoiling domains along X Chromosomes. In contrast to transcription-coupled negative supercoiling, single-strandedness, which is most pronounced at transcript end sites, is dependent on high AT content and symmetrically positioned nucleosomes. We propose that sharp transitions in sequence composition at functional genomic elements constitute a common regulatory code and that DNA structure and propagation of torsional stress at regulatory elements are critical parameters in shaping important developmental events.

Single-molecule micromanipulation studies of methylated DNA

Cytosine methylated at the five-carbon position is the most widely studied reversible DNA modification. Prior findings indicate that methylation can alter mechanical properties. However, those findings were qualitative and sometimes contradictory, leaving many aspects unclear. By applying single-molecule magnetic force spectroscopy techniques allowing for direct manipulation and dynamic observation of DNA mechanics and mechanically driven strand separation, we investigated how CpG and non-CpG cytosine methylation affects DNA micromechanical properties. We quantitatively characterized DNA stiffness using persistence length measurements from force-extension curves in the nanoscale length regime and demonstrated that cytosine methylation results in longer contour length and increased DNA flexibility (i.e., decreased persistence length). In addition, we observed the preferential formation of plectonemes over unwound single-stranded "bubbles" of DNA under physiologically relevant stretching forces and supercoiling densities. The flexibility and high structural stability of methylated DNA is likely to have significant consequences on the recruitment of proteins recognizing cytosine methylation and DNA packaging.

Combined Effects of Methylated Cytosine and Molecular Crowding on the Thermodynamic Stability of DNA Duplexes

Methylated cytosine within CpG dinucleotides is a key factor for epigenetic gene regulation. It has been revealed that methylated cytosine decreases DNA backbone flexibility and increases the thermal stability of DNA. Although the molecular environment is an important factor for the structure, thermodynamics, and function of biomolecules, there are few reports on the effects of methylated cytosine under a cell-mimicking molecular environment. Here, we systematically investigated the effects of methylated cytosine on the thermodynamics of DNA duplexes under molecular crowding conditions, which is a critical difference between the molecular environment in cells and test tubes. Thermodynamic parameters quantitatively demonstrated that the methylation effect and molecular crowding effect on DNA duplexes are independent and additive, in which the degree of the stabilization is the sum of the methylation effect and molecular crowding effect. Furthermore, the effects of methylation and molecular crowding correlate with the hydration states of DNA duplexes. The stabilization effect of methylation was due to the favorable enthalpic contribution, suggesting that direct interactions of the methyl group with adjacent bases and adjacent methyl groups play a role in determining the flexibility and thermodynamics of DNA duplexes. These results are useful to predict the properties of DNA duplexes with methylation in cell-mimicking conditions.

Demonstration of fibrinogen-FcRn binding at acidic pH by means of Fluorescence Correlation Spectroscopy

The neonatal Fc receptor (FcRn) interacts with IgG and albumin at acidic pH within endosomes, thus protecting these plasma proteins from degradation. Recently, we proposed fibrinogen as a new binding partner of FcRn. This work was aimed at providing a direct demonstration of FcRn-fibrinogen binding at acidic pH by Fluorescence Correlation Spectroscopy. The increase in diffusion time between free and fibrinogen-bound FITC-labelled FcRn was assumed as the binding indicator. We observed that, at acidic pH (pH = 5.3), FcRn diffusion time shifted from ≈730 μs (FITC-labelled FcRn alone) to >1200 μs (FITC-labelled FcRn added with fibrinogen). A similar trend was exhibited by albumin, a known FcRn interactor, while no significant variations in diffusion time were observed upon incubation with catalase as negative control. Our results demonstrate a binding interaction between fibrinogen, one of the most abundant plasma proteins, and FcRn, a receptor involved in the regulation of the levels of IgG and albumin. This interaction is likely responsible for fibrinogen protection from intracellular degradation and recycling in plasma. Fibrinogen is crucial not only in haemostasis but also in acute inflammatory response and in some pathological conditions. The interaction with FcRn can influence not only the levels of fibrinogen in plasma and other tissues, but also the levels of other FcRn binding partners, among which are some plasma proteins of clinical relevance.

Solubilization of the chlorin TPCS2a in the presence of Pluronic® F127/Tween 80 mixtures

The efficacy of surfactant mixtures of Pluronic^® F127 and Tween 80 at overall concentration in the micromolar range and molar ratio 1:1, 1:10, and 10:1 in inhibiting aggregation of the photosensitizer meso-tetraphenyl chlorin disulphonate (TPCS_2a) was investigated in aqueous media at pH 2.9 by means of steady-state absorption and fluorescence emission spectroscopy as well as time-resolved fluorescence analysis. Corresponding experiments were performed at pH 7.4 in the absence of surfactants to determine the spectroscopic properties of a monomeric sample. Aggregation resulted in a red shift of the Soret absorption band and in substantial fluorescence quenching. The fluorescence lifetime of TPCS_2a was a particularly sensitive indicator of the aggregation state, as the monomer at pH 7.4 decayed with a ∼ 10 ns time constant, while aggregation resulted in subnanosecond decay. The critical micelle concentration (CMC) of the surfactant mixtures was determined spectrophotometrically in the presence of TPCS_2a. The ability of the surfactant mixtures to prevent aggregation at acidic pH was evaluated at overall surfactant concentration below and above CMC. Solubilization of TPCS_2a in Pluronic^® F127/Tween 80 mixtures prevented aggregation of the photosensitizer at overall surfactant concentrations much lower than those needed for both pure Pluronic^® F127 and pure Tween 80.

Optical Trapping Nanometry of Hypermethylated CPG-Island DNA

Cytosine methylation is a key mechanism of epigenetic regulation. CpG-dense loci, called "CpG islands", play a particularly important role in modulating gene expression. Methylation has long been suspected to alter the physical properties of DNA, but the full spectrum of the evoked changes is unknown. Here we measured the methylation-induced nanomechanical changes in a DNA molecule with the sequence of a CpG island. For the molecule under tension, contour length, bending rigidity and intrinsic stiffness decreased in hypermethylated dsDNA, pointing at structural compaction which may facilitate DNA packaging in vivo. Intriguingly, increased forces were required to convert hypermethylated dsDNA into an extended S-form configuration. The reduction of force hysteresis during mechanical relaxation indicated that methylation generates a barrier against strand unpeeling and melting-bubble formation. The high structural stability is likely to have significant consequences on the recognition, replication, transcription, and reparation of hypermethylated genetic regions.

Genome-Wide Discriminatory Information Patterns of Cytosine DNA Methylation

Cytosine DNA methylation (CDM) is a highly abundant, heritable but reversible chemical modification to the genome. Herein, a machine learning approach was applied to analyze the accumulation of epigenetic marks in methylomes of 152 ecotypes and 85 silencing mutants of Arabidopsis thaliana. In an information-thermodynamics framework, two measurements were used: (1) the amount of information gained/lost with the CDM changes I R and (2) the uncertainty of not observing a SNP L C R . We hypothesize that epigenetic marks are chromosomal footprints accounting for different ontogenetic and phylogenetic histories of individual populations. A machine learning approach is proposed to verify this hypothesis. Results support the hypothesis by the existence of discriminatory information (DI) patterns of CDM able to discriminate between individuals and between individual subpopulations. The statistical analyses revealed a strong association between the topologies of the structured population of Arabidopsis ecotypes based on I R and on LCR, respectively. A statistical-physical relationship between I R and L C R was also found. Results to date imply that the genome-wide distribution of CDM changes is not only part of the biological signal created by the methylation regulatory machinery, but ensures the stability of the DNA molecule, preserving the integrity of the genetic message under continuous stress from thermal fluctuations in the cell environment.

Information Thermodynamics of Cytosine DNA Methylation

Cytosine DNA methylation (CDM) is a stable epigenetic modification to the genome and a widespread regulatory process in living organisms that involves multicomponent molecular machines. Genome-wide cytosine methylation patterning participates in the epigenetic reprogramming of a cell, suggesting that the biological information contained within methylation positions may be amenable to decoding. Adaptation to a new cellular or organismal environment also implies the potential for genome-wide redistribution of CDM changes that will ensure the stability of DNA molecules. This raises the question of whether or not we would be able to sort out the regulatory methylation signals from the CDM background (“noise”) induced by thermal fluctuations. Here, we propose a novel statistical and information thermodynamic description of the CDM changes to address the last question. The physical basis of our statistical mechanical model was evaluated in two respects: 1) the adherence to Landauer’s principle, according to which molecular machines must dissipate a minimum energy ε = k_BT ln2 at each logic operation, where k_B is the Boltzmann constant, and T is the absolute temperature and 2) whether or not the binary stretch of methylation marks on the DNA molecule comprise a language of sorts, properly constrained by thermodynamic principles. The study was performed for genome-wide methylation data from 152 ecotypes and 40 trans-generational variations of Arabidopsis thaliana and 93 human tissues. The DNA persistence length, a basic mechanical property altered by CDM, was estimated with values from 39 to 66.9 nm. Classical methylome analysis can be retrieved by applying information thermodynamic modelling, which is able to discriminate signal from noise. Our finding suggests that the CDM signal comprises a language scheme properly constrained by molecular thermodynamic principles, which is part of an epigenomic communication system that obeys the same thermodynamic rules as do current human communication systems.