Fisher Information and the Dynamics of Multicellular Ageing

TL;DR

This paper uses Fisher information from information theory to model how multicellular organisms age through the balance of cell competition and cooperation.

Contribution

A novel theoretical model of multicellular ageing using Fisher information to explain dynamics of cell competition and cooperation.

Findings

Fisher information parametrizes the tension between cell reproduction and multicellular homeostasis through non-monotonic dynamics.

Optimal ageing dynamics depend on a balance of competition and cooperation among cells.

Fisher information suggests biological age inference evolves through complex dynamics over an organism's lifespan.

Abstract

Information theory has long been integrated into the study of biological ageing, for example, in examining the roles of genetic and epigenetic fidelity in cellular and organismal longevity. Here, we introduce a theoretical model that interprets ageing in multicellular systems through the lens of Fisher information. Previous theories have suggested that the ageing of multicellular organisms is an inevitable consequence of the inherent tension between individual cell reproduction and the homeostasis of the multicellular system. Utilising concepts from information theory and statistical mechanics, we show that Fisher information parametrises the dynamics of this tension through non-monotonic behaviour, which depends on an optimal balance of competition and cooperation between cells. Moreover, Fisher information suggests that the ability to infer true biological age from a sample evolves through complex dynamics over an organism’s lifespan.

Figures

Full text

1. Introduction

Multicellular life faces a fundamental conflict between individual cell competition and intercellular cooperation, whereby cells devote some degree of their individual resources to the common good [1,2]. Cooperative strategies can be summarised in five groups which benefit the whole but at a cost to individual cells: inhibition of cell proliferation, transport of resources, programmed cell death, division of labour, and creation and maintenance of the extracellular environment [2,3,4,5,6,7]. These processes are necessary and advantageous to multicellular life, requiring that cells sacrifice some measure of resources to promote the functioning and maintenance of the organism. However, cooperation between cells creates networks of interdependence characterised by shared resources and information, which in turn engenders potential for damage accumulation and ageing. Previous work has shown that interdependence within generalised networks of failure-prone components leads to characteristic ageing dynamics and, eventually, system collapse in both living and engineered systems [8,9,10,11,12,13]. An ubiquitous solution to mitigate ageing-induced degradation involves renewal and repair. In multicellular organisms, one such type of maintenance involves controlled cell proliferation. This mechanism, however, opens opportunities for deleterious somatic mutations, which can lead to a gradual deterioration of cellular functions and ultimately culminate in cellular senescence [14,15,16]. Senescent cell accumulation can be harmful to the vitality of the organism and has been linked with various age-related disorders such as cardiovascular and neurodegenerative diseases [16,17].

Competition between cells offers a mechanism to eliminate unfit senescent cells, potentially extending organismal vitality [2,16,17,18,19,20,21,22]. However, excessive intercellular competition can promote the proliferation of non-cooperative, “cheater” cells, such as cancer cells, which ultimately undermines organismal vitality [1,2,23,24,25,26,27]. These “cheater” cells enjoy a selective advantage over cooperative cells by eschewing investments in cooperative traits that benefit the multicellular organism but reduce their individual fitness [19,28]. Given these dynamics, the emergence of such cheaters and the consequent decline in cellular cooperation appear inescapable [14,29,30,31,32].

2. Ageing as Optimal Competition

A recent study [14] explores cancer and senescence as inevitable consequences of the dilemma of competition, and acts as our conceptual basis for modelling ageing. Essentially, competition is a double-bind where either cancer or senescence is favoured depending on the strength of competition. Organisms have therefore evolved an optimal level of competition which eliminates senescent cells but delays the onset of cancer, temporarily prolonging overall vitality. This interplay is a general model of ageing which relates competition, cancer, and ageing, where ageing is an intrinsic and inevitable characteristic of multicellular life. It is also compatible with existing explanations where ageing is a failure to select alleles which preserve the organism in late life; however, in this case, ageing is inevitable even if selection is perfect [14].

In our previous work [33], we introduced a minimal analytic framework using master equations to encapsulate these dynamics of multicellular ageing, where ageing is a problem of optimal competition. This is essentially a toy model which expresses key dynamics in the problem of optimal competition. Here, we show that this model can be reformulated in Lagrangian form. The resulting Lagrangian corresponds to Fisher information, thus providing an interpretation of this model of multicellular ageing in terms of information theory. Specifically, we introduce Fisher information to the study of ageing-as-optimal-competition, introducing a novel kind of information in ageing, to which biological relevance may be attributed.

3. Master Equation for the Dynamics of Multicellular Ageing

We consider a multicellular organism with cells characterised by two traits based on the notion that cells use resources to maintain their own fitness, but also sacrifice resources to cooperate: vigour v and cooperation c [14,33]. We use vigour v as a general measure of the resources available to the cell or cellular function. Healthy cells have high vigour, while senescent cells have low vigour. Cooperation c measures the fraction of resources that a cell devotes to cooperative activities necessary for the functioning of the organism (division of labour, resource allocation, creation and maintenance of the extracellular environment, proliferation inhibition and controlled cell death, etc. [1,2]). Healthy cells are highly cooperative, while cancer cells have low cooperation. We assume that vigour and cooperation take discrete values $[eqn]$ and $[eqn]$ , defining a two-dimensional coordinate system of cell types $[eqn]$ (Figure 1). On this coordinate system, “healthy” cells (“h”) are highly functioning and cooperative, thus have high vigour ( $[eqn]$ ) and high cooperation ( $[eqn]$ ). Cooperative cells with low vigour are “senescent” (“s”), while vigorous cells with low cooperation are “cancerous” (“c”). Cells with low vigour and cooperation are classified as both senescent and cancerous (“b”).

We describe the system dynamics on this two-dimensional coordinate system as a the result of (i) cell proliferation (via intercellular competition) with replication rate $[eqn]$ , which depends on state $[eqn]$ according to Equation (3), and (ii) somatic mutations that affect vigour or cooperation. We assume that mutations change either v or c in single steps with rates $[eqn]$ and $[eqn]$ for senescence- and cancer-causing mutations, respectively, as in our previous work, where we found that results from a continuous model were not qualitatively different [33]. Moreover, mutations in v and c are independent from each other, and we do not account for complex genetic interactions, e.g., epistatic interactions.

We use a master equation to describe the time evolution of the fraction $[eqn]$ of cells in state $[eqn]$ at time t [33]:

$$[eqn]$$

where

$$[eqn]$$

is the population averaged competition. The terms on the first line of Equation (1) describe the effect of selection via intercellular competition: the fraction of cells with $[eqn]$ larger than the average value $[eqn]$ will increase over time, while cells with $[eqn]$ will be purged away by selection. The terms on the second and third lines of Equation (1) describe the effect of somatic mutations that lower vigour and cooperation, respectively, with mutation rates $[eqn]$ and $[eqn]$ . We choose $[eqn]$ values to be much greater than $[eqn]$ , reflecting the fact that only about one percent of human genes are linked to cancer risk [14].

To close the system of equations, we need to specify a relationship for the replication rate $[eqn]$ as a function of $[eqn]$ which reflects the competitive strength of different phenotypes. We choose a linear fitness function

$$[eqn]$$

where k is the cost of cooperation and $[eqn]$ is the highest fitness. With this choice, cancer cells with $[eqn]$ and $[eqn]$ have the highest fitness, $[eqn]$ . Senescent cells with $[eqn]$ and $[eqn]$ have zero replication rate, $[eqn]$ , while healthy cells with $[eqn]$ and $[eqn]$ have intermediate fitness, $[eqn]$ . This choice effectively models the selective advantage of cells that either have access to high levels of resources (high vigour) or minimise resource investment in cooperative functions (low cooperation), reflecting the selection pressure due to competition that favours less cooperative “cheater” cells. The fitness landscape is therefore skewed towards cancerous phenotypes, highlighting that cheating provides a benefit in competition.

Time Evolution of Cell Fractions and the Inevitability of Multicellular Ageing

The population dynamic model Equation (1) describes multicellular ageing as an inevitable consequence of the trade-offs between cell competition and cooperation, where the selective pressures for individual cell proliferation conflict with those for maintaining organismal vitality, an idea discussed in [14]. Figure 1a shows an example of the time evolution of cell fractions for a four-state model, where $[eqn]$ . In this case, cells can be in one of four states: healthy (‘h’), senescent (‘s’), cancerous (‘c’), and both senescent and cancerous (‘b’). The dynamics of the four-state cell model unfold in two distinct phases due to the separation of timescales between mutation rates affecting senescence and cancer development. Initially, we observe a first phase where senescent cells accumulate, driven by senescence-causing mutations. Subsequently, the system transitions into a slower second phase, where senescent cells are incrementally eliminated due to the onset of cancer-causing mutations and the selection pressures from non-cooperative cancer cells. As a result, the fraction of cancerous cells begins to increase, displaying a sigmoidal growth curve that eventually plateaus. At the end, healthy and senescent cells always tend to zero, leaving a steady-state of only cells in states c and b.

4. Lagrangian Formulation of Master Equation

The master Equation (1) takes the form of a generalised replicator equation with mutation terms that change cell state. We now demonstrate that this equation can be derived from a Lagrangian function that admits an interpretation in terms of the Fisher information.

For this discussion, it is convenient to label two-dimensional cell states $[eqn]$ using a single index i in a $[eqn]$ -long state vector with components $[eqn]$ (Figure 1b). Using this notation, Equation (1) can be written as

$$[eqn]$$

where $[eqn]$ are the entries of a $[eqn]$ stochastic matrix $[eqn]$ with

$$[eqn]$$

The off-diagonal terms in $[eqn]$ describe mutation transitions from state j to state i, whereas the diagonal terms account for conservation of probability (Equation (5)). For example, for the four-state model, $[eqn]$ is a $[eqn]$ -matrix with entries

$$[eqn]$$

Populations are ordered h, s, c, b (Figure 1b), such that $[eqn]$ is the mutation rate from h to s, etc.

In order to derive the Lagrangian form of Equation (1), we introduce a composite proliferation rate $[eqn]$ combining the effect of intercellular competition $[eqn]$ and mutations $[eqn]$ as [34,35]

$$[eqn]$$

This transformation allows us to rewrite Equation (1) as an effective replicator equation with replication rate $[eqn]$ (see Appendix A):

$$[eqn]$$

This result is consistent with previous work that showed that any dynamical system on the $[eqn]$ dimensional probability simplex

$$[eqn]$$

can be formulated as replicator dynamics [36,37]. Using this transformation and after some algebra (see Appendix B), we can reformulate Equation (8) as the Euler–Lagrange equation for the following Lagrangian [34]:

$$[eqn]$$

subject to the constraint $[eqn]$ (conservation of probability). This Lagrangian takes the classical form $[eqn]$ , where

$$[eqn]$$

is a kinetic energy term and

$$[eqn]$$

is a potential energy term.

The kinetic energy term T admits an interesting geometric interpretation. It can be written as

$$[eqn]$$

where

$$[eqn]$$

is the length of a vector $[eqn]$ in the Fisher–Rao metric, also known as the Fisher information metric, a fundamental concept in information geometry used to measure distances on the probability simplex $[eqn]$ [38,39]. Moreover, information rate takes the distance, i.e., the relative entropy, between two sequential states of a time-dependent p.d.f., effectively measuring the rate of change of the p.d.f. in a path-dependent manner. The information rate $[eqn]$ measures how quickly new information is revealed as a p.d.f. evolves in time [40], and it can be written as

$$[eqn]$$ $$[eqn]$$

where J is the Jensen divergence. For more detail, refer to [40].

The dimension of $[eqn]$ is $[eqn]$ , so, by using units where the length is dimensionless, the information rate can be seen as the kinetic energy per unit mass. Furthermore, $[eqn]$ can be seen as the Fisher information if time is interpreted as an external control parameter [40]. In this context, information is revealed as cell populations evolve on a characteristic timescale $[eqn]$ [40].

Similarly, the potential energy term can be written as

$$[eqn]$$

where

$$[eqn]$$

denotes the average of the vector $[eqn]$ over the probability distribution $[eqn]$ . The potential energy V therefore measures the variance of the combined replication rate $[eqn]$ over the cell state distribution $[eqn]$ .

5. Link to Fisher Information

We now discuss the link between the Lagrangian Equation (10) and Fisher information. As well as information geometry, Fisher information is a key concept in statistical estimation theory, quantifying the amount of information that an observable random variable carries about an unknown parameter upon which the probability distribution of the random variable depends [38,39]. For a parametric model where the probability distribution $[eqn]$ of states $[eqn]$ varies smoothly with a control parameter $[eqn]$ , the Fisher information is defined as

$$[eqn]$$

where $[eqn]$ is the average of X over the probability distribution $[eqn]$ .

Interestingly, the kinetic energy term T (Equation (11)) corresponds to the Fisher information associated with the cell state distribution $[eqn]$ and the control parameter $[eqn]$ representing time $[eqn]$ :

$$[eqn]$$

Similarly, using the master equation $[eqn]$ , we can rewrite the potential energy term (Equation (12)) as

$$[eqn]$$

The Lagrangian Equation (10) is therefore recognised as twice the Fisher information (Figure 2b):

$$[eqn]$$

such that Fisher information is a representation of model dynamics. Fisher information is therefore a representation of optimal competition. In the problem of ageing, organisms have evolved an optimal level of competition which prolongs vitality. Depending on the cost of cooperation k, there is an optimal strength of competition $[eqn]$ , which maximises organismal vitality at that time, denoted as scenario O. This shown by the black line in Figure 3a at time $[eqn]$ . If cooperation is more costly, then the organism cannot tolerate higher levels of competition. The region towards higher competition is dominated by cancer, denoted as scenario C, and the region towards lower competition is dominated by senescence, denoted as scenario S. From Figure 3b, we see that Fisher information evolves in a non-monotonic manner, expressing dynamics dependent on which scenario we consider. Notably, Fisher information in the optimal case (scenario O) displays this same delay before the onset of cancer and before cell demographics reach a steady state. This transient feature is characteristic of our model and is also expressed by Fisher information. We now turn back to Fisher information in the context of statistical estimation.

6. Fisher Information and Statistical Inference

A key application of Fisher information in statistical estimation is the Cramér–Rao inequality. This states that if $[eqn]$ is an unbiased estimator of $[eqn]$ , then the variance of the estimation is greater than or equal to the inverse of Fisher information:

$$[eqn]$$

More information about the parameter $[eqn]$ means a lower bound on the variance of the estimator $[eqn]$ , and thus higher precision. In our model, Fisher information quantifies how much information a sample of cell fractions $[eqn]$ provides about current time t (Figure 4). Figure 2b shows the time evolution of inverse Fisher information $[eqn]$ for the four-state model. The evolution of $[eqn]$ displays features which represent key phases in the model: initial accumulation of senescent cells, a transition from senescence to cancer, and late-stage proliferation of cancer at late times. At $[eqn]$ , we have complete knowledge about the initial condition where all cells are healthy; thus, the inverse Fisher information vanishes. A detection of cancerous or senescent cells would indicate that cell populations have evolved beyond this initial state. Therefore, as senescent cells start to accumulate, the inverse Fisher information increases. A transient decrease in inverse Fisher information then follows during the transition from senescence- to cancer-domination, when the rate of population change is the greatest. According to Equation (23), this minimum in $[eqn]$ suggests an optimal inference of t. Finally, inverse Fisher information diverges as cancer cells proliferate and, eventually, we reach a steady state between c and b cells. This can be explained by cell demographics reaching a steady-state (thus, zero information rate), and also by the inability to infer true time during the steady state based on a measurement (Fisher information and statistical inference).

7. Conclusions

Life, in its essence, is information. It is therefore unsurprising that ageing has been extensively studied through the lens of information theory [41,42,43,44,45,46]. For example, the “Information Theory of Ageing” posits that this loss of information is primarily in the form of epigenetic information critical for maintaining cellular function and identity [41,47]. Here, we have built upon a model of multicellular ageing based on the principle of optimal intercellular cooperation and competition [14,33] and shown that its dynamics can be understood in terms of Fisher information. Fisher information introduces an information-geometric interpretation which has been previously applied to replicator dynamics, natural selection, and the abstract Price equation [34,39,48], which we extended here to multicellular ageing. Identifying Fisher information allows us to consider inference and sampling via the Cramer–Rao inequality in the context of multicellular ageing; we are not suggesting any causal relationship. Conventionally, we understand the role of information in ageing by the principle that the disorder of the genome and epigenome only increases with age [46,49]. This “drift” can be used to estimate a metric of biological age, which may differ from an organism’s chronological age based on congenital conditions, lifestyle, and environment [50,51,52]. Lost information in our model relates to the ensemble of cell states composing a tissue, where uncertainty is inherent to statistical sampling, rather than damage to biological code. Our work provides new insights into the role of information in ageing, which could inspire potential strategies to extend healthy lifespans in multicellular organisms.