Analysis of a mechanistic Markov model for gene duplicates evolving under subfunctionalization

Area of Science:

• Evolutionary biology
• Genomics
• Mathematical modeling

Background:

• Gene duplication is a primary driver of functional innovation in genomes.
• Subfunctionalization, where gene duplicates divide ancestral functions, is a key evolutionary model.
• Previous analyses of subfunctionalization were often coarse-grained and lacked mechanistic depth.

Purpose of the Study:

• To develop and analyze a mechanistic mathematical model for gene duplicate evolution via subfunctionalization.
• To derive testable predictions and estimate biologically relevant parameters from genomic data.
• To compare model predictions against empirical data from multiple species.

Main Methods:

• Developed a mathematical model based on subfunctionalization mechanics and Poisson mutation rates.
• Utilized Phase-Type distribution theory for exact analytical solutions.
• Fitted the model's survival function to genomic data from *Homo sapiens*, *Mus musculus*, *Rattus norvegicus*, and *Canis familiaris*.

Main Results:

• The mechanistic model provides testable predictions and allows estimation of parameters like regulatory region counts and mutation rates.
• Model fits to genome data suggest duplicates typically possess few regulatory regions.
• Estimated mutation rates indicate coding regions evolve 5-10 times faster than regulatory regions in duplicates.

Conclusions:

• The subfunctionalization model shows strong agreement with empirical genomic data.
• This mechanistic approach offers a consistent explanation for the evolution of numerous gene duplicates.
• Provides the first model-based estimates for the number of regulatory regions in gene duplicates.