Mode and tempo of microsatellite length change in a malaria parasite mutation accumulation experiment

Microsatellite sequences are widely assumed to evolve neutrally, but also play an important role in bacterial pathogenesis, human disease and transcript abundance. The malaria parasite Plasmodium falciparum genome is extraordinarily AT-rich, containing 132,449 microsatellites-stretches of perfect 1-9 bp repeats between 10-1000bp, which comprise 10.74% of the 23 Mb genome. This project was designed to determine the mode and tempo of microsatellite mutations in malaria parasites. We maintained 31 parasite lines derived from a single 3D7 parasite cell for 114-267 days, with frequent bottlenecking to a single cell to minimize effective population size, allowing us to measure mutations accumulated over ~13,207 mitotic divisions. We Illumina sequenced the genomes of both progenitor and end-point mutation accumulation (MA) parasite lines in duplicate to validate stringent calling parameters. Calls were 99.89% (GATK), 99.99% (freeBayes) and 99.96% (HipSTR) concordant in duplicate sequence runs from independent sequence libraries. We observed 98 microsatellite mutations, giving rates of 2.11 × 10 -7 - 1.46 × 10 -8 /cell division that were strongly influenced by repeat motif and array length. Mutation rate was low relative to other organisms. However, despite this, in a single infection (1011 parasites) there will be 1.46 × 10 3 - 2.11 × 10 4 independent mutations at any single microsatellite locus. Given that many microsatellites are found in promotors, introns, within or close to coding sequences, we suggest that they may be important regulators of transcriptional and phenotypic variation in this pathogen.