Long-time behaviour of a multidimensional age-dependent branching process with a singular jump kernel modelling telomere shortening
Jules Olayé, Milica Tomasevic
TL;DR
This work analyzes the long-time behavior of a multidimensional age-structured branching process with a singular, non-absolute continuous jump kernel modeling telomere shortening. By combining stochastic comparisons with Bellman-Harris renewal theory and a weak form of Harnack inequality, the authors establish exponential ergodicity of the first moment semigroup after an exponential normalisation, despite an unbounded birth rate and a non-compact kernel. A weighted-normalised semigroup and an auxiliary Markov process are constructed to obtain a precise stationary profile $(N,oldsymbol{ u},oldsymbol{ ho})$, along with a density representation that separates telomere-length and age components. The paper provides concrete, verifiable conditions under which the ergodic result holds and illustrates two models where the assumptions are satisfied, linking the mathematics closely to telomere biology and suggesting applicability to other age-structured branching systems. Overall, the results offer rigorous verification of stationary telomere-length distributions in complex, biologically realistic branching processes and supply a framework for exploring further Lyapunov-type conditions and large-number limits in such systems.
Abstract
In this article, we investigate the ergodic behaviour of a multidimensional age-dependent branching process with a singular jump kernel, motivated by studying the phenomenon of telomere shortening in cell populations. Our model tracks individuals evolving within a continuous-time framework indexed by a binary tree, characterised by age and a multidimensional trait. Branching events occur with rates dependent on age, where offspring inherit traits from their parent with random increase or decrease in some coordinates, while the most of them are left unchanged. Exponential ergodicity is obtained at the cost of an exponential normalisation, despite the fact that we have an unbounded age-dependent birth rate that may depend on the multidimensional trait, and a non-compact transition kernel. These two difficulties are respectively treated by stochastically comparing our model to Bellman-Harris processes, and by using a weak form of a Harnack inequality. We conclude this study by giving examples where the assumptions of our main result are verified.
