  In our previous studies, we developed discrete-space Birth, Death and
Innovation Models (BDIM) of genome evolution. These models explain the origin
of the characteristic Pareto distribution of paralogous gene family sizes in
genomes, and model parameters that provide for the evolution of these
distributions within a realistic timeframe have been identified. Here we
develop the diffusion version of BDIM whose dynamics is described by the
Fokker-Plank equation and the stationary solution could be any specified Pareto
function. The diffusion models have time-dependent solutions of a special kind,
namely, the generalized self-similar solutions, which describe the transition
from one stationary distribution of the system to another; this provides for
the possibility of examining the temporal dynamics of genome evolution.
Analysis of the generalized self-similar solutions of the diffusion BDIM
reveals a biphasic curve of genome growth in which the initial, relatively
short, self-accelerating phase is followed by a prolonged phase of slow
deceleration. In biological terms, this regime of evolution can be tentatively
interpreted as a punctuated-equilibrium-like phenomenon such that whereby
evolutionary transitions are accompanied by rapid gene amplification and
innovation, followed by slow relaxation to a new stationary state.
