A population genetic model to infer allotetraploid speciation and long-term evolution applied to two yarrow species

Allotetraploid speciation, that is, the generation of a hybrid tetraploid species from two diploid species, and the long-term evolution of tetraploid populations and species are important in plants. We developed a population genetic model to infer population genetic parameters of tetraploid populations from data of the progenitor and descendant species. Two yarrow species, Achillea alpina-4x and A. wilsoniana-4x, arose by allotetraploidization from the diploid progenitors, A. acuminata-2x and A. asiatica-2x. Yet, the population genetic process has not been studied in detail. We applied the model to sequences of three nuclear genes in populations of the four yarrow species and compared their pattern of variability with that in four plastid regions. The plastid data indicated that the two tetraploid species probably originated from multiple independent allopolyploidization events and have accumulated many mutations since. With the nuclear data, we found a low rate of homeologous recombination or gene conversion and a reduction in diversity relative to the level of both diploid species combined. The present analysis with a novel probabilistic model suggests a genetic bottleneck during tetraploid speciation, that the two tetraploid species have a long evolutionary history, and that they have a small amount of genetic exchange between the homeologous genomes.