# Bernstein problem in mathematical genetics

Bernshtein problem

Let $\Delta ^ { n - 1 }$ be the simplex in ${\bf R} ^ { n }$ spanned by the canonical basis $\{ e _ { i } \} _ { 1 } ^ { n }$. Any set of numbers $p _ { i k , j } \geq 0$ ($1 \leq i , k , j \leq n$) such that $\sum _ { j } p _ { i k,j } = 1$ and $p _ { i k ,\, j} = p _ { k i ,\, j}$ defines a stochastic quadratic mapping $V:\Delta ^ { n - 1 } \rightarrow \Delta ^ { n - 1 }$ by the formulas

\begin{equation*} x _ { j } ^ { \prime } = \sum _ { i , k } p _ { i k,j } x_i x _ { k } , \quad x _ { i } \geq 0 , \sum _ { i } x _ { i } = 1. \end{equation*}

This mapping is called Bernstein (or stationary) if $V ^ { 2 } = V$. The Bernstein problem is to explicitly describe all such mappings. This problem was posed by S.N. Bernshtein [a1] in order to create a mathematical foundation of population genetics. For $n = 3$, this problem has been solved in [a2]. (For $n \leq 2$ the problem is trivial.)

The classical Mendel mechanism of heredity defines a mapping

\begin{equation*} x _ { 1 } ^ { \prime } = p ^ { 2 } , x _ { 2 } ^ { \prime } = q ^ { 2 } , x _ { 3 } ^ { \prime } = 2 p q, \end{equation*}

\begin{equation*} p = x _ { 1 } + \frac { 1 } { 2 } x _ { 3 } , \quad q = x _ { 2 } + \frac { 1 } { 2 } x _ { 3 } \end{equation*}

(the Hardy–Weinberg laws, cf. [a5]). Biologically, $p$ and $q$ are the probabilities of an alternating pair of genes, say $A$ and $a$ respectively, in a population where the individuals may be of genotypes $A A$, $aa$ and $A a$. Then $x_{1}$, $x _ { 2 }$, $x _ { 3 }$ are the probabilities of these genotypes in a generation. If the next generation is formed by random mating, then the probabilities turn into $x _ { 1 } ^ { \prime }$, $x _ { 2 } ^ { \prime }$, $x _ { 3 } ^ { \prime }$. As a result, $p ^ { \prime } = p$, $q ^ { \prime } = q$ and then $x _ { i } ^ { \prime \prime } = x _ { i } ^ { \prime }$ ($1 \leq i \leq 3$), i.e. the Hardy–Weinberg mapping is stationary. Conversely, if for $n = 3$ a stationary mapping $V$ is such that $p _ { 12,3 } = 1$ and all quadratic forms $x _ { j } ^ { \prime } \not\equiv 0$, then $V$ is a Hardy–Weinberg mapping (see [a2], [a6]). Thus, the only Mendelian heredity is stationary and such that all offsprings for the parental couple $( A A , a a )$ are $A a$ (and, in addition, such that all genotypes are present in the next generation).

For any stochastic quadratic mapping $V$, the linear form $f = \sum _ { j } a _ { j} x_j$ is called invariant if $f ^ { \prime } = f$. The mapping $V$ is called regular if there exists a family $\{ f _ { i } \} _ { 1 } ^ { m }$ of invariant linear forms such that $x _ { j } ^ { \prime } = \sum _ { i , k } c _ { i k } f _ { i } f _ { k }$ for certain constant coefficients $c _ { i k }$. The Hardy–Weinberg mapping is regular. Another interesting example is the quadrille mapping (see [a2], [a6]):

\begin{equation*} x _ { 1 } ^ { \prime } = p _ { 1 } q _ { 1 } ,\, x _ { 2 } ^ { \prime } = p _ { 1 } q _ { 2 }, \end{equation*}

\begin{equation*} x _ { 3 } ^ { \prime } = p _ { 2 } q _ { 1 } , x _ { 4 } ^ { \prime } = p _ { 2 } q _ { 2 } \end{equation*}

where

\begin{equation*} p _ { 1 } = x _ { 1 } + x _ { 2 } , \quad p _ { 2 } = x _ { 3 } + x _ { 4 }, \end{equation*}

\begin{equation*} q _ { 1 } = x _ { 1 } + x _ { 3 } , \quad q _ { 2 } = x _ { 2 } + x _ { 4 }. \end{equation*}

For the regular case, the Bernstein problem has been solved in [a6], [a8], [a9]. This is precisely the case when the stationarity is based on a system of genes (see [a7]). The genes correspond to the extremal rays of the cone of non-negative invariant linear forms. After a normalization, these forms are just the probabilities of the genes.

A standard genetical interpretation also requires $V$ to be normal in the sense that:

1) all $x _ { i } ^ { \prime } \neq 0$;

2) for any pair $i$, $k$ the quadratic forms $x _ { i } ^ { \prime }$ and $x _ { k} ^ { \prime }$ are not proportional;

3) there is no pair $i$, $k$ such that all $x _ { j } ^ { \prime }$ are functions of $x _ { i } + x _ { k }$ and of the remaining $x _ { l }$. If $V$ is normal together with its restrictions to all invariant faces of $\Delta ^ { n - 1 }$, then $V$ is called ultranormal. All stationary ultranormal mappings are regular [a13].

A non-regular stochastic Bernstein mapping appears already for $n = 3$:

\begin{equation*} x _ { 1 } ^ { \prime } = x _ { 1 } ( s + v ), \end{equation*}

\begin{equation*} x _ { 2 } ^ { \prime } = x _ { 3 } ^ { \prime } = \frac { 1 } { 2 } [ ( x _ { 1 } + x _ { 2 } ) s - x _ { 1 } v ], \end{equation*}

where $s = x _ { 1 } + x _ { 2 } + x _ { 3 }$, $v = x_3 - x_2$. Here, all invariant linear forms are trivial, i.e. proportional to $s$.

For the non-regular case there are some partial results for the Bernstein problem. In particular, the results cover the low dimensions $n = 4,5,6$ (cf. [a3], [a4], [a10]).

In the course of these investigations, Bernstein algebras were introduced as a powerful tool (see [a6], [a8] or Bernstein algebra). The theory of Bernstein algebras was subsequently developed by itself.

In the algebraic context, the Bernstein problem is to explicitly describe those Bernstein algebras which are stochastic with respect to the basis $\{ e _ { i } \} _ { 1 } ^ { n }$. The latter means that the product of every pair of basis vectors belongs to $\Delta ^ { n - 1 }$.

The Bernstein algebra corresponding to a regular mapping $V$ is regular by definition. This class is the most important from the genetics point of view.

Another important tool in the study of the Bernstein problem is a topological structure on the set of essential faces of $\Delta ^ { n - 1 }$, the faces such that their intersections with the image of $V$ are non-empty [a11].

See [a12] for a systematic presentation of the results and methods regarding the Bernstein problem up to the middle of the 1980s.

How to Cite This Entry:
Bernstein problem in mathematical genetics. Encyclopedia of Mathematics. URL: http://encyclopediaofmath.org/index.php?title=Bernstein_problem_in_mathematical_genetics&oldid=50460
This article was adapted from an original article by Yu.I. Lyubich (originator), which appeared in Encyclopedia of Mathematics - ISBN 1402006098. See original article