Abstract:We estimate density of defects frozen into a biological Turing pattern which was turned on at a finite rate. A self-locking of gene expression in individual cells, which makes the Turing transition discontinuous, stabilizes the pattern together with its defects. A defect-free pattern can be obtained by spatially inhomogeneous activation of the genes.
Biological Physics,Condensed Matter,Pattern Formation and Solitons,Chemical Physics,Quantitative Biology
What problem does this paper attempt to address?
The problem that this paper attempts to solve is the density estimation of defects being frozen into permanent patterns during the formation of dynamic Turing patterns. Specifically, the paper explores how defects form in patterns when Turing patterns are initiated at a finite rate, and these defects become stable due to the self - locking mechanism of gene expression. In addition, the paper also proposes the possibility of obtaining a perfect defect - free uniform pattern by spatially non - uniformly activating genes.
### Background and Research Motivation
- **Turing Patterns**: Turing pointed out that simple reaction - diffusion systems can explain the formation of biological patterns. However, Turing patterns in continuous space are usually smooth and non - permanent, while actual biological patterns are often permanent, and the functions of adjacent cells and the expressed genes may be significantly different.
- **Defect Problem**: In Turing patterns, the existence of defects is an important problem. For some biological structures, such as the formation of organs, a single defect may be fatal. Therefore, it is crucial to understand the origin of defects and control methods.
### Research Methods
- **Model Description**: The paper uses a simplified toy model, which contains two genes \(A\) and \(B\), which are strongly mutually inhibitory. Gene expression is controlled by a common activator \(a\). When the concentration of \(a\) exceeds a critical value \(a_c\), genes \(A\) and \(B\) are simultaneously activated in cells.
- **Dynamical Equations**: The dynamics of the model are described by stochastic differential equations, taking into account protein degradation and the switching process of gene expression. In the adiabatic limit, it can be simplified to mean - field equations.
- **Phase - transition Analysis**: The paper analyzes the attractor structure of the system and finds that there are two critical activation levels \(a_{c1}\) and \(a_{c2}\), and the system exhibits different phases and symmetry breaking between these two critical points.
### Main Results
- **Defect Density Estimation**: When the activator \(a\) increases at a finite rate, the system will form finite correlation domains, which independently select the final phase (\(\alpha\beta\) or \(\beta\alpha\)). Due to the critical slowing - down effect, a certain finite correlation length \(\hat{\xi}\) will be "frozen" into fluctuations, which determines the density of defects in Turing patterns.
- **Kibble - Zurek Scenario**: This effect is similar to the Kibble - Zurek scenario in cosmology and condensed - matter physics, where the density of defects has a power - law relationship \(\hat{\xi}\sim\tau^{1/4}\) with the transition rate \(\tau\).
- **Defect - free Pattern Generation**: The paper also proposes a method of generating defect - free uniform patterns by spatially non - uniformly activating genes. By activating \(a\) at the initial point and then gradually diffusing, the initially selected phase (such as \(\alpha\beta\)) can gradually expand to the entire system through the coupling between cells, thereby avoiding the formation of defects.
### Conclusion
Through theoretical analysis and numerical simulation, the paper reveals the mechanism of defect formation in Turing patterns and proposes an effective method for generating defect - free biological patterns. These results are of great significance for understanding pattern formation in biological development processes.