Abstract:Consider an arbitrary large population at the present time, originated at an unspecified arbitrary large time in the past, where individuals in the same generation reproduce independently, forward in time, with the same offspring distribution but potentially changing among generations. In other words, the reproduction is driven by a Galton-Watson process in a varying environment. The genealogy of the current generation backwards in time is uniquely determined by the coalescent point process $(A_i, i\geq 1)$, where $A_i$ is the coalescent time between individuals $i$ and $i+1$. In general, this process is not Markov. In constant environment, Lambert and Popovic (2013) proposed a Markov process of point measures to reconstruct the coalescent point process. We present a counterexample where we show that their process does not have the Markov property. The main contribution of this work is to propose a vector valued Markov process $(B_i,i\geq 1)$, that reach the goal to reconstruct the genealogy, with finite information for every $i$. Additionally, when the offspring distributions are lineal fractional, we show that the variables $(A_i, i\geq 1)$ are independent and identically distributed.
What problem does this paper attempt to address?
### The problems the paper attempts to solve
This paper attempts to solve the problem of reconstructing the population genealogy driven by the Galton - Watson process in a changing environment. Specifically, the author focuses on how to reconstruct the genealogy of individuals in the current generation with limited information in a changing environment.
### Background and motivation
In genetics and evolutionary biology, understanding the genealogy of a population is very important. Traditional research methods usually assume that the environment is constant, but in reality, the environment is often changing. Therefore, a new method is needed to deal with the problem of genealogy reconstruction in such a changing environment.
### Main contributions
1. **Propose a new Markov process**:
- The author proposes a new vector - valued Markov process \((B_i, i\geq0)\), which can use limited information to reconstruct the genealogy of the population.
- This process has good properties in a changing environment, especially in the case of linear - fraction distribution, \((A_i, i\geq1)\) are independent and identically distributed random variables.
2. **Correct predecessors' mistakes**:
- The author points out that the process \((\tilde{B}_i, i\geq0)\) proposed by Lambert and Popovic (2013) lacks the Markov property and provides a counterexample.
- By introducing the new process \((B_i, i\geq0)\), the author successfully solves this problem.
### Technical details
- **Galton - Watson process**:
- In a changing environment, the Galton - Watson process describes the development of the number of individuals, where each individual reproduces independently, sharing the same offspring distribution, but these distributions may be different between different generations.
- The author defines the generating functions \(f_n\) and \(f_{m,n}\) and discusses their related properties.
- **Genealogy and the coalescent point process**:
- The author defines the coalescent point process \((A_i, i\geq1)\), where \(A_i\) represents the coalescent time between the \(i\)-th individual and the \((i + 1)\)-th individual.
- By introducing the auxiliary Markov process \((D_i, i\geq0)\), the author shows how to extract the information of \(A_i\) from \(D_i\).
- **The new Markov process \((B_i, i\geq0)\)**:
- The author defines \(B_i\) by restricting the first \(C_{1,i + 1}\) entries of \(D_i\), where \(C_{1,i+1}\) represents the number of generations when the first individual on the \(i\)-th ridge becomes an ancestor of the individuals in the current generation.
- The author proves that \((B_i, i\geq0)\) is a Markov chain and gives its transition probabilities.
### Conclusion
This paper successfully solves the problem of genealogy reconstruction in a changing environment by proposing a new vector - valued Markov process \((B_i, i\geq0)\). This method is not only of great theoretical significance but also of high value in practical applications, especially in the field of evolutionary biology.