Retracing and rewriting the evolutionary trajectories of mammalian developmental enhancers

) with massively parallel reporter assays (MPRAs). This phylogenetic dissection revealed diverse modes of evolution, from lineage-restricted activity to deep functional conservation despite extensive sequence divergence. To pinpoint causal changes, we developed a model-driven reconstitution strategy that uses deep learning-based predictions of chromatin accessibility to re-introduce a succession of mutations into ancestral orthologs; this revealed critical transcription factor binding site (TFBS) changes and pervasive context-dependent epistasis, including instances where mutational effects were strongly contingent on the order of their introduction. When we extended this strategy to tune the activity of extant orthologs, we found that ablation of enhancer function required as few as one to seven mutations, whereas enhancement was constrained by element-specific activity ceilings-a striking asymmetry in the predictability of model-guided enhancer editing. Together, these results shed light on how the plasticity of mammalian enhancers intersects with their evolution, and advance a framework for reprogramming the activity of endogenous CREs at nucleotide resolution.