Figure 3
The effect of clonal decay on Tajima's D, Fu and Li's D*. Individual based simulation (upper panel) assumed a population of 500 individuals fed by new clones at rate c, deleterious mutations accumulating at rate 0.5, each with a selection coefficient s. Neutral model assumes s = 0. To demonstrate the effect of varying U, the black square and black dot symbols indicate the means values for U = 0.1 and U = 0.05, respectively. Selection coefficient s = 0.5 in both cases. Coalescent simulation (lower panel) assumed two sexual demes of a total size of Nc = 25,000 individuals interconnected by migration rate mc and asexual population of Nc = 25,000 individuals composed of 10 clonal lineages. Each clone had a constant probability Pc of extinction per generation or it increased at rate (1 - (1 - f')α'). We sampled 5 individuals per deme/clone. Mutations in the neutral locus accumulated at rate 10-8 in both cases. In the coalescent simulation, data are estimated according to mutation overlaid along the total pedigree (curves called 'Total') or along the branches evolved during the asexual phase of sequence evolution (curves called 'Pruned'). Black colour below or above the grey parts of bars indicates the proportions of cases where simulated values of D, D* were lower or higher than the 95% CI of the neutral expectation.