Return to search

On the Diffusion Approximation of Wright–Fisher Models with Several Alleles and Loci and its Geometry / Die Diffusionsnäherung von Wright-Fisher-Modellen mit mehreren Allelen und Loci und ihre Geometrie

The present thesis is located within the context of the diffusion approximation of Wright–Fisher models and the Kolmogorov equations describing their evolution. On the one hand, a full account of recombinational Wright–Fisher model is developed as well as their enhancement by other evolutionary mechanisms, including some information geometrical analysis. On the other hand, the thesis addresses several issues arising in the context of analytical solution schemes for such Kolmogorov equations, namely the inclusion of the entire boundary of the state space. For this, a hierarchical extension scheme is developed, both for the forward and the backward evolution, and the uniqueness of such extensions is proven.

First, a systematic approach to the diffusion approximation of recombinational two- or more loci Wright–Fisher models is presented. As a point of departure a specific Kolmogorov backward equation for the diffusion approximation of a recombinational two-loci Wright–Fisher model is chosen, to which – with the help of some information geometrical methods, i. e. by calculating the sectional curvatures of the corresponding statistical manifold (which is the domain equipped with the corresponding Fisher metric) – one succeeds to identify the underlying Wright–Fisher model. Accompanying this, for all methods and tools involved a suitable introduction is presented. Furthermore, the considerations span a separate analysis for the two most common underlying models (RUZ and RUG) as well as a comparison of the two models. Finally, transferring corresponding results for a simpler model described by Antonelli and Strobeck, solutions of the Kolmogorov equations are contrasted with Brownian motion in the same domain.

Furthermore, the perspective of the diffusion approximation of recombinational Wright–Fisher models is widened as the model underlying the Ohta–Kimura formula is subsequently extended by an integration of the concepts of natural fitness and mutation. Simultaneously, the corresponding extensions of the Ohta–Kimura formula are stated. Crucial for this is the development of a suitable fitness scheme, which is accomplished by a multiplicative aggregation of fitness values for pairs of gametes/zygotes. Furthermore, the model is generalised to have an arbitrary number of alleles and – in the following step – an arbitrary number of loci respectively. The latter involves an increased number of recombination modes, for which the concept of recombination masks is also implemented into the model. Another generalisation in terms of coarse-graining is performed via an application of schemata; this also affects the previously introduced concepts, specifically mask recombination, which are adapted accordingly. Eventually, a geometric analysis of linkage equilibrium states of the multi-loci Wright–Fisher models is carried out, relating to the concept of hierarchical probability distributions in information geometry, which concludes the considerations of recombinational Wright–Fisher models and their extensions.

Subsequently, the discussion of analytical solution schemes for the Kolmogorov equations corresponding to the diffusion approximation of Wright–Fisher models is ushered in, which represents the second part of the thesis. This is started with the simplest setting of a 1-dimensional Wright–Fisher model, for which the solution strategy for the corresponding Kolmogorov forward equation given by M. Kimura is recalled. From this, one may construct a unique extended solution which also accounts for the dynamics of the model on lower-dimensional entities of the state space, i. e. configurations of the model where one of the alleles no longer exists in the population, utilising the concept of (boundary) flux of a solution; a discussion of the moments of the distribution confirms the findings. A similar treatment is then carried out for the corresponding Kolmogorov backward equation, yielding analogous results of existence and uniqueness for an extended solution. For the latter in particular, a corresponding account of the configuration on the boundary turns out to be crucial, which is also reflected in the probabilistic interpretation of the backward solution. Additionally, the long-term behaviour of solutions is analysed, and a comparison between such solutions of the forward and the backward equation is made.

Next, it is basically aimed to transfer the results obtained in the previous chapter to the subsequent increasingly complicated setting of a Wright–Fisher model with 1 locus and an arbitrary number of alleles: With solution schemes for the interior of the state space (i. e. not encompassing the boundary) already existing in the literature, an extension scheme for a successive determination of the solution on lower-dimensional entities of the domain is developed. This scheme, again, makes use of the concept of the (boundary) flux of solutions, and one may therefore show that this extended solution fulfils additional properties regarding the completeness of the diffusion approximation with respect to the boundary. These properties may be formulated in terms of the moments of the distribution, and their connection to the underlying Wright–Fisher model is illustrated. Altogether, stipulating such a moments condition, existence and uniqueness of an extended solution on the entire domain are shown. Furthermore, the corresponding Kolmogorov backward equation is examined, for which similarly a (backward) extension scheme is presented, which allows extending a solution in a domain (perceived as a boundary instance of a larger domain) to all adjacent higher-dimensional entities of the larger domain along a certain path. This generalises the integration of boundary data observed in the previous chapter; in total, the existence of a solution of the Kolmogorov backward equation in the entire domain is shown for arbitrary boundary data.

Of particular interest to the discussion are stationary solutions of the Kolmogorov backward equation as they describe eventual hit probabilities for a certain target set of the model (in accordance with the probabilistic interpretation of solutions of the backward equation). The presented backward extension scheme allows the construction of solutions for all relevant cases, reconfirming some results by R. A. Littler for the stationary case, but now providing a previously missing systematic derivation. Eventually, the hitherto missing uniqueness assertion for this type of solutions is established by means of a specific iterated transformation which resolves the critical incompatibilities of solutions by a successive blow-up while the domain is converted from a simplex into a cube. Then – under certain additional assumptions on the regularity of the transformed solution – the uniqueness directly follows from general principles. Lastly, several other aspects of the blow-up scheme are discussed; in particular, it is illustrated in what way the required extra regularity relates to reasonable additional properties of the underlying Wright–Fisher model.

Identiferoai:union.ndltd.org:DRESDEN/oai:qucosa.de:bsz:15-qucosa-150016
Date25 September 2014
CreatorsHofrichter, Julian
ContributorsUniversität Leipzig, Fakultät für Mathematik und Informatik, Prof. Dr. Jürgen Jost, Prof. Dr. Jürgen Jost, Prof. Dr. Angela Stevens
PublisherUniversitätsbibliothek Leipzig
Source SetsHochschulschriftenserver (HSSS) der SLUB Dresden
LanguageEnglish
Detected LanguageEnglish
Typedoc-type:doctoralThesis
Formatapplication/pdf

Page generated in 0.0026 seconds