Related Experiment Video
Updated: May 10, 2026
06:55
Synthesis of Cyclic Polymers and Characterization of Their Diffusive Motion in the Melt State at the Single Molecule Level
Published on: September 26, 2016
Complete numerical solution of the diffusion equation of random genetic drift
Lei Zhao1, Xingye Yue, David Waxman
1Centre for Computational Systems Biology, Fudan University, Shanghai 20433, People's Republic of China.
Genetics
|June 11, 2013
Summary
A new numerical method accurately models random genetic drift by solving the diffusion equation, ensuring probability conservation for complete solutions. This approach unifies calculations for allele frequencies, fixation, and loss, even with changing population sizes.
Area of Science:
- Population Genetics
- Computational Biology
- Mathematical Biology
Background:
- Random genetic drift is a key evolutionary mechanism influencing allele frequencies.
- Existing methods for solving the diffusion equation for genetic drift have limitations in probability conservation and handling boundary conditions.
- A unified computational framework is needed for comprehensive drift analysis.
Purpose of the Study:
- To present a novel numerical method for solving the diffusion equation of genetic drift at a single locus with two alleles.
- To develop a probability-conserving method that yields complete solutions, including fixation and loss probabilities.
- To create a versatile framework applicable to various genetic drift scenarios, including time-dependent parameters.
Main Methods:
- Developed a numerical method to solve the diffusion equation for allele frequency changes due to random genetic drift.
- Ensured the numerical solution conserves probability, resulting in a complete probability distribution (total probability of unity).
- Implemented the method to handle both internal allele frequencies and boundary conditions (fixation and loss).
Main Results:
- The numerical method produces complete solutions that inherently include fixation and loss probabilities.
- Demonstrated the method's ability to accurately model allele frequencies and fixation probabilities under neutrality, selection, and demographic changes.
- Validated the framework's straightforward implementation for diverse genetic drift calculations.
Conclusions:
- The presented numerical method offers a unified and robust framework for studying random genetic drift.
- The probability-conserving approach ensures accurate representation of allele frequency dynamics, including fixation and loss.
- This method facilitates broader applications in population genetics, especially for scenarios with time-varying parameters.
Related Concept Videos
03:33Genetic Drift
Natural selection—probably the most well-known evolutionary mechanism—increases the prevalence of traits that enhance survival and reproduction. However, evolution does not merely propagate favorable traits, nor does it always benefit populations.
01:09Mutation, Gene Flow, and Genetic Drift
In a population that is not at Hardy-Weinberg equilibrium, the frequency of alleles changes over time. Therefore, any deviations from the five conditions of Hardy-Weinberg equilibrium can alter the genetic variation of a given population. Conditions that change the genetic variability of a population include mutations, natural selection, non-random mating, gene flow, and genetic drift (small population size).
01:49Hardy-Weinberg Principle
Diploid organisms have two alleles of each gene, one from each parent, in their somatic cells. Therefore, each individual contributes two alleles to the gene pool of the population. The gene pool of a population is the sum of every allele of all genes within that population and has some degree of variation. Genetic variation is typically expressed as a relative frequency, which is the percentage of the total population that has a given allele, genotype or phenotype.
01:25Modeling with Differential Equations
Population dynamics can be described mathematically by considering the population size P(t) as a function of time. The rate of change of the population is then represented by the derivative of P(t). A simple assumption is that the rate of growth is proportional to the size of the population itself. This leads to an exponential growth model, where the population increases rapidly without bound. While this is a useful first approximation, it does not reflect realistic long-term...
01:29Mechanistic Models: Compartment Models in Algorithms for Numerical Problem Solving
Mechanistic models play a crucial role in algorithms for numerical problem-solving, particularly in nonlinear mixed effects modeling (NMEM). These models aim to minimize specific objective functions by evaluating various parameter estimates, leading to the development of systematic algorithms. In some cases, linearization techniques approximate the model using linear equations.
In individual population analyses, different algorithms are employed, such as Cauchy's method, which uses a...
In individual population analyses, different algorithms are employed, such as Cauchy's method, which uses a...
01:27Linear Differential Equations
The integrating factor method provides a systematic way to solve first-order linear differential equations, especially those that cannot be handled by separation of variables. This method is particularly useful in modeling time-dependent physical systems influenced by both constant inputs and resistive forces. A common example is the motion of a car subjected to a constant engine force while experiencing air resistance proportional to its velocity.In such scenarios, Newton’s second law yields a...