An inverse problem in cell dynamics: Recovering an initial distribution of telomere lengths from measurements of senescence times

Telomeres are repetitive sequences situated at both ends of the chromosomes of eukaryotic cells. At each cell division, they are eroded until they reach a critical length that triggers a state in which the cell stops to divide: the senescent state. In this work, we are interested in the link between the initial distribution of telomere lengths and the distribution of senescence times. We propose a method to retrieve the initial distribution of telomere lengths, using only measurements of senescence times. Our approach relies on approximating our models with transport equations, which provide natural estimators for the initial telomere lengths distribution. We investigate this method from a theoretical point of view by providing bounds on the errors of our estimators, pointwise and in all Lebesgue spaces. We also illustrate it with estimations on simulations, and discuss its limitations related to the curse of dimensionality.