Selection constrains high rates of satellite DNA mutation in Daphnia pulex

A long-standing evolutionary puzzle is that all eukaryotic genomes contain large amounts of tandemly-repeated satellite DNA whose composition varies greatly among even closely related species. To elucidate the evolutionary forces governing satellite dynamics, quantification of the rates and patterns of mutations in satellite DNA copy number and tests of its selective neutrality are necessary. Here we used whole-genome sequences of 28 mutation accumulation (MA) lines of Daphnia pulex in addition to six isolates from a non-MA population originating from the same progenitor to both estimate mutation rates of abundances of satellite sequences and evaluate the selective regime acting upon them. We found that mutation rates of individual satellite sequence \"kmers\" were both high and highly variable, ranging from additions/deletions of 0.29 - 105 copies per generation (reflecting changes of 0.12 - 0.80 percent per generation). Our results also provide evidence that new kmer sequences are often formed from existing ones. The non-MA population isolates showed a signal of either purifying or stabilizing selection, with 33 % lower variation in kmer abundance on average than the MA lines, although the level of selective constraint was not evenly distributed across all kmers. The changes between many pairs of kmers were correlated, and the pattern of correlations was significantly different between the MA lines and the non-MA population. Our study demonstrates that kmer sequences can experience extremely rapid evolution in abundance, which can lead to high levels of divergence in genome-wide satellite DNA composition between closely related species.