PhyloSim

PhyloSim simulates sequence evolution using Monte Carlo methods to model substitutions, insertions, deletions, rate variation, and selective constraints for molecular evolutionary studies.

Key Features:

• Extensibility and Object-Oriented Framework: Implements an extensible object-oriented framework in R to enable incorporation of additional model complexity into simulations.
• Gillespie Algorithm Integration: Employs the Gillespie algorithm to simulate multiple concurrent processes such as substitutions, insertions, and deletions.
• Complex Rate Variation and Indel Processes: Simulates arbitrarily complex patterns of rate variation and multiple indel (insertion–deletion) processes.
• Selective Constraints on Indels: Imposes selective constraints on indel events to model the influence of selection on insertion–deletion dynamics.
• User-Defined Mutation and Selection Patterns: Supports user-defined mutation and selection patterns for custom simulation scenarios.
• Integration with R: Implemented in R for direct integration with R-based statistical and graphical tools.

Scientific Applications:

• Phylogenetic inference benchmarking: Generates simulated sequence data to assess the performance of phylogenetic inference methods.
• Sequence alignment evaluation: Produces alignments under controlled evolutionary scenarios to evaluate sequence alignment algorithms.
• Modeling evolutionary scenarios: Enables exploration of realistic evolutionary dynamics, including complex rate variation and indel behavior, for studies in molecular evolution and phylogenetics.

Methodology:

Monte Carlo simulation of sequence evolution implemented in R using an extensible object-oriented framework; uses the Gillespie algorithm to manage concurrent substitutions, insertions, and deletions; simulates complex rate variation, multiple indel processes, selective constraints on indels, and user-defined mutation and selection patterns.