Coalescent simulation in continuous space: algorithms for large neighbourhood size.

Area of Science:

• Population genetics
• Spatial modeling
• Evolutionary biology

Background:

• Continuous spatial distributions are common in nature, but traditional population models use artificial grids (demes).
• Existing models like isolation by distance have technical flaws, hindering accurate simulation.
• A new 2D extinction-recolonization model offers a robust framework for spatially continuous populations.

Purpose of the Study:

• To develop efficient and exact simulation algorithms for the 2D coalescent process.
• To address the inefficiency of direct simulation for large neighborhood sizes in spatial population models.

Main Methods:

• Development of novel algorithms for simulating the coalescent process in a spatial continuum.
• Analysis of the efficiency and accuracy of these new simulation algorithms.
• Application to arbitrary sample sizes and multiple genetic loci.

Main Results:

• Efficient and exact algorithms for simulating spatial coalescent processes were successfully developed.
• The new algorithms overcome the computational limitations of previous methods for large neighborhood sizes.
• The simulation approach is applicable to various sample sizes and numbers of loci.

Conclusions:

• The developed algorithms provide a computationally tractable method for studying populations with continuous distributions.
• This work establishes a rigorous foundation for analyzing spatial population genetics using a 2D extinction-recolonization model.
• The findings facilitate more accurate modeling of evolutionary processes in geographically widespread species.