You are here : Home > News > A mathematical model applicable to nuclear processes

Scientific result | DNA | Cellular mechanisms

A mathematical model applicable to nuclear processes at the genome scale, and at very high spatial resolution


The protocols applied for in vivo mapping of nuclear processes, on a genome-wide scale and at very high spatial resolution do not allow for high temporal resolution. Thus, data sets usually contain only a few samples. A team from I2BC has shown that such an understanding can be wrong, using a mathematical model. Their study published in PLoS Computational Biology analyzes sequencing signals as the superposition of signals from independent cells. Their model can be applied to any process that can be modeled by a transition between two states.

Published on 13 December 2022

A potentially erroneous understanding

New technologies for rapid sequencing of large samples of DNA and RNA (NGS - Next Generation Sequencing), have made it possible to map nuclear processes in vivo, on a genome-wide scale and at very high spatial resolution. This is particularly the case for gene transcription and lesion repair.

While the protocols applied are efficient, they do not allow for high temporal resolution. Thus, the data sets generally contain only a few samples over several hours, and often the transition from one point to another is implicitly considered as direct. The "Transcriptional regulation of the genome" team of the I2BC therefore wanted to demonstrate that such an understanding could be wrong

A mathematical model for complete analysis

The team decided to submit previous sequencing data of cyclobutanic pyrimidine dimers (CPD) 1 obtained in the yeast Saccharomyces cerevisiae to a mathematical model. CPDs are one of the DNA damages induced by UVs, so their mapping provides information on the repair process. The model applied is a probabilistic model, which comes from physics 2, established in the 1940s to describe changes of state in solids.

Compared to other studies that consider sequencing signals as an average behavior, their study, published in PLoS Computational Biology, analyzes them as the superposition of signals from independent cells. Their repair model thus constructed leads to consider that the change observed in the sequencing signals of CPDs is not constant over time. Moreover, it allowed them to establish correlations with other nuclear processes/parameters such as transcription rate and nucleosomal density.

With this model, repair positions on the genome appear as "patterns" that develop within a population of cells. It thus allows a complete analysis of nuclear processes at the population level. As their model is validated by comparing their repair dynamics with independent sequencing data, it can be applied to any process that can be modeled by a transition between two states, such as repaired and damaged.

1 : by mapping cyclobutanic pyrimidine dimers (CPD)

2 : KJMA model for Kolmogorov, Johnson Mehl and Avrami

Contact Joliot : 

Julie Soutourina (julie.soutourina@cea.fr)


Top page