Tracking somatic drift reveals small effective population size of stem cells and high somatic mutation rate in asexual planaria

Planarian flatworms have emerged as highly promising models of body regeneration due to the many stem cells scattered through their bodies. Currently, there is no consensus as to the number of stem cells active in each cycle of regeneration or the equality of their relative contributions. We approached this problem with a population genetic model of somatic genetic drift. We modeled the fissiparous life cycle of asexual planarians as an asexual population of cells that goes through repeated events of splitting into two subpopulations followed by population growth to restore the original size. We sampled a pedigree of obligate asexual clones of Girardia cf. tigrina at multiple time points encompassing 14 generations. Effective population size of stem cells was inferred from the magnitude of temporal fluctuations in the frequency of somatic variants and under most of the examined scenarios was estimated to be in the range of a few hundreds. Average genomic nucleotide diversity was 0.00398. Assuming neutral evolution and mutation-drift equilibrium, the somatic mutation rate was estimated in the 10-5 - 10-7 range. Alternatively, we estimated Ne and somatic from temporal changes in nucleotide diversity{pi} without the assumption of equilibrium. This second method suggested even smaller Ne and larger . A key unknown parameter in our model on which estimates of Ne and depend is g, the ratio of cellular to organismal generations determined by tissue turnover rate. Small effective number of propagating stem cells might contribute to reducing reproductive conflicts in clonal organisms.