Neutral evolution of duplicated DNA: an evolutionary stick-breaking process causes scale-invariant behavior

Area of Science:

• Genomics
• Evolutionary Biology
• Bioinformatics

Background:

• Identical matching sequences are increasingly found in eukaryotic genomes.
• The length distribution of these sequences follows a power law, prompting investigation into evolutionary mechanisms or functional constraints.
• Understanding the origin of these sequence patterns is crucial for genome evolution studies.

Purpose of the Study:

• To introduce a simple, evolutionarily neutral model for the formation of identical matching sequences.
• To explain the observed power law distribution in sequence lengths.
• To determine if the power law exponent is universal.

Main Methods:

• Development of a neutral evolutionary model incorporating point mutations and segmental duplications.
• Analysis of statistical features generated by the model.
• Extension of a mathematical model for random stick breaking.

Main Results:

• The proposed model successfully replicates the statistical features of genomic data, including the power law distribution.
• Analytical results show the exponent of the power law tail is -3.
• The exponent is demonstrated to be universal, independent of model specifics.

Conclusions:

• A simple model of point mutations and segmental duplications can explain the observed enrichment and length distribution of identical matching sequences in eukaryotic genomes.
• The power law exponent of -3 is a universal feature, suggesting fundamental evolutionary principles at play.
• This finding provides insights into the evolutionary dynamics shaping genome sequences.